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ABSTRACT 


The  flowfield  about  a  Vertically-Launched  Surface-to-Air  Missile  model 
at  an  angle  of  attack  of  50°  and  a  Reynolds  number  of  1.1x10s  was  investi¬ 
gated  in  a  low  speed  wind  tunnel  at  the  Naval  Postgraduate  School.  The  goal 
of  this  thesis  is  to  determine  the  location  and  intensity  of  the  asymmetric 
vortices  in  the  wake  of  the  VLSAM  model  and  to  display  these  vortices  by 
velocity  mapping  and  pressure  contours.  The  two  model  configurations  tested 
were  for  a  cruciform  missile  with  wings  and  tails;  one  at  0°  rc.ll  angle  ("plus" 
aspect)  and  the  other  at  a  45°  roll  angle  ("cross"  aspect).  Two  flowfield 
conditions  were  treated:  the  nominal  ambient  wind  tunnel  condition  and  a 
condition  with  a  grid-generated  turbulence  of  length  scale  1.08  inches  and 
1.88%  turbulence  intensity.  The  turbulence  length  scale  is  61.7%  of  the  model 
diameter  and  4.7%  of  the  model  length.  The  following  conclusions  were 
reached:  1)  An  increase  in  turbulence  intensity  tended  to  reduce  the  strength 
of  the  asymmetric  nose-generated  vortices;  2)  the  two  asymmetric  vortices 
remained  in  approximately  the  same  position  for  an  increase  in  turbulence;  3) 

"cross"  aspect  vortices  were  more  diffused,  slightly  larger  and  centered 
further  away  from  the  model  surface  than  those  of  the  "plus"  aspect  body 
configuration,  which  correlates  with  the  differences  in  induced  side  forces  for 
these  configurations  observed  \  Rabang;  4)  the  top  vortex  of  the  two 
asymmetric  vortices  was  closer  .  .he  model  surface  and  appeared  to  be 
stronger  for  both  configurations;  and  5)  the  addition  of  wings  and  tails  did  not 

greatly  alter  the  vortex  pattern  around  the  nose  of  the  missile  model.  □ 

d  □ 

I  Jojuc  icoU on _ _ 


By -  1 

Dlstr 

Aval 

i but  Ion/ 

lability  Cods a 

Diet 

Avail 

Spec 

mmmmi 

and/or 

lal 

TABLE  OF  CONTENTS 


I.  INTRODUCTION  .  1 

A.  BACKGROUND .  1 

B.  HIGH  ANGLE  OF  ATTACK  AERODYNAMICS . 3 

1.  Aerodynamic  Regimes . 3 

2.  Asymmetric  Vortex  Theory  . 5 

3.  Two-Dimensional  Crossflow . 6 

4.  Three-Dimensional  Vortices . 11 

C.  TURBULENCE . 13 

D.  LIFTING  SURFACE  EFFECTS  . 14 

E.  VLSAM  LAUNCH  ENVIRONMENT . 16 

1.  Marine  Environment . 16 

2.  Launch  and  Cross  wind  Velocities . 17 

3.  Additional  Launch  Considerations . 18 

II.  EXPERIMENT  AND  PROCEDURES . 19 

A.  PURPOSE . 19 

B.  APPARATUS . 21 

1.  Wind  Tunnel . . 21 

2.  Turbulence- Generating  Grids . 23 

3.  VLSAM  Model  and  Support  Equipment . 29 

4.  Velmex  8300  3-D  Traverser . 29 

5.  5-Hole  Pressure  Probe . 32 

6.  Scanivalve  and  HP  Data  Acquisition  System  . 32 


iv 


C.  EXPERIMENTAL  CONDITIONS . 34 

D.  SOFTWARE  AND  PROCEDURES . 37 

1.  PPROBE  Program . 38 

2.  CALP  Program . 40 

3.  CONVERT  Program . 42 

4.  TECPLOT  System . 43 

E.  PRELIMINARY  TESTS  . 43 

1 .  Dynamic  Pressure  Calibration . 43 

2.  Yaw  and  Pitch  Angle  Corrections . 44 

III.  RESULTS . - . 45 

i 

A.  CONFIGURATION  0A  ("PLUS"  WITHOUT  TURBULENCE) . 45 

B.  CONFIGURATION  3A  ("PLUS  '  WITH  TURBULENCE  GRID  3) .  48 

C.  CONFIGURATION  0C  ("CROSS"  WITHOUT  TURBULENCE) . 52 

D.  CONFIGURATION  3C  ("CROSS"  WITH  TURBULENCE  GRID  3) .  57 

E.  COMPARISONS . - . 57 

1.  Between  Body  Configurations  (A  and  C)  .  57 

2.  Between  Turbulence  Levels  (0  and  3)  .  61 

3.  With  Body-Only  Configuration  (B)  . 62 

IV.  CONCLUSIONS  AND  RECOMMENDATIONS . 66 

APPENDIX  A.  PPROBE  Program . 68 

APPENDIX  B.  CALP  Program  . 79 

APPENDIX  C.  CONVERT  Program  . 84 

APPENDIX  D.  RESULT  00.DAT  . 87 

APPENDIX  E.  RESULT  OA.DAT . - . 90 

APPENDIX  F.  RESULT  3A.DAT  — - . 100 


v 


APPENDIX  G.  RESULT  OC.DAT . 110 

APPENDIX  H.  RESULT  3C.DAT . 120 

LIST  OF  REFERENCES . 130 

INITIAL  DISTRIBUTION  LIST . 133 


LIST  OF  FIGURES 


Figure 

1 .  Airflow  Regimes . 4 

2.  Vortex  Flow  About  a  Slender  Nose  Cylinder . 6 

3.  Two-Dimensional  Crossflow  About  a  Cylinder . 7 

4.  Side  Force  to  Normal  Force  Ratio . 9 

5.  Effect  of  Reynold's  Number  on  Maximum  Side  Force  at  «==55° .  10 

6.  Side  Force  Variations  with  Nose  Roll  Angle  . 11 

7.  The  Planar  Survey  Grid . 20 

8.  Naval  Postgraduate  School  Wind  Tunnel . 22 

9.  Planview  of  VLSAM  Model  With  Pressure  Probe  and 

Grid  In  the  Test  Section  of  the  Wind  Tunnel  . 24 

10.  Square-mesh  Turbulence-Generating  Grid . 25 

11.  Grid  Turbulence  Intensities . 26 

12.  Grid  Turbulence  Length  Scales . 27 

13.  Turbulence-Generating  Grids  . 28 

14.  Specifications  of  VLSAM  Model . 30 

15.  Velmex  8300  Traversing  Assembly . 31 

16.  5-Hole  Pressure  Probe  and  Measuring  Tip . 33 

17.  HP  Data  Acquisition  System . 35 

18.  VLSAM  Model  Body  Configurations . 36 

19.  Program/Data  File  Flowchart  . 37 

20.  Calibration  Manometer . 41 

21.  Velocity  Vector  Plot  -  Configuration  OA . 46 

vii 


22.  Total  Pressure  Coefficient  -  Configuration  OA 

23.  Static  Pressure  Coefficient  -  Configuration  OA 

24.  Velocity  Vector  Plot  -  Configuration  3A . 

25.  Total  Pressure  Coefficient  -  Configuration  3A 

26.  Static  Pressure  Coefficient  -  Configuration  3A 

27.  Velocity  Vector  Plot  -  Configuration  OC - 

28.  Total  Pressure  Coefficient  -  Configuration  OC 

29.  Static  Pressure  Coefficient  -  Configuration  OC 

30.  Velocity  Vector  Plot  -  Configuration  3C - 

31.  Total  Pressure  Coefficient  -  Configuration  3C 

32.  Static  Pressure  Coefficient  -  Configuration  3C 

33.  Velocity  Vector  Plot  -  Configuration  OB . 

34.  Total  Pressure  Coefficient  -  Configuration  OB 

35.  Static  Pressure  Coefficient  -  Configuration  OB 


LIST  OF  TABLES 


Table 

1.  Grid  Specifications . 25 

2.  Grid  Turbulence  Parameters  . 27 


NOMENCLATURE 


Cps  Static  Pressure  Coefficient 

Cpx  Total  Pressure  Coefficient 

d  Base  diameter  of  missile  body 

K  Wind  tunnel  calibration  factor 

In  Nose  length 

Ld  Missile  diameter  scale 

Lj  Missile  length  scale 

Lu  Longitudinal  turbulence  length  scale 

Q  Freestream  dynamic  pressure 

Ps  Freestream  static  pressure 

Pt  Freestream  total  pressure 

PSL  Local  static  pressure 

PtL  Local  total  pressure 

Pi  Total  pressure 

P2  Lateral  static  pressure 

P3  Lateral  static  pressure 

P4  Pitch  angle  pressure 

P5  Pitch  angle  pressaure 

Re  Reynold's  number 

Tu  Turbulence  intensity 

u'  Root-mean- square  (rms)  velocity  fluctuation 

Um  Measured  velocity 


x 


U°o  Longitudinal  mean  velocity 

Z  Roughness  length 

«  Angle  of  attack  (AOA) 

«av  AOA  for  the  formation  of  asymmetric  vortices 

e  Blockage  correction  factor 

0a  Nose  semi-vertex  angle 

<}>  Angle  from  crossflow 

0R  Roll  angle 

p  Air  density 


xi 


ACKNOWLEDGEMENTS 


Sincere  graditude  goes  to  my  thesis  advisor,  Dr.  Richard  M.  Howard,  for 
his  guidance  and  encouragement  throughout  the  course  of  my  study.  I  have 
been  provided  with  a  rewarding  learning  experience  by  being  afforded  the 
chance  to  conduct  wind  tunnel  experiments. 

I  would  also  like  to  thank  the  following  Naval  Postgraduate  School  per¬ 
sonnel  for  their  support  in  this  research: 

Mr.  Jack  King,  Aerolab 
Mr.  John  Moulton,  Aeronautics  Shop 
Capt.  Dan  Johnson,  USA,  NPS  Student 
Lt.  Jim  Pinaire,  USN,  NPS  Student. 


I.  INTRODUCTION 


A.  BACKGROUND 

The  development  of  the  Vertical  Launch  Surface -to- Air  Missile  (VLSAM) 
has  provided  greater  weapon  system  reliability,  availability  and  flexibility 
over  its  predecessors.  Earlier  systems,  such  as  trainable  launchers  and  box 
launchers,  were  rather  cumbersome  and  required  excessive  deck  space,  thus 
limiting  the  number  of  missiles  per  ship.  Additionally,  these  systems  were 
slow  to  reload  since  trainable  launchers  required  cycling  of  rounds  to  get  the 
desired  one  in  position  for  launch,  while  box  launchers  were  reloaded  by 
hand.  The  vertical  launcher  provides  for  quick  access  to  any  round  without 
cycling,  rapid  reloading  and,  due  to  its  design  and  the  fact  that  the  missile 
blast  is  kept  in  a  concentrated  area  near  the  launcher,  it  requires  less  deck 
space.  [Refs.  1,  2  and  3] 

The  VLSAM’s  trajectory  allows  it  to  point  to  its  target  after  launch  ami 
subsequently  guide  itself  to  the  correct  heading.  Its  aerodynamic  characteris¬ 
tics  may  significantly  change  as  it  operates  from  subsonic  to  supersonic 
speeds  during  the  launch,  midcourse  and  terminal  phases  of  its  flight.  [Ref.  2] 
The  varying  flight  control  requirements  and  the  airflow  about  the  missile 
during  these  phases  provide  for  different  flight  regimes. 

When  launched,  the  VLSAM  enters  an  open  ocean  environment  and  is 
subject  to  potentially  significant  cross-winds,  the  result  of  which  is  a  missile 
flying  at  relatively  low  velocities  at  a  high  angle  of  attack.  [Refs.  1  and  2]  In 
particular,  an  example  of  this  low  velocity/high  angle  of  attack  condition  is 
the  Standard  Missile  2-Block  4,  or  AEGIS  Extended  Range  version,  which 
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has  an  added  sustainer  section  and  leaves  its  launcher  at  a  much  lower  speed 
than  the  unboosted  SM-2  version.  A  missile  flying  at  relatively  low  velocities 
at  a  high  angle  of  attack  may  bring  with  it  the  formation  of  asymmetric 
vortices  about  the  missile  nose  and  afterbody  which  can  generate  unpre¬ 
dictable  and  undesired  side  forces,  thereby  affecting  overall  missile  flight  sta¬ 
bility.  The  airflow  around  a  ship’s  superstructure,  the  ocean  surface  condi¬ 
tions  and  the  atmospheric  boundary  layer  conditions  during  launch  may  also 
provide  turbulent  flow  over  the  missile.  [Ref.  2]  The  advent  of  the  VLSAM 
and  the  desire  to  have  highly  maneuverable  supersonic  missiles  have 
increased  the  need  for  further  studies  in  high  angle  of  attack  research  under 
various  flowfield  conditions. 

This  thesis  continues  experiments  and  research  conducted  to  date  at  the 
Naval  Postgraduate  School  (NPS)  to  understand  the  effects  of  turbulence  on 
the  VLSAM  aerodynamic  characteristics.  Previously,  Roane  [Ref.  1]  developed 
a  system  model  to  measure  flowfield  turbulence  through  the  use  of  four  dif¬ 
ferent  grids  which  generate  varying  turbulence  levels  in  the  NPS  low  speed 
wind  tunnel.  Rabang  [Ref.  2]  studied  the  effects  of  this  turbulence  on  the 
asymmetric  vortex  forces  on  the  missile  model.  Lung  [Ref.  3]  determined  the 
location  and  intensity  of  the  asymmetric  vortices  in  the  wake  of  the  model  by 
experimental  flowfield  measurements  about  a  body-only  missile  configuration 
with  and  without  freestream  turbulence.  Similarly,  the  goal  of  this  study  is  to 
determine  vortex  locations  and  intensities,  through  a  series  of  wind  tunnel 
experiments,  for  two  missile  configurations,  both  with  and  without  turbu¬ 
lence.  The  configurations  considered  are  for  a  cruciform  missile  with  wings 
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and  tails;  one  at  0°  roll  angle  ("plus"  aspect)  and  the  other  at  a  45°  role  angle 
("cross"  aspect). 

B.  HIGH  ANGLE  OF  ATTACK  AERODYNAMICS 

In  high  angle  of  attack  aerodynamics,  flow  separation  from  the  body, 
wing  and  tail  surfaces  is  important  due  to  the  strong  normal  and  side  forces 
which  may  be  generated.  Major  factors  which  influence  flow  separation  are 
nose  shape,  angle  of  attack,  crossflow  Reynolds  number  and  nose  fineness 
ratio.  Other  factors  include  roll  angle  and  roll  rate,  free  stream  turbulence, 
surface  roughness,  acoustic  environment  and  VLSAM  model  vibrations 
[Ref.  4]. 

1.  Aerodynamic  Regimes 

As  the  angle  of  attack,  «,  of  a  slender  body  of  revolution  ranges 
from  0°  to  90°,  there  are  at  least  four  distinct  airflow  regimes  through  which 
the  missile  body  transitions.  [Ref.  4] 

(1)  Regime  I  (0°<  <  5°):  At  very  low  angles  of  attack,  the  axial  flow  domi¬ 

nates  and  there  is  no  discemable  boundary  layer  separation  (flow  is 
attached). 

(2)  Regime  II  (5°<  «  <  20°);  At  intermediate  angles  of  attack,  boundary 
layer  separation  occurs  on  the  lee  side  of  the  body.  This  becomes  a  free 
shear  layer  which  rolls  up  into  a  symmetric  vortex  pair  that  is  steady 
with  time.  No  side  forces  or  yawing  moments  are  induced. 

(3)  Regime  III  (20°<  <  60°);  At  high  angles  of  attack,  crossflow  effects 

start  to  dominate  and  vortices  are  now  shed  asymmetrically.  These 
vortices  induce  side  forces  (out-of-plane  forces)  and  yawing  moments. 
The  more  asymmetric  the  vortex,  the  greater  the  side  force  magnitude. 
There  are  some  flow  instabilities  toward  the  higher  end  of  this  angle  of 
attack  range. 
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(4)  Regime  IV  (60°<  «  <  90°):  At  very  high  angles  of  attack,  the  crossflow 
completely  dominates  and  flow  separation  is  unsteady.  The  Reynolds 
number,  Mach  number  and  geometry  determine  whether  the  boundary 
layer  is  shed  as  a  von  Karman  vortex  street  or  a  random  wake-like 
flow.  [Refs.  1,  2  and  3] 

Figure  1  shows  sketches  of  all  four  vortex-generation  regimes. 
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Figure  1.  Airflow  Regimes  [Ref.  2] 


2.  Asymmetric  Vortex  Theory 

The  principal  cause  of  the  formation  of  asymmetric  vortices  is  still 
not  completely  understood  and  may  be  attributed  to  many  factors.  One  idea  is 
that  boundary-layer-induced  asymmetry  in  the  location  of  flow  separation 
causes  the  vortex  flowfield  to  become  asymmetric.  These  boundary  layer 
asymmetries  may  result  from  transition  and  separation  differences  on  each 
side  of  the  missile  body.  Another  proposition  is  that  a  hydrodynamic  (inviscid) 
instability  in  the  pair  of  initially  symmetric  vortices  causes  the  asymmetry. 
[Refs.  5,  6  and  7]  These  vortices,  which  increase  in  strength  with  angle  of 
attack,  interact  with  the  surrounding  potential  flowfield  to  provide  the 
asymmetric  configuration.  A  vortex-switching  phenomenon  has  also  been 
observed  in  which  the  vortex  pattern  rapidly  switches  from  an  almost  sym¬ 
metric  to  a  highly  asymmetric  configuration,  which  may  possibly  relate  to  a 
second  inviscid  solution  in  the  leeward  flowfield.  [Refs.  8  and  9]  At  any  rate, 
even  though  their  major  cause  has  not  been  determined,  the  behavior  of 
asymmetric  vortices  has  been  well  documented  for  a  large  number  of  models 
and  shapes. 

Nose-generated  asymmetric  vortices  appear  in  the  flowfield  around  an 
ogive-nosed,  slender,  cylindrical  body.  The  vortex  formation  occurs  along  the 
entire  body  length  and  induces  a  significant  side  force  on  the  body.  Figure  2 
shows  this  vortex  flow  along  the  length  of  an  unyawed,  slender  nose  cylinder. 
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Figure  2.  Vortex  Flow  About  a  Slender  Nose  Cylinder  [Ref.  1] 


3.  Two-Dimensional  Crossflow 

Airflow  over  a  slender  body  can  be  divided  into  normal  and  axial 
components.  Axial  flow  is  along  the  missile  body  length  while  crossflow  is 
essentially  a  two-dimensional  flow  normal  to  a  cylinder.  The  crossflow 
analogy  provides  information  for  cylinder  lift  and  drag  which  act  in  the 
crossflow  direction.  Depending  on  the  type  of  flow  separation  on  either  side  of 
the  cylinder,  side  forces  (at  right  angles  to  the  crossflow)  may  exist.  [Ref.  2] 

Boundary  layer  transition  and  separation  mechanisms  may  provide 
an  explanation  for  flow  separation  and  subsequent  asymmetric  vortex  genera¬ 
tion.  The  primary  factor  which  influences  the  separation  location  of  the 
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boundary  layer  is  the  crossflow  Reynolds  number.  Other  factors  include  sur¬ 
face  roughness  and  turbulence. 

Flow  around  a  cylinder  in  incompressible  flow  can  be  classified  into 
four  regions,  represented  by  differing  flow  separation  and  drag  behavior,  as 
depicted  in  Figure  3.  [Refs.  10,  11  and  2]  In  the  subcritica!  range,  the  bound¬ 
ary  layer  is  laminar  and  flow  separation  occurs  close  to  the  lateral  meridian 


Figure  3.  Two-Dimensional  Crossflow  About  a  Cylinder  [Ref.  2] 
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where  the  angle  (0)  from  the  crossflow  direction  varies  from  80°  to  90°. 
[Ref.  12]  In  the  critical  Reynolds  number  region,  a  drag  bucket  is  produced. 
The  laminar  boundary  layer  separates  from  the  body  at  <j)  =  90°  and  reat¬ 
taches  as  a  turbulent  boundary  layer  which  is  more  energetic.  Separation  is 
delayed  to  0  «  140°,  resulting  in  a  reduction  of  the  drag.  [Ref.  3]  A  laminar 
bubble  is  formed  between  the  laminar  separation  and  the  turbulent  reat¬ 
tachment.  At  this  point,  the  flow  separation  may  easily  fluctuate  from  critical 
to  subcritical  for  small  changes  in  Reynolds  n limber.  From  Figure  4  [Ref.  10], 
when  critical  separation  exists  on  one  side  of  the  body  while  subcritical  sepa¬ 
ration  is  on  the  other,  a  large  difference  in  0  is  possible.  Therefore,  vortices 
will  be  at  maximum  asymmetry  and  the  side  force  at  the  highest  magnitude. 
Since  maximum  side  force  occurs  at  the  critical  Reynolds  number,  it  is  a 
noteworthy  parameter.  [Ref.  13] 

As  the  Reynolds  number  is  further  increased,  turbulent  separation 
moves  forward  to  0  <  140°  and  the  laminar  bubble  no  longer  exists.  At 
0  >  140°,  asymmetric  vortices  are  ineffective  at  producing  a  significant  side 
force,  thus  the  sudden  decrease  in  magnitude  (Cy/Cn)  as  shown  in  Figure  4. 
For  supercritical  and  transcritical  Reynolds  numbers,  the  laminar  transition 
point  moves  towards  0  =  0°  and  turbulent  separation  occurs  at  0  =  100°.  The 
asymmetric  transcritical  separation  point  moves  towards  the  lateral  merid- 
ian;  where  the  vortices  once  again  produce  a  significant  side  force.  [Refs.  5 
and  2]  The  Reynolds  number  provides  the  greatest  influence  on  the  normal 
force  and  drag  characteristics,  particularly  within  the  critical  range  where 
the  maximum  normalized  side  force  and  maximum  vortex  asymmetry  occur. 
[Ref.  12] 
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Figure  4.  Side  Force  to  Normal  Force  Ratio  [Ref.  5] 


Another  study  by  Lamont  [Ref.  14]  describes  a  different  effect  of 
Reynolds  number  on  the  maximum  side  force  as  illustrated  by  Figure  5, 
where  the  side  force  at  an  angle  of  attack  of  55°  is  plotted.  The  maximum  side 
force  falls  from  a  high  value  at  laminar  separation  to  almost  zero  in  the 
middle  of  the  transition  region  before  climbing  again  to  higher  levels  at  fully 
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Figure  5.  Effect  of  Reynold's  Number  of  Maximum 
Side  Force  at  =<  =  55°  [Ref.  14] 


turbulent  separation.  Thus,  there  appears  to  be  two  different  mechanisms  for 
producing  asymmetric  flow  and,  hence,  a  side  force  on  an  ogive-cylinder.  One 
mechanism  operates  in  both  the  laminar  and  the  fully  turbulent  separation 
regimes,  in  which  the  side  force  results  from  asymmetric  vortex  patterns  in 
the  wake  of  a  body.  The  other  mechanism  occurs  only  in  the  transitional 
separation  regime.  Here,  the  Reynolds  number  at  which  the  near-zero  side 
forces  were  recorded,  is  the  same  range  of  Reynolds  number  in  which  the 
minimum  drag  coefficient  on  a  2-D  cylinder  occurs  and  in  which  no  coherent 
vortex  shedding  can  be  detected. 
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4.  Three-Dimensional  Vortices 


The  missile  nose  geometry  is  an  important  factor  in  vortex  genera¬ 
tion  and  disposition  since  vortices  shed  at  the  nose  tend  to  dominate  other 
vortices  along  the  body  length,  [Ref.  15,  8  and  16] 

Nose-generated  vortices  are  sensitive  to  the  nose  roll  angle  due  to 
surface  imperfections  and  nose  geometric  deviations  [Ref.  17].  Rabang  varied 
the  ’•oil  angle  in  45°  increments  and  investigated  the  resulting  side  force  coef¬ 
ficients,  shown  in  Figure  6.  The  vortex  system  generated  by  the  nose  domi¬ 
nates  afterbody  vortices  for  body  configurations  with  and  without  wings 
regardless  of  the  turbulence  conditions  [Refs.  2,  18,  and  19], 


Figure  6.  Side  Force  Variations  With  Nose  Roll  Angle  [Ref.  2] 
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Nose  geometry  may  be  pointed  or  blunt  for  cones  and  ogives.  For 
pointed  noses,  angle  of  attack  for  the  onset  of  asymmetric  vortices  («av)  is  a 
function  of  the  semi-vertex  angle  (0a).  Asymmetric  vortices  start  at  the  nose 
and  are  rapidly  shed,  yielding  unsteady  side  forces.  [Ref.  3]  At  all  Mach 
numbers,  asymmetric  vortex  shedding  starts  when  «av  is  greater  than  0a- 
For  a  conical  nose: 

eA  =s  oc  Av/  2  (approximation)  (1) 


For  a  tangent  ogive  nose: 


0 .  =  tan 
A 


-1 


Vd 


l.:/d-0.25 

N 


(2) 


where  In  is  the  nose  length  and  d  is  the  base  diameter,  or 

0a  =  Wd  (approximation)  (3) 

Nose  fineness  ratio  also  affects  the  asymmetric  vortex  induced  side 
forces  in  that  as  this  ratio  increases,  the  side  force  also  increases.  With  an 
increasing  ratio,  both  the  nose  apex  angle  (0a)  and  the  angle  of  attack  for  the 
onset  of  asymmetric  vortices  («av)  will  decrease,  making  the  missile  more 
susceptible  to  induced  side  forces  at  lower  angles  of  attack.  [Refs.  15  and  20] 
Decreasing  nose  fineness  ratio  has  been  found  to  be  more  beneficial  in  reduc¬ 
ing  side  forces  than  blunting  the  nose  [Ref.  15].  Side  force  decreases  with  an 
increasing  Mach  number  to  nearly  zero  at  supersonic  speeds  [Refs.  18  and 
21]. 
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C.  TURBULENCE 

Turbulence  denotes  the  presence  of  random,  short  duration  variations  in 
a  flowfield  with  a  given  mean  velocity.  When  calculating  turbulence  effects  on 
a  body  in  the  flowfield,  a  comparison  between  the  scale  of  the  body  and  that  of 
the  turbulence  must  be  made.  The  energy  in  the  turbulence  flowfield  should 
also  be  considered.  [Ref.  1] 

Turbulence  intensity,  Tu,  is  the  measure  of  the  relative  magnitude  of 
velocity  fluctuations  in  the  flowfield.  For  a  horizontal  flowfield  or  cross  wind, 
it  is  the  ratio  of  the  root-mean-square  (rms)  velocity  fluctuation,  u’,  to  the 
mean  velocity  component  in  the  flowfield,  U°°. 

Tu  =  u'/TJoo  (4) 

Turbulence  length  scales  describe  the  time-averaged  measure  of  the  size 
of  the  constantly  changing  fluid  disturbance  eddies.  An  increase  in  the  spatial 
length  of  the  turbulence  corresponds  to  an  increase  in  the  time  the  body  is 
exposed  to  the  fluctuation.  Large  and  small  scale  turbulence  length  scales  are 
both  found  in  a  flowfield. 

From  a  single  source,  the  cascade  effect  produces  turbulence  eddies  of 
different  length  scales.  This  "cascade”  effect  is  caused  by  a  strain  in  one  direc¬ 
tion  (x,  y,  z  plane)  which  affects  the  orthogonal  components  due  to  the  conser¬ 
vation  of  angular  momentum.  For  example,  an  increase  in  the  x  and  y  velocity 
components  of  a  vortex  rotating  in  the  x-y  plane  will  have  an  effect  on  the  ve¬ 
locity  and  length  scales  of  the  y-z  components.  [Ref.  1]  Cascading  continues 
until  the  smaller  eddies  disappear  due  to  the  viscosity.  As  turbulence 
decreases,  the  energy  transfer  decreases  and  the  individual  intensities  of 
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each  eddy  will  decrease  at  a  faster  rate.  [Ref.  22]  Thus,  the  larger  scale  turbu¬ 
lence  predominates. 

The  length  scale  to  body  size  ratio  may  determine  the  manner  in  which 
turbulence  affects  the  VLSAM  flowfield.  It  may  be  compared  to  missile 
length,  Lu:Li,  or  missile  diameter  Lu:Ld.  [Refs.  1  and  2] 

For  length  scales  much  greater  than  the  body,  LU»L],  the  effect  is  like  a 
steady-state  flowfield,  where  deviations  in  speed  and  direction  would  be  of 
long  duration.  The  flowfield  effects  on  vortex  development  are  mainly  domi¬ 
nated  by  the  same  factors  and  conditions  as  for  a  two-dimensional  cylinder. 

In  contrast,  unwanted  rolling,  pitching  and  yawing  motion  of  the  body  is 
primarily  caused  when  the  turbulence  length  scale  is  comparable  to  the  body 
length,  Lu=Li.  [Ref.  23]  The  flowfield  is  distinctly  non-steady  for  this  case. 

When  the  length  scale  is  of  a  dimension  much  smaller  than  the  body, 
most  significantly,  when  it  is  smaller  than  the  missile  diameter,  Lu«Ld,  it 
has  a  magnitude  comparable  to  the  boundary  layer  thickness  on  the  missile 
surface.  Thus,  boundary  layer  development  and  flow  separation  over  the  body 
may  be  affected  by  the  presence  of  small  scale  turbulence.  An  increase  in  tur¬ 
bulence  intensity  with  a  length  scale  on  the  order  of  the  boundary-layer  scale 
tends  to  reduce  the  magnitude  of  induced  side  forces.  [Refs.  2,  24,  and  25]  A 
goal  of  the  current  investigation  is  to  determine  the  effect  of  vortex-scale 
turbulence  on  the  asymmetric  vortices  and  resulting  induced  side  forces. 

D.  LIFTING  SURFACE  EFFECTS 

The  complete  vortex  structure  of  the  missile  is  a  net  result  of  the  indi¬ 
vidual  contributions  from  the  body,  wings,  strakes  and  tails.  In  general,  mis¬ 
siles  use  low  aspect  ratio  wings  (when  compared  with  aircraft).  Since  some 
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missiles  have  wing  spans  that  approach  body  diameter,  it  is  important  to 
consider  the  joint  effects  from  a  wing-body  combination.  Nose  vortices  dictate 
flow  behavior  over  a  missile  body  at  high  angles  of  attack  and,  consequently, 
these  vortices  may  also  be  felt  by  the  wings.  Nose  and  body  vortices  move 
away  from  a  missile  body  without  wings  but,  when  wings  are  added,  they 
move  closer  to  the  body  This  result  is  comparable  to  increasing  the  effective 
angle  of  attack  causing  unsteady  asymmetric  vortices.  For  wings  with  low 
aspect  ratio,  a  major  portion  of  the  lift  produced  by  the  wing  will  be  a  result  of 
vortex  lift.  The  net  effect  of  the  wing-body  combination  appears  to  be  a  reduc¬ 
tion  in  the  effective  angle  of  attack  for  the  onset  of  asymmetric  vortices  and 
side  forces.  [Refs.  2,  3,  and  26] 

Vortex  lift  effects  are  improved  by  incorporating  strakes  with  low  aspect 
ratio  wings.  The  strakes  produce  additional  strong  vortices.  Some  researchers 
have  found  that  placing  long  strakes  on  a  missile  would  induce  interference 
with  the  crossflow  component  around  the  body,  thus  decreasing  the  effect  of 
the  forces  and  moments  generated  by  asymmetric  body  vortices.  [Ref.  6] 
Rabang  has  shown  that  the  addition  of  typical  VLSAM  wings  and  strakes 
tend  to  preserve  the  induced  side  force  for  all  levels  of  turbulence  intensities 
and  length  scales.  [Ref.  2] 

The  addition  of  tails  has  a  minor  influence  on  forebody  flowfields  and 
maximum  side  forces,  in  particular,  at  low  angles  of  attack.  At  higher  angles 
of  attack,  nose  and  wing  vortices  may  have  a  slight  effect  on  the  tailflow, 
depending  on  the  wing  placement  and  afterbody  length  between  the  wings 
and  tails. 


15 


E.  VLSAM  LAUNCH  ENVIRONMENT 
1.  Marine  Environment 

Turbulence  conditions  which  exist  within  the  atmospheric  boundary 
layer  (ABL)  may  significantly  impact  the  VLSAM.  The  ABL  is  the  lowest  por¬ 
tion  of  the  atmosphere  and  is  formed  by  its  interaction  with  the  surface  over 
which  it  flows.  Turbulence  in  this  layer  is  the  result  of  the  transfer  of  heat, 
momentum  and  mass. 

The  surface  layer,  the  lowest  segment  of  the  ABL,  can  vary  in 
height  from  5  to  200  meters  but  is  typically  on  the  order  of  50  meters.  It  is 
also  characterized  by  mechanically  produced,  small-scale  turbulence  resulting 
from  surface  roughness  or  friction  from  waves  on  the  ocean  surface.  This 
small-scale  turbulence  is  larger  than  the  missile  length.  This  region  is 
described  by  variations  in  wind  speed,  nearly  vertical  heat  and  mass  fluxes, 
and  other  meteorological  fluctuations  with  height.  [Refs.  2?  and  23]  Further¬ 
more,  the  majority  of  the  flow  in  the  surface  layer  itself  can  be  considered 
horizontally  homogeneous  and  two-dimensional  [Ref.  23]. 

A  measure  of  the  roughness  of  the  surface  is  called  the  roughness 
length,  Z0,  which  is  a  function  of  the  mean  wind  velocity  at  various  heights 
above  the  surface.  By  combining  the  roughness  length  with  the  elevation  and 
wind  speed,  both  the  turbulence  intensity  and  length  scale  can  be  empirically 
determined.  [Refs.  2,  3  and  28]  Typical  open  ocean  surface  roughness  lengths 
are  in  the  range  of  0.001<Zo<0.01,  with  Zq  in  meters.  For  a  10  meter  elevation 
and  a  mean  wind  speed  of  25  m/sec  in  a  neutral  atmosphere,  turbulence 
intensities  may  range  from  13  to  17  percent.  [Ref.  28]  The  longitudinal 
turbulence  length  scale  would  then  range  from  262<LU<295  feet  (80<Lu<  90 
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meters).  Therefore,  for  a  typical  missile  with  a  1.1  foot  diameter,  the 
turbulence  length  scale  to  missile  diameter  length  scale  ratio  is  Lu:Ld«280:l. 
[Ref.  28]  This  represents  a  length  scale  very  much  larger  than  a  conventional 
missile  length  and,  therefore,  would  have  little  effect  on  its  boundary  layer 
development.  However,  the  cascade  effect  from  large  scale  turbulence  and 
from  crosswind  interaction  with  a  ship's  superstructure  allows  length  scales, 
initially  much  larger  (85  meters)  than  the  dimension  of  a  missile,  to  decrease 
(cascade)  to  scales  where  they  could  affect  the  missile  boundary  layer 
development  and  the  development  of  vortices  from  the  missile  nose.  The 
actual  amount  of  such  small-scale  turbulence  present  in  the  marine 
atmosphere  is  largely  unknown,  however.  [Ref.  23] 

2.  Launch  and  Crosswind  Velocities 

A  typical  VLSAM  at  launch  (vertical  velocity  =  164  ft/sec)  is  still 
well  within  the  surface  layer  environment  and  is  subject  to  both  crosswind 
and  turbulence  effects.  [Ref.  29]  Crosswind  velocity  depends  on  both  the 
ambient  wind  speed  and  the  speed  of  the  launch  platform.  A  ship  speed  of  20 
knots  with  a  head  wind  of  20m/sec,  combined  with  the  VLSAM  launch 
velocity,  results  in  an  effective  angle  of  attack  of  31°  at  191  ft/sec,  which 
places  the  missile  in  the  asymmetric  vortex  region  (Regime  III)  almost  imme¬ 
diately  after  clearing  the  exit  plane  of  the  launcher,  0.2  seconds  after  launch. 
[Refs.  1  and  2] 

When  the  missile  first  leaves  its  launcher  it  will  experience  an  even 
greater  effective  angle  of  attack  due  to  its  slower  velocity.  A  ship's  hull  and 
superstructure  can  dictate  changing  flow  fields  and  turbulence  at  this  initial 
launch  altitude.  The  ship  airwake  may  increase  crosswind  velocities  and 
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cause  significant  crosswind  gradients  in  the  flowfield,  thus  increasing  turbu¬ 
lence  intensities  while  decreasing  turbulence  length  scales.  [Ref.  2] 

Later  in  its  launch  profile,  when  the  VLSAM  pitches  over  towards  a 
target,  it  may  reach  effective  angles  of  attack  of  up  to  50°.  [Refs.  2  and  30] 
Thus,  during  the  launch  phase  there  exists  a  definite  possibility  of  asymmet¬ 
ric  vortex  induced  side  forces  on  a  VLSAM. 

3.  Additional  Launch  Considerations 

During  launch,  a  missile  is  influenced  by  many  factors  which  are 
dependent  on  missile  design,  ship's  orientation,  launcher  mechanics  and 
missile  flight  control.  Shipboard  roll,  pitch  and  yaw  are  directly  transmitted 
to  the  launch  platform  and  must  be  taken  into  account.  Inherent  factors,  such 
as  the  plume  (jet)  effect  of  the  missile's  engine  and  blast  effects  of  the  vented 
exhaust  gases,  can  also  affect  VLSAM  aerodynamics.  Exhaust  gases  can 
impinge  directly  on  the  missile  surfaces  or  they  can  impact  the  flowfield  in 
which  the  missile  is  launched,  especially  if  the  gases  are  vented  upward  into 
the  vicinity  of  the  accelerated  missile.  [Ref.  1]  The  manner  in  which  control 
systems  respond  to  missile  orientation  changes  is  another  factor.  Should  the 
missile  change  its  flight  attitude,  the  flowfield  around  it  will  also  be  altered. 
Obviously,  the  considerations  discussed  in  this  section  are  not  all-encompass¬ 
ing.  There  are  many  other  factors  that  affect  missile  flight  behavior  during 
launch,  but  they  are  beyond  the  scope  of  this  thesis  and  will  not  be  included. 
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II.  EXPERIMENT  AND  PROCEDURES 


A.  PURPOSE 

In  this  study,  the  location  and  intensity  of  asymmetric  vortices  in  the 
wake  of  the  VLSAM  model  were  determined  for  varying  levels  of  turbulence. 
This  turbulence  was  generated  by  the  placement  of  a  series  of  grids  in  the 
wind  tunnel.  The  vortices  were  displayed  by  velocity  mapping  and  pressure 
contours.  To  accomplish  this,  wind  tunnel  flowfield  pressure  measurements 
were  taken  for  a  specified  survey  grid  using  a  scamvalve/probe  data  acquisi¬ 
tion  system. 

FiBure  7  [Ref.  3]  shows  the  planar  survey  grid,  the  x-y  plane,  which  was 
perpendicular  to  the  freestream  velocity  and  located  10.5  inches  downstream 
from  the  missile  model's  nose.  The  model  body  was  centered  on  the  y  dimen¬ 
sion.  For  the  actual  data  acquisition  runs,  23  points  were  measured  in  the  y 
direction.  There  were  11  points  above  and  below  the  model  centerline,  with 
point  12  directly  at  the  centerline.  Along  the  x  axis,  13  points  were  measured, 
the  x-y  dimension  for  the  experiment  was  3x5.5  inches,  with  a  step  distance 
(increment)  of  0.25  inch  (22x0.25=5.5  and  12x0.25=3  inches).  This  dimension 
covered  the  main  portion  of  asymmetric  vortices. 

Pressure  measurements  were  obtained  by  the  5-hole  probe  throughout 
the  survey  plane.  The  data  from  the  pressure  probe  was  reduced  through  the 
use  of  computer  programs  to  obtain  isobars  of  total  pressure  coefficient  and 
static  pressure  coefficient,  and  to  map  the  crossflow  velocity  vectors.  These 


19 


Figure  7.  The  Planar  Survey  Grid  [Ref.  3] 
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results  were  correlated  with  the  force  measurements  of  Rabang  and  with  the 
previous  experiments  by  Lung  to  provide  a  greater  understanding  of  the 
vortex  flowfield.  The  following  sections  further  discuss  the  equipment  and 
software  used,  and  the  experimental  procedures  followed. 

B.  APPARATUS 

Information  about  the  construction,  specifications  and  configurations  of 
the  major  pieces  of  equipment  used  in  this  study  is  described  in  this  section. 

1.  Wind  Tunnel 

The  low-speed,  single  return,  horizontal-flow  tunnel  located  in 
Halligan  Hall  at  NPS  was  utilized.  (Figure  8,  [Ref.  31])  .  It  is  powered  by  a 
100  horsepower  electric  motor  coupled  to  a  three-blade  variable-pitch  fan  via 
a  4-speed  Dodge  truck  transmission.  Aft  of  the  fan  blades  are  a  set  of  stator 
blades  which  help  straighten  flow.  Two  fine  wire  mesh  screens  located  up¬ 
stream  of  the  settling  chamber  plus  turning  vanes  at  all  four  comers  reduce 
turbulence.  A  heavy  wire  screen  behind  the  test  section  prevents  foreign 
object  damage  to  the  fan  blades  [Refs.  1,  2  and  31]  The  tunnel  is  64  feet  long 
and  ranges  from  21.5  to  25.5  feet  wide. 

The  wind  tunnel  test  section  measures  45  inches  by  32  inches.  The 
walls  diverge  slightly  to  prevent  reduction  in  freestream  pressure  due  to 
boundary  layer  growth.  The  settling  chamber  area  to  test  section  area 
contraction  ratio  is  approximately  10:1.  Comer  fillets,  which  house  the  light¬ 
ing,  reduce  the  section  area  from  10  ft2  to  9.88  ft2.  Fillets  are  found  at  wall 
intersections  throughout  the  tunnel  to  help  reduce  boundary  layer  effects. 
[Ref.  31] 
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Figure  8.  Naval  Postgraduate  School  Wind  Tunnel  [Ref.  31] 


A  reflection  plaue  installed  in  the  test  section  reduces  the  available 
height  to  28  inches.  A  flush-mounted  turntable  allows  for  changes  in  model 
pitch  angle  or  angle  of  attack  via  a  remotely  controlled  electric  motor  beneath 
the  tunnel.  Since  the  test  section  operates  at  atmospheric  pressure,  breather 
slots  are  installed  around  the  tunnel  perimeter  to  replenish  air  lost  through 
leaks  and  to  ensure  a  uniform  test  section  pressure.  The  tunnel  was  designed 
to  provide  test  section  velocities  of  up  to  290  fl/sec.  [Ref  31] 
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Wind  tunnel  temperature  is  measured  by  a  dial  thermometer 
extending  into  the  settling  chamber.  Dynamic  pressure  is  measured  by  the 
static  pressure  difference  between  the  test  section  and  the  settling  chamber 
using  a  water  filled  manometer.  The  static  pressure  is  measured  by  four  pres¬ 
sure  taps  located  upstream  from  the  test  section  to  avoid  interference  from 
the  model.  These  taps  are  connected  via  a  common  manifold  prior  to  feeding 
into  the  manometer.  Pressure  differences  are  measured  in  centimeters  of 
water.  Equation  (5)  is  used  to  convert  to  the  actual  wind  tunnel  velocity. 
[Ref.  3] 

‘  (2X2.0475)  (P  H_0)  1^ 

U  =  - —  (5) 

m  (K)(p) 


where: 

Um  =  measured  velocity  (ft/sec) 

2.0475  =  conversion  factor 
P cm  H2O  =  manometer  reading 
K  =  calibration  factor  (for  specific  grid) 
p  =  air  density  (lb/ft2) 

2.  Turbulence-Generating  Grids 

Four  grids  are  used  to  create  turbulence  of  varying  intensities  and 
length  scales.  Each  is  mounted  in  a  wooden  frame  and  placed  73  inches  for¬ 
ward  of  the  pivot  axis  of  the  model  support  system  (see  Figure  9).  Three  of  the 
grids  are  constructed  from  wood  and  the  fourth  is  made  of  wire.  Grid  specifi¬ 
cations  are  listed  in  Table  1  and  are  shown  in  Figure  10.  They  are 
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Figure  9.  Planview  of  VLSAM  Model  With  Pressure  Probe  and  Grid  In  the 
Test  Section  of  the  Wind  Tunnel  (not  drawn  to  scale)  [Ref.  2] 
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TABLE  1.  GRID  SPECIFICATIONS  [REF.  2] 


Grid 

Mesh  Width  (in.) 

Bar  Diameter  (in.) 

Mesh/Diameter 

Material 

One 

5.00 

1.00 

5  j 

Wood 

Two 

3.75 

0.75 

5 

Wood 

Three 

2.50 

5 

Wood 

Four 

1.00 

0.0625 

16 

Wire 

Figure  10.  Square-Mesh  Turbulence-Generating  Grid  [Ref.  3] 
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square-mesh  square-bar  biplanar  grids  which  generate  nearly  isotropic  homo¬ 
geneous  turbulence.  [Ref.  29J  Roane  measured  turbulence  intensities  and 
estimated  length  scales,  shown  in  Figures  11  and  12.  [Ref.  1]  The  grid  turbu¬ 
lence  parameters  taken  by  Roane  are  summarized  in  Table  2.  Grid  turbulence 
effects,  with  respect  to  changing  length  scales  at  constant  intensity  or  con¬ 
stant  length  scales  with  changing  intensities,  can  not  be  investigated  with  the 
present  grid  geometries.  [Ref.  3]  Figure  13  shows  photographs  of  the  grids. 


Figure  11.  Grid  Turbulence  Intensities  (dashed  line  indicates 
model  pivot  axis)  [Ref.  2] 
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Figure  12.  Grid  Turbulence  Length  Scales  [Ref.  2] 


TABLE  2.  GRID  TURBULENCE  PARAMETERS 
(AT  MODEL  PIVOT  AXIS)  [REF.  1] 


Grid 

Intensity 

(percent) 

Length  Scale 
(in.) 

Turbulence/ 
Model  Dia. 

Dynamic  Pressure 
(Ib/ft1) 

One 

3.31 

1.84 

1.05 

15.35 

Two 

2.78 

1.56 

0.89 

14.88 

Three 

1.88 

1.08 

0.62 

16.38 

Four 

0.47 

0.27 

0.15 

15.61 

None 

0.23 

- 

'  - 

15.85 

Figure  13.  Turbulence-Generating  Grids 
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3.  VLSAM  Model  and  Support  Equipment 

The  model  was  designed  to  represent  a  current  cruciform  tail-con¬ 
trol  missile  with  very  low  aspect  ratio  wings  (long  dorsal  fins).  It  was  con¬ 
structed  from  6061  and  2024  aluminum  alloy  by  NPS  personnel.  [Ref.  1] 

The  hollow  cylinder  body  section  contains  locating  pin  attachment 
points  for  the  balance,  sleeve,  wings  and  tails.  The  machined  sleeve  provides 
a  close  tolerance  fit  between  the  balance  gage  and  the  interior  of  the  model. 
[Ref.  3]  Both  body  roll  angle  and  nose  roll  angle  may  be  varied  in  45°  incre¬ 
ments.  The  wings  with  strakes  and  the  tail  control  fins  are  rigidly  connected 
to  the  model  body  by  countersunk  screws.  Figure  14  depicts  the  dimensions 
and  specifications  of  the  VLSAM  model.  [Ref.  2]  The  model's  surface  is 
polished  and  free  of  protruberances. 

The  model  support  is  rigidly  fixed  in  the  test  section  by  the  reflec¬ 
tion  plane  turntable  at  the  base  and  an  aluminum  reinforced  clear  plexiglass 
section  at  the  top.  The  pivot  point  of  this  rotating  support  coincides  with  the 
anp-oximate  center  of  the  VLSAM  model.  The  plexiglass  has  three  slots  (7-,  8- 
and  10-inches  long)  cut  in  it,  each  5/4  inches  wide.  These  slots  correspond  to 
the  positions  of  model  length  to  diameter  ratios  of  3,  6  and  9;  i.e.,  5.25,  10.5 
and  15.75  inches  from  the  nose.  [Refs.  2  and  3] 

4.  Veimex  8300  3-D  Traverser 

The  Veimex  8300  is  composed  of  a  motor  controller  assembly  and  a 
traversing  assembly,  and  uses  three  microcomputer-controlled  stepping 
motors  (one  for  each  axis  of  movement).  The  motor  controller  assembly  is 
capable  of  interpreting  motor  movement  commands  from  a  host  computer, 
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Total  length 
Base  diameter 
Length/diameter  ratio 
Ogive  nose  length 
Ogive/diameter  ratio 
Wing  span, 'root  chord 
Tail  span.'root  chord 
Center  of  pressure 


22,85  in. 

1.75  in. 

13.06 
4.0  in. 

2.29 

3.13  in./ 13.55  in. 

5.50  in. /1. 70  in. 

13.5  inches  aft  of  nose  tip  (approx.) 


NAVAL  POSTGRADUATE  SCHOOL 
SURFACE-TO-AIR 
MISSILE  MODEL 


Figure  14.  Specifications  of  VLSAM  Model  [Ref.  2] 


programmable  control  or  terminal.  Software  commands  allow  the  operator  to 
select  motor  variables  such  as  velocity,  acceleration,  increment  distance  and 
units  (motor  steps  or  inches).  [Ref.  32] 
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The  stainless  steel  and  aluminum  traverser  assembly  (Figure  15) 
consists  of  three  separate  motor/jackscrew  assemblies. 


Figure  15.  Velmex  8300  Traversing  Assembly 

The  traverser  was  mounted  to  existing  hardware  on  top  of  the  tun¬ 
nel  so  as  to  minimize  tunnel-induced  vibrations.  A  5-hole  pressure  probe, 
attached  to  the  8300  control  drive,  can  be  accurately  and  effectively  moved 
through  the  test  section. 
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5.  5-Hole  Pressure  Probe 

The  three-dimensional  5-hole  probe,  Figure  16  [Ref.  33],  is  made  of 
corrosion  resistant  non-magnetic  stainless  steel.  It  is  0.125  inch  in  diameter 
and  24  inches  in  total  length  with  22  inches  of  reinforcement  tubing.  The 
probe  has  five  measuring  holes  located  on  its  prism-shaped  tip.  A  centrally 
located  hole  (Pi)  measures  total  pressure,  while  two  lateral  pressure  holes 
(P2,  P3)  are  used  to  determine  yaw  angle  of  flow.  Pitch  angle  is  determined  by 
pressure  holes  (P4,  P5)  located  above  and  below  the  total  pressure  hole.  The 
probe  is  usable  for  speeds  up  to  Mach  0.7. 

The  speed  of  reading  depends  on  the  length  and  diameter  of  the 
pressure  passage  inside  the  probe,  the  size  of  the  pressure  tubes  to  the 
manometer,  and  the  displacement  volume  of  the  manometer.  [Ref.  33]  For 
smaller  diameter  tubes,  the  time  constant  increases  rapidly.  For  this  experi¬ 
ment,  the  tube  diameter  was  1/4-inch  O.D.  and  the  tube  lengths  were  three 
feet,  so  the  time  delay  was  about  0.15-0.26  second.  [Ref.  3] 

6.  Scanivalve  and  HP  Data  Acquisition  System 

One  48-port  scanivalve  was  used  to  measure  each  of  the  5-hole 
probe  pressures.  The  Hewlett-Packard  (HP)  data  acquisition  system  consists 
of  a  combination  of  hardware  and  software  that  enables  the  IBM  PC-AT  com¬ 
puter  to  act  as  a  fully  automated  instrumentation  system.  [Ref.  34]  Individual 
instruments  include  the  Relay  Multiplexer,  Digital  Multimeter  and  Relay 
Actuator. 
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Figure  16.  The  5-Hole  Pressure  Probe  and  Measuring  Tip  [Ref.  33] 


The  scanivalve  mechanism  puts  out  a  7-bit  binary  coded  decimal 
(BCD)  signal  that  corresponds  to  the  port  (1-48)  currently  connected  to  the 
scanivalve  transducer.  This  allows  remote  electronic  monitoring  of  the  port 
assembly  configuration.  [Ref  3]  The  scanivalve  consists  of  a  transducer, 
motor  drive,  port  assembly  and  solenoid  controller,  which  actually  regulates 
the  scanivalve.  Two  commands  allowed  by  this  solenoid  are  STEP,  which 
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moves  the  scanivalve  to  one  port  location,  and  HOME,  which  sends  the  scani- 
valve  to  port  number  48.  The  Relay  Actuator  is  used  solely  in  controlling  the 
scanivalve  to  STEP  or  HOME.  [Ref.  3] 

The  HP  Data  Acquisition  System  is  shown  in  Figure  17.  The  scanivalve 
signal,  containing  probe  port  pressure  information,  passes  to  the  Relay  Multi¬ 
plexer  which  provides  one  common  output  channel  for  the  Digital  Multimeter 
(DMM).  The  signal  is  conditioned  by  a  low  pass  filter  prior  to  being  measured 
by  the  multimeter.  The  DMM  automatically  converts  input  analog  voltage 
signals  into  a  digital  (binary)  form  which  can  be  read  by  the  computer. 

C.  EXPERIMENTAL  CONDITIONS 

To  facilitate  data  correlation,  the  conditions  for  this  experiment  were 
similar  to  the  conditions  of  the  previous  studies  by  Lung  and  Rabang. 

(1)  Test  section  reference  dynamic  pressures  were  set  at  7.2  cm  H2O  for  the 
no  grid  run  and  at  10.0  cm  H2O  for  grid  #3,  which  yielded  a  subcritical 
Reynolds  number  of  Re=l.lxl02 3 * 5.  These  reference  pressures  are  the 
same  as  those  used  in  the  turbulence  mapping  by  Roane  [Ref.  1]  and 
were  duplicated  in  this  study  in  order  to  ensure  comparable  test 
section  velocities  and  turbulence  grid  length  scales  and  intensities. 

(2)  The  VLSAM  model  nose  geometry  was  held  fixed  at  nose  position  eight, 
which  Rabang  showed  gave  the  maximum  side  force  magnitude. 
[Ref.  2] 

(3)  Afterbody  roll  angles  were  as  follows: 

•  Body  A:  wings  and  tails  at  roll  angle  4>r=0°  in  a  "+"  configuration 

•  Body  B:  no  wings  or  tails  at  roll  angle  4>r=45° 

•  Body  C:  wings  and  tails  at  roll  angle  4>r=45°  in  a  "x"  configuration. 

Figure  18  shows  these  three  configurations.  Only  Bodies  A  and  C  were 
tested  in  this  study. 
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Figure  18.  VLSAM  Model  Body  Configurations  [Ref.  2] 


(4)  Model  blockage  factor  corrections  were  calculated  by  Rabang  for  each 
body  configuration.  These  factors  are  a  function  of  the  model  angle  of 
attack.  For  this  experiment,  the  angle  of  attack  was  fixed  at  50°,  where 
the  total  blockage  correction  £  equalled  0.0123.  This  factor  was  imple¬ 
mented  in  data  conversion  programs. 

(5)  The  longitudinal  position  for  data  acquisition  was  at  a  length/body 
diameter  ratio  of  6,  which  was  10.5  inches  from  the  nose  of  the  missile 
model. 

(6)  Wind  tunnel  temperatures  were  not  allowed  to  vary  by  more  than  20°F 
from  the  beginning  to  the  end  of  a  run.  Wind  tunnel  settling  chamber 
temperatures  tended  to  rise  quickly  due  to  air  friction,  particularly 
when  the  grid  was  added.  When  temperatures  were  excessive,  the  tests 
were  stopped  and  the  air  in  the  tunnel  was  circulated  until  it  cooled 
down  sufficiently  before  tests  were  continued. 


D.  SOFTWARE  AND  PROCEDURES 


In  order  to  correlate  results  with  previous  data,  the  computer  programs 
used  (or  developed)  by  Lung  to  acquire  and  reduce  data  were  also  used  for 
this  study.  Figure  19  [Ref.  3]  is  a  schematic  flowchart  of  the  various  programs 
and  their  resulting  data  files.  The  following  sections  provide  further  elabora¬ 
tion  on  these  programs. 


Figure  19.  Program/Data  File  Flowchart  [Ref.  3] 
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1.  PPROBE  Program  (Data  Acquisition) 

The  BASICA  application  program  which  runs  the  VLSAM  experi¬ 
ment  is  comprised  of  STATEFILE,  PGMSHEL  and  traverser  programs. 
STATEFILE  gives  the  computer  configurations  of  the  data  acquisition 
instruments,  while  PGMSHEL  informs  the  computer  of  all  the  functions 
available  at  e»cb  of  these  instruments.  STATEFILE  and  PGMSHEL  perform 
initialization  chores  and  allow  communication  between  the  HP  instruments 
and  the  IBM  computer.  [Ref.  34] 

The  traverser  program  is  the  actual  application  code  which  allows 
the  operator  to  precisely  control  pressure  probe  movement  by  either  manual 
or  computer-controlled  input.  The  traverser  program  was  written  by 
Kindelspire  [Ref.  35]  in  the  Advanced  Basic  Language.  The  PPROBE  program 
is  shown  in  Appendix  A. 

Manual  control  was  used  to  initialize  the  pressure  probe  position 
prior  to  collecting  data  for  the  actual  run,  which  utilized  computer-controlled 
movements.  While  in  the  manual  mode,  the  program  asked  a  series  of  ques¬ 
tions  which  enabled  the  operator  to  test  motor  movements  and  set  traverser 
motor  velocity  and  acceleration  default  values.  Through  manual  inputs,  the 
probe  was  positioned  such  that  the  Pi  (total  pressure)  hole  was  centered  on 
the  lengthwise  axis  of  the  VLSAM  model  body  and  placed  as  close  as  possible 
to  it.  From  this  point,  the  probe  was  moved  vertically  downward  to  the  posi¬ 
tion  of  the  first  point  in  the  data  collection  plane.  For  this  experiment,  the 
origin  was  located  2.75  inches  below  the  model  axis,  and  the  data  field  dimen¬ 
sions  were  2.75  inches  above  and  below  and  3  inches  outward  from  the  model. 
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After  this  initial  probe  position  was  set,  the  computer-controlled 
motor  movement  option  was  selected.  The  field  dimensions  (x,  y  coordinates), 
the  traverser  motor  step  distance,  and  the  input  file  name  were  entered  into 
the  program.  The  data  plane  was  3.0"  by  5.5"  with  a  0.25-inch  step  distance 
for  this  study.  From  this  input,  PPROBE  then  reiterated  the  total  number  of 
points  to  be  measured  (299  this  case)  and  assigned  filenames  for  each  column 
of  data.  An  example  of  how  data  files  were  named  is  as  follows: 
example:  R0A1A3 

where: 

R  =  run 

0  =  grid  number  (type) 

A  =  VLSAM  model  configuration 
1  =  field  dimension  (3,  5.5) 

A  =  step  distance  (0.25) 

3  =  test  number 


Thus,  from  the  example  above,  the  program  assigned  filenames 
R0A1A301.DAT  through  R0A1A313.DAT,  which  represented  the  13  columns 
of  data  (23  points  per  column). 

The  5-hole  pressure  probe  scale  wheel  was  then  adjusted  until  the 
P2  and  P3  (lateral)  pressures  were  nearly  equal  (nulling),  as  measured  by  a 
portable  digital  manometer/calibrator.  The  measured  yaw  angle  was  read  off 
the  wheel  and  typed  into  the  computer.  Once  this  was  accomplished, 
PPROBE  moved  the  scanivalve  from  port  1  to  port  4  via  the  Relay  Actuator. 
There  was  a  one-second  delay  to  allow  pressure  equalization  before  the 
Digital  Multimeter  sampled  the  output  voltage  from  the  scanivalve  trans¬ 
ducer  via  the  Relay  Multiplexer.  After  ten  samples  were  taken  at  port  4,  the 
Relay  Actuator  stepped  the  scanivalve  to  the  next  port  (5),  where  another  ten 
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samples  were  taken.  This  process  was  repeated  until  all  five  channel  pres¬ 
sures  (ports  4  through  8)  were  measured.  [Ref.  3]  Note  that  scanivalve  port  4 
represents  probe  pressure  Pi,  port  5  is  P2,  port  6  is  P3,  port  7  is  P4,  and  port  8 
is  P5. 

The  Relay  Actuator  homed  the  scanivalve  to  port  48  after  all  the 
pressures  were  measured,  and  then  PPROBE  displayed  the  measurements 
and  average  values  for  each  channel  on  the  computer  screen.  The  program 
either  moved  the  probe  upward  one  step  (0.25")  or  remeasured  the  same 
point,  depending  on  whether  the  data  was  within  tolerance  or  not,  as  deter¬ 
mined  by  the  operator.  For  this  study,  the  following  tolerances  were  used:  P2 
and  P3  differed  by  0.09  psf  or  less,  and  Pi  was  a  positive  number  (or  a  very 
small  negative  number  on  the  order  of  -0  5  psf).  Once  a  column  of  data  was 
measured,  PPROBE  would  store  the  average  values  for  each  point  in  a  file  (23 
pts)  and  move  the  traverser  to  the  next  column  position.  The  data  acquisition 
process  was  continued  until  all  299  points  were  completed. 

2.  CALP  Program  (Scanivalve  Transducer  Calibration) 

The  CALP  program  (Appendix  B)  was  the  other  data  acquisition 
program  utilized  in  this  study.  It  was  run  both  before  and  after  the  actual 
test  (PPROBE)  to  account  for  any  change  in  experimental  conditions  which 
might  have  occurred  over  the  8-10  hour  period  it  took  to  run  PPROBE. 

The  transducer  voltage  was  first  adjusted  to  approximately  zero 
millivolts.  Calibration  manometer  (Figure  20)  readings  were  then  entered 
into  the  computer.  The  manometer  provided  a  known  pressure  source  for 
scanivalve  calibration.  From  the  transducer  output  voltage  and  pressure  data 
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Figure  20.  Calibration  Manometer  [Ref.  3] 


provided  by  CALP,  two  calibration  curve  equations  were  calculated  (before/ 
after  PPROBE  runs).  An  averaged  slope  equation  was  then  used  in  data 
reduction  computations  to  minimize  errors.  This  slope  equation  was  used  by 
the  CONVERT  program  to  change  voltage  data  into  dynamic  pressure  data. 
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3.  CONVERT  Program 

This  program  (Appendix  C)  was  used  to  reduce  data.  It  read  the 
PPROBE  data  files  (ex:  R0A1A303.DAT)  containing  average  pressures  (PI¬ 
PS)  and  converted  them  into  x-y  coordinates,  velocity,  yaw  angle,  pitch  angle 
(alpha),  total  local  pressure,  total  pressure  coefficient,  local  static  pressure, 
and  static  pressure  coefficient.  The  x,  y  coordinates  and  yaw  angle  were  sim¬ 
ply  read  from  the  PPROBE  data  files  and  input  into  CONVERT.  The  pitch 
angle  coefficient  was  determined  by  a  ratio  of  the  pressures  measured  by  the 
5-hole  probe.  This  coefficient  was  then  used  in  an  equation  developed  by  Lung 
[Ref.  3]  using  commercial  curve-fitting  software  in  order  to  find  the  corre¬ 
sponding  pitch  angle.  Calibration  curves  provided  by  the  probe  manufacturer 
were  used  for  the  curve-fit.  The  program  also  obtained  velocity  pressure  coef¬ 
ficients  for  particular  pitch  angles,  from  which  the  local  velocity  was  calcu¬ 
lated  after  addition  of  a  wind  tunnel  calibration  factor,  K,  and  other  unit  con¬ 
version  factors.  Additionally,  the  total  and  static  pressure  coefficients  were 
obtained  for  different  pitch  angle  regions. 

Room  ambient  pressure  was  used  as  the  reference  pressure  for  these 
coefficients,  which  are  functions  of  total  and  static  pressures  and  are  non- 
dimensionalized  by  the  tunnel  dynamic  pressure: 

Cps  =  (PsL-PS)/Q  (6) 

Cpr  =  (PtL  — Pt)/Q  (7) 

where 

Cps  =  Static  pressure  coefficient 
Cpr  =  Total  pressure  coefficient 
Q  =  Freestream  dynamic  pressure 

Ps  =  Freestream  static  pressure 

Pt  =  Freestream  total  pressure 

PSL  =  Local  static  pressure 

PtL  =  Local  total  pressure 

The  actual  dynamic  pressure  is  (nearly)  the  same  for  the  two  cases. 
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For  this  experiment,  the  reference  dynamic  pressure  values  were  7.2  cm 
H2O  for  grid  0  and  10.0  cm  H2O  for  grid  3.  Temperature  input  for  the 
CONVERT  program  was  an  average  of  the  initial  and  final  wind  tunnel 
temperatures  for  the  entire  run  time.  Similarly,  barometric  pressure  values 
were  recorded  before  and  after  each  run.  CONVERT  also  added  yaw  and 
pitch  angle  (alpha)  corrections  to  the  data.  These  factors  are  +5.0°  for  yaw 
and  -17.942°  for  pitch.  They  were  determined  from  a  preliminary  run 
conducted  with  no  grid  and  no  missile  model  in  the  wind  tunnel.  Further 
explanation  of  this  run  is  discussed  in  the  Preliminary  Tests  section.  The 
output  of  the  CONVERT  program  was  stored  in  a  file  named  RESULT.DAT, 
which,  in  turn,  was  used  as  input  to  the  TECPLOT  system. 

4.  TECPLOT 

The  commercial  TECPLOT  software  system  was  used  to  generate 
crossplane  velocity  vector  plots  and  pressure  contour  plots.  These  plots  could 
be  tailored  in  many  different  ways  by  choosing  scale  factors,  arrowhead 
wedge  angles,  contour  levels  and  spacing,  and  many  other  parameters.  A 
Hewlett-Packard  7470A  x-y  pen  plotter  was  utilized  in  conjunction  with 
TECPLOT  to  provide  both  the  vector  and  contour  plots. 

E.  PRELIMINARY  TESTS 

1.  Dynamic  Pressure  Calibration 

All  of  the  turbulence  grids  were  previously  calibrated  in  the  wind 
tunnel.  Readings  from  the  tunnel  calibration  manometer  and  from  a  pitot- 
static  tube  inserted  in  the  center  of  the  test  section  were  recorded  over  a 
speed  range  and  wind  tunnel  calibration  factors  were  obtained.  These  factors 
were  used  to  adjust  the  tunnel  flow  velocity  to  the  expected  experimental 
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condition  for  the  different  grids.  [Ref.  3]  The  calibration  factors  K  are  0.8391, 
1.5084,  1.6487,  1.6545  and  1.1167  for  no  grid,  and  grids  1  through  4  respec¬ 
tively.  These  values  were  used  in  the  CONVERT  program  to  calculate  the 
pressure  and  velocity  in  the  test  section. 

2.  Yaw  and  Pitch  Angle  Corrections 

A  test  was  conducted  to  find  correction  factors  for  yaw  and  pitch 
angle.  From  previous  arrow  plot  data  by  Lung  [Ref.  3],  inconsistent  crossflow 
velocity  magnitudes  and  directions  were  noted  toward  the  outer  boundaries  of 
the  body-only  missile  configuration  run.  These  outer  boundaries  represented 
the  wind  tunnel  freestream  region,  where  crossflow'  velocity  is  expected  to 
reach  zero.  Therefore,  to  duplicate  just  the  freestream  region,  this  prelimi¬ 
nary  test  consisted  of  placing  the  pressure  probe  in  the  tunnel  with  no 
VLSAM  model  and  no  grid.  Thus,  the  expected  pitch  and  yaw  angles  should 
both  have  been  zero.  This  was  not  the  case  however. 

Results  of  the  preliminary  run  (R001A2)  are  listed  in  Appendix  D 
(Result  00.DAT  File).  The  yaw  angles  measured  ranged  from  -4.00°  to  -7.00°, 
with  an  average  of  -5.00°.  The  pitch  angles  (alpha)  ranged  from  approxi¬ 
mately  +  17.4°  to  approximately  +18.4°,  with  the  average  +17.942°.  Thus,  to 
correct  for  these  errors,  +5.00°  was  added  to  the  yaw  angle  and  -17.942°  was 
subtracted  from  the  alpha  values  in  the  output  file  of  the  CONVERT  program 
(RESULT.DAT). 

Though  the  exact  cause  for  the  yaw  and  pitch  angle  errors  was  not 
known,  one  possibility  might  have  been  a  slight  bend  which  was  noted  in  the 
5-hole  pressure  probe.  Other  causes  may  have  been  improper  calibration  of 
the  probe  or  a  misalignment  of  the  traverser  assembly.  The  corrected  errors 
only  effect  the  velocity  vector  plots,  and  should  have  no  effect  on  the  pressure 
contour  plots. 
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III.  RESULTS 


The  following  sections  discuss  the  velocity  vector  plots  and  the  total  and 
static  pressure  coefficient  contour  plots  for  the  VLSAM  model  configurations 
A  (plus)  and  C  (cross),  both  with  and  without  turbulence.  All  plots  depict  the 
3"  by  5.5"  data  acquisition  field  and  the  position  of  the  missile  model  (nose 
aspect)  relative  to  the  field.  Vortex  sizes  and  vortex  distances  from  the  model 
surface  are  referenced  to  the  model  base  diameter  d  (1.75"). 

A.  CONFIGURATION  OA  CTLUS"  WITHOUT  TURBULENCE) 

For  the  plus  configuration,  the  swirling  patterns  of  the  velocity  vector 
plot  (Figure  21)  clearly  denote  the  two  asymmetric  vortices,  which  form 
circles  that  rotate  in  opposite  directions.  Although  the  bottom  vortex  center  is 
evident,  this  plot  fails  to  show  the  center  of  the  top  vortex.  The  vortex 
strength  is  a  maximum  on  the  outer  edges  of  the  vortex  cores,  denoted  by  the 
large  vector  arrows,  where  the  velocities  flow  back  toward  the  missile  body. 
Towards  the  outer  boundaries  of  the  data  acquisition  field,  where  the  vortex 
strengths  are  minimal,  the  vectors  plot  as  points. 

The  total  pressure  coefficient  (Cpt)  contour  plot  (Figure  22)  shows  that 
the  extent  of  the  top  vortex  is  approximately  0.72d  at  a  distance  of  0.33d  from 
the  missile  surface.  The  bottom  vortex  is  0.83d  at  0.67d  from  the  body.  There 
are  more  changes  in  the  pressure  gradient  within  a  smaller  area  for  the  top 
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Figuie  21.  Velocity  Vector  Plot  -  Configuration  OA 
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vortex,  which  varies  in  Cpt  from  -0.14  to  -1.35  within  an  inch,  indicating 
greater  vortex  strength.  The  bottom  vortex  is  slightly  more  diffused.  The  Cpx 
plot  shows  that  the  distance  between  the  two  vortex  centers  is  approximately 
0.77d. 

The  static  pressure  coefficient  (Cps)  contour  plot  (Figure  23)  shows  that 
the  top  vortex  extends  approximately  l.Od  at  a  distance  of  0.49d  from  the 
body.  The  bottom  vortex  extends  about  the  same  ( 1. Id)  at  a  distance  of  0.55d. 
Again,  the  top  vortex  appears  to  be  stronger  than  the  bottom  vortex.  Cps  for 
the  top  vortex  varies  from  -0.55  to  -2.70  within  1.2  inches  while  the  bottom 
vortex  varies  from  -0.55  to  about  -2.18  in  1.6  inches.  The  distance  between  the 
vortices  is  approximately  0.70d  on  the  Cps  plot. 

B.  CONFIGURATION  3A  ("PLUS"  WITH  TURBULENCE) 

With  added  turbulence  (grid  3),  the  velocity  vector  plot  (Figure  24)  still 
indicates  vortex  asymmetry,  but  it  also  indicates  that  the  vortices  have  less 
strength  (smaller  vector  arrows)  than  for  the  no  grid  condition.  Again  the 
vortex  strength  is  maximized  on  the  outer  edges  of  the  vortex  cores. 

From  the  Cpr  contour  plot  (Figure  25),  the  top  vortex  extends  to  approx¬ 
imately  0.72d,  centered  at  about  0.44d  from  the  model.  The  bottom  vortex  is 
slightly  larger  (0.83d)  and  is  located  0.55d  from  the  model  body.  The  pressure 
gradient  of  the  top  vortex  is  steeper  than  for  the  bottom  vortex,  indicating 
greater  strength.  Cpx  varies  from  -0.14  to  -1.5  within  an  inch  for  the  top 
vortex  and  from  -0.14  to  -1.35  within  1.4  inches  for  the  bottom  vortex.  The 
distance  between  the  two  vortices  as  measured  on  the  Cpr  plot  is  approxi¬ 
mately  0.72d. 
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Examination  of  the  Cps  contour  plot  (Figure  26)  reveals  that  both 
vortices  are  roughly  0.83d  wide.  The  top  vortex  is  closer  to  the  missile  body 
(0.39d)  than  the  bottom  vortex  (0.55d).  More  contour  levels  in  a  smaller  area 
and,  consequently,  a  higher  gradient  and  stronger  vortex  exists  for  the  top 
vortex.  Cps  for  the  top  vortex  varies  from  -0.27  to  -1.60  within  0.8  inch  while 
the  bottom  vortex  Cps  ranges  from  -0.27  to  -1.36  within  1.0  inch.  On  the  Cps 
plot,  the  two  vortices  are  about  0.66d  apart. 

C.  CONFIGURATION  0C  ("CROSS"  WITHOUT  TURBULENCE) 

For  the  model  C  configuration  (cross),  the  velocity  vector  plot  (Figure  27) 
still  clearly  displays  the  two  asymmetric  vortices.  As  with  configuration  A, 
the  strengths  of  the  vortices  are  maximized  at  the  edges  of  the  circular  swirls, 
where  the  flow  is  back  toward  the  missile  body. 

Figure  28,  the  total  pressure  coefficient  plot,  displays  a  0.77d  wide  top 
vortex  located  at  about  0.5d  from  the  model  surface.  The  bottom  vortex 
extends  to  approximately  0.88d  and  is  0.66d  from  the  model.  The  bottom 
vortex  is  more  tightly  wrapped  (i.e.,  more  contour  levels  per  area)  closer  to  its 
core  than  is  the  top  vortex.  The  Cpp  range  for  the  bottom  vortex  (-0.14  to 
-1.23)  is  slightly  greater  than  the  range  for  thp  top  vortex  (-0.14  to  -1.09).  The 
vortex  centers  are  roughly  0.66d  apart. 

From  the  static  pressure  coefficient  contour  (Figure  29),  both  vortices 
extend  to  0.95d  with  the  top  vortex  slightly  closer  to  the  model  body  (0.56d) 
than  the  bottom  vortex  (0.66d).  The  relative  strengths  of  the  vortices  is  diffi¬ 
cult  to  interpret  from  the  Cps  plot,  but  the  top  vortex  appears  to  be  a  little 
less  diffused  than  the  bottom  vortex  and  therefore  stronger.  As  with  the  total 
pressure  coefficient  plot,  ti  e  vortices  are  0.66d  apart. 


52 


Y-AXIS  (IN) 


53 


Y-AXIS  (IN) 


-2.0  -1.0  0.0  1.0  2. 

X-AXIS  (IN) 


Figure  27.  Velocity  Vector  Plot  -  Configuration  0C 
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Figure  28.  Total  Pressure  Coefficient  -  Configuration  OC 
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D.  CONFIGURATION  3C  ("CROSS"  WITH  TURBULENCE) 

The  velocity  vector  plot  (Figure  30)  shows  that  the  addition  of  turbulence 
(grid  3)  diffuses  the  asymmetric  vortices  for  model  configuration  C.  The  vector 
arrows  plot  smaller  for  this  turbulent  case  than  for  the  no  grid  case. 

The  Cpr  contour  plot  (Figure  31)  reveals  a  0.92d  bottom  vortex  at  a  dis¬ 
tance  of  0.66d  from  the  model  body  and  a  0.83d  wide  top  vortex  at  0.55d  from 
the  model.  The  bottom  vortex  appears  to  be  slightly  more  tightly  wrapped 
towards  the  vortex  center  than  the  top  vortex.  Cpr  levels  for  both  vortices 
range  from  about  0.0  to  -1.23.  The  distance  between  the  vortices  is  approxi¬ 
mately  0.72d. 

Figure  32,  the  static  pressure  coefficient  plot,  like  the  Cpt  plot,  also 
shows  the  distance  between  the  vortices  to  be  0.72d.  The  top  vortex  (extent 
0.77d)  is  roughly  0.55d  from  the  model's  surface  while  the  bottom  vortex 
(extend  0.88d)  is  0.66d.  The  vortices  appear  to  have  the  same  strength  on  the 
Cps  plot,  in  which  the  coefficients  range  from  0.0  to  approximately  -1.1. 

E.  COMPARISONS 

1.  Between  Body  Configurations  (A  and  C) 

The  following  observations  were  noted  when  making  comparisons 
between  the  two  missile  configurations  for  both  the  turbulent  (grid  3)  and 
non-turbulent  (no  grid)  runs. 

(1)  From  the  total  pressure  coefficient  (Cpr)  contour  plots,  the  vortices 
are  slightly  larger  for  configuration  C.  This  is  more  pronounced  for 
the  run  conducted  with  added  turbulence  (grid  3). 
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Figure  30.  Velocity  Vector  Plot  -  Configuration  3C 
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(2)  The  vortices  of  the  C  configuration  appear  to  be  centered  a  little  fur¬ 
ther  away  from  the  model  surface  than  do  the  A  configuration 
vortices.  This  shift  is  on  the  order  of  0.2  inches,  and  is  apparent  in 
both  the  no-grid  and  grid  3  runs. 

(3)  The  relative  distance  between  vortex  centers  does  not  seem  to  be 
affected  when  the  configurations  are  changed. 

(4)  The  bottom  vortex  is  more  tightly  wrapped  (stronger)  for  configura¬ 
tion  C  based  on  the  Cpt  plots.  This  is  more  noticeable  for  the  run 
conducted  without  turbulence. 

2.  Between  Turbulence  Levels  (0  and  3) 

The  following  observations  were  noted  when  making  comparisons 
between  the  two  turbulence  levels  tested  for  each  body  configuration. 
Turbulence  level  0  was  the  no-grid  condition  and  level  3  was  the  run  con¬ 
ducted  with  a  grid  inserted. 

(1)  The  velocity  vector  plots  indicate  smaller  arrows,  therefore  weaker 
vortices  (less  crossflow),  when  turbulence  is  added.  This  trend  holds 
for  both  body  configurations. 

(2)  The  two  vortices  remain  in  approximately  the  same  positions  as  indi¬ 
cated  on  both  the  vector  plots  and  the  static  pressure  coefficient  (Cps) 
plots). 

(3)  From  the  total  press  coefficient  (Cpr)  plots,  the  center  of  the  bottom 
vortex  for  configuration  A  shifts  slightly  closer  to  the  model  body 
when  the  turbulence  increases.  (-0.2  to  0.0  on  the  x-axis). 

(4)  Also  on  the  Cpt  plots,  with  added  turbulence,  the  top  vortex  core  for 
configuration  C  seems  to  shift  away  from  the  VLSAM  body. 

(5)  The  Cpr  plots  for  both  configurations  are  slightly  more  diffused  with 
added  turbulence.  Specifically,  the  vortex  centers  are  less  tightly 
wrapped,  thus  indicating  weaker  vortices  for  the  more  turbulent 
condition. 

(6)  There  are  no  noticeable  differences  in  the  static  pressure  coefficient 
contours  between  the  turbulent  and  non-turbulent  runs. 
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3.  With  Body-Only  Configuration  (B) 

The  following  observations  were  made  in  order  to  correlate  the 
results  of  this  study  with  the  previous  experiment  conducted  by  Lung  for  a 
body-only  VLSAM  model  configuration.  [Ref.  3]  The  purpose  is  to  note  how 
the  addition  of  wings  and  strakes  might  affect  the  vortex  flow  pattern  around 
the  missile.  The  velocity  vector,  total  pressure  and  static  pressure  coefficient 
plots  for  the  body-only  run  (OB)  are  displayed  in  Figures  33,  34,  and  35 
respectively.  These  comparisons  were  made  for  the  nominal  ambient  (no-grid) 
flowfield  condition. 

(1)  From  the  velocity  vector  plots,  the  vortex  pattern  for  configuration  B 
is  similar  to  the  one  for  body  C  in  that  the  vortex  asymmetry  is  more 
pronounced  for  these  cases.  Body  A  vortices  are  closer  to  the  model 
surface  and  asymmetry  is  less  pronounced. 

(2)  From  the  total  pressure  coefficient  contours,  the  vortices  appear  to  be 
located  much  closer  to  the  missile  body  for  body  A  than  for  either 
body  C  or  B.  Apparently  the  strake/wing-generated  vortices  for  this 
configuration  act  to  hold  the  nose-generated  vortices  near  the  body, 
resulting  in  the  higher  induced  side  forces  for  this  configuration 
observed  by  Rabang. 

(3)  Also  from  the  Cpt  plots,  the  strengths  of  the  bottom  vortices  appear 
to  be  roughly  the  same  for  bodies  A  and  B,  while  the  plot  for  body  C 
seems  to  be  more  tightly  wrapped  (stronger)  near  its  center. 

(4)  From  the  static  pressure  contours,  the  body-only  vortices  are  more 
diffused  and  larger  than  either  configuration  with  wings  (A  and  C). 

(5)  The  relative  distances  between  the  top  and  bottom  vortices  remains 
the  same  on  all  three  plots  for  all  body  configurations. 

In  general,  the  vortex  pattern  around  the  nose  of  the  missile  model  with 
wings  (in  either  configuration)  resembles  the  vortex  pattern  for  a  body-only 
configuration.  Though  there  are  subtle  differences  as  previously  noted,  the 
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relative  strengths,  sizes  and  positions  of  the  asymmetric  vortices  were  com¬ 
parable  in  all  cases. 


Figure  33.  Velocity  Vector  Plot  -  Configuration  OB 
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Figure  34.  Total  Pressure  1 
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-  Configuration  OB 


Figure  35.  Static  Pressure  Coefficient  -  Configuration  OB 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  flowfield  about  a  Vertically-Launched  Surface-to-Air  Missile 
(VLSAM)  model  at  high  angle  of  attack  was  investigated  in  the  wind  tunnel 
of  the  Naval  Postgraduate  School.  Missile  "plus”  and  "cross"  body  configura¬ 
tions  (A  and  C  respectively)  were  both  tested.  The  angle  of  attack  was  set  at 
50°  with  a  Reynolds  number  of  1.1x10s  for  all  data  runs.Two  flowfield  condi¬ 
tions  were  treated:  the  nominal  ambient  wind  tunnel  condition  (no  grid)  and 
a  condition  with  a  grid-generated  turbulence  of  length  scale  1.08  inches  and 
1.88%  turbulence  intensity  (see  Table  2).  The  following  conclusions  were 
reached: 

(1)  An  increase  in  turbulence  intensity  tended  to  reduce  the  strength  of 
the  asymmetric  nose-generated  vortices. 

(2)  The  two  asymmetric  vortices  remained  in  approximately  the  same 
position  for  an  increase  in  turbulence. 

(3)  The  top  vortex  was  closer  to  the  model  surface  and  appeared  to  be 
stronger  foi  both  body  configurations.  This  condition  was  more  pro¬ 
nounced  for  configuration  A. 

(4)  The  wing/strake  arrangement  of  configuration  C  caused  the  vortices 
to  be  centered  further  away  from  the  model  surface  than  those  of  con¬ 
figuration  A,  correlating  with  the  differences  in  induced  side  forces 
for  those  configurations  observed  by  Rabang. 

(5)  Configuration  C  vortices  are  more  diffused  and  larger.  This  was  more 
apparent  when  turbulence  was  added. 

(6)  Crossflow  velocity  vector  plots  agreed  with  the  behavior  denoted  by 
the  total  and  static  pressure  coefficient  contours  for  both  body  config¬ 
urations  and  turbulence  levels. 
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(7)  Though  subtle  differences  exist,  the  addition  of  wings  and  tails  did 
not  greatly  alter  the  vortex  pattern  around  the  nose  of  the  missile 
model,  when  compared  to  the  body-only  configuration  tested  by  Lung. 


Recommendations  for  a  continued  study  of  the  behavior  of  asymmetric 
vortices  under  varying  flowfield  conditions  are  suggested  as  follows: 

(1)  Examine  the  vortices  at  positions  further  back  along  the  model  body, 
where  effects  from  the  wings  might  better  be  seen.  (Such  as  at  a 
length/diameter  ratio  of  9.) 

(2)  Investigate  asymmetric  vortex  behavior  for  just  a  body-only  configu¬ 
ration,  at  the  position  described  above,  to  provide  comparisons. 

(3)  Continue  to  study  vortex  behavior  at  various  angles  of  attack  and 
turbulence  with  varying  intensities  and  length  scales  in  order  to  pro¬ 
vide  a  large  data  base,  which  can  be  used  to  calculate  vorticity 
contours. 
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APPENDIX  A.  PPROBE  PROGRAM 


1  JEF  SEG: CLEAR  , &HFEOO : GOTO  4  'Begin  PCIB  Program  Shell 

2  GOTO  1000  1  User  program 

3  GOTO  900  '  Error  handling 

4  I=&HFE00  '  Copvright  Hewlett-Packard  1984,1985 

5  PCIB. DIRS=ENVIRON$( "PCIB" ) 

6  I*=PCIB.DIRS+"\PCIBILC. BLD" 

7  BLOAD  1* , I 

8  CALL  I  (PCIB.  DIRS,  r'.,J%)  :  PC  I B  .  SEG=  I*/. 

9  IF  J%=0  THEN  GOTO  13 

10  PRINT  "Unable  to  load."; 

11  PRINT  "  (Error 

12  END 

13  ’ 

14  DF.F  SEG  =  PCIB . SEG : O . S=  5 : C . S= 10 :  I  .  V  =  1 5 

15  I ,C=20:L.P=25:LD.FILE=30 

10  GET . MEM- 35 : L . S- 40 : PANELS  =  4  5 : DEF . EPR=  50 

17  PCIB. ERRS  =  STRINGS (  64,32)  :  PCI B . NAMES  =  STRI NGS ( 1 6 . 3 2 ) 

38  CALL  DEF. ERRtPCIB. ERR. PCIB. ERRS, PCIB. NAMES, PCIB. GLBERU)  :  PC  1 1) .  BASERR=  2  5  5 

19  ON  ERROR  GOTO  3 

20  J  =  -  1 

21  I*=PCJB.DIRS+"\PCIB.SYN" 

22  CALL  O.S( IS  ) 

2  3  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

24  1=0 

25  CALL  I  .  V  (  I, READ. REGISTER, READ. SF.LF1D, DEFINE,  INITIALIZE. SVSTEM) 

26  IF  PCIB.ERROO  THEN  FRPOP  CCIR  PASSER 

27  CALL  I ,V( I .ENABLE. SYSTEM, DISABLE. SYSTEM, INITIALIZE. POWER. ON) 

28  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

29  CALL  I . V ( I , MEASURE , OUTPUT . START , HALT ) 

30  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASEFR 

31  CALL  I . V ( I .ENABLE. I  NT . TRIGGER , DISABLE . INT . TRIGGER , ENABLE .  (OUTPUT , DISABLE . OUTPU 
T) 

32  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

33  CALL  I . V ( I , CHECK. DOME. GET. STATUS, SET. FUNCTION. SET .RANGE ) 

34  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

35  CALL  I .V( I , SET. MODE, WRITE. CAL, READ. CAL, STORE. CAL) 

36  IF  FC1B.EFRO0  THEM  ERROR  PCIB. BASERR 

3  7  CALL  I . V( I .DELAY. SAVE. SYSTEM, J, J  ) 

3  8  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

39  1  =  1 

40  CALL  I  .V(  I  ,  SET.  GATF.TIME,  SET.  SAMPLES,  SET.  SLOPE,  SET.  SOURCE! 

41  IF  PCIB.ERROO  THEN  ERROR  PCI  B .  FASERR 

42  CALL  I . C ( I , FREQUENCY , AUTO . FREQ , PERIOD , AUTO . PER ) 

43  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

44  CALL  I ,C( I , INTERVAL, RATIO, TOTALIZE, R100MILLI ) 

4  5  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

46  CALL  I.C( I.R1,R10,R100,R1KILO) 

47  IF  PCIB.ERROO  THEM  ERROR  PCIB. BASERR, 

48  CALL  I . C ( I , R IOMEGA , R100MEGA , CH AN . A , CHAN . B ) 

4  9  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

50  CALL  I .C( I .POSITIVE, NEGATIVE, COMN , SEPARATE) 

51  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

52  1=2 
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5  9  CALL  I  .  V  (  I  ,  ZERO  .  OHMS  .  SET  .  SPEED  .  J  ,  J  ) 

55  IF  PCIB.ERROO  THEN  ERROR  PC I  R  .  BASERR 

56  CALL  I . C ( I , DC VOLTS . AC VOLTS , OHMS , R200MI LLI ) 

5  7  IF  PC1B.ERROO  THEN  ERROR  PCI  B.  BASERR 

58  CALL  I .C( I ,R2,R20,R200, R2KIL0 ) 

59  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

60  CALL  I . C ( I , R20KILO, R200KILO , R2MEGA , R 2 OMEGA ) 

61  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

62  CALL  I . C ( I , AUTOM ,R2.5,R12.5,J) 

6  3  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 
69  1=9 

65  CALL  I . V< I , SET. COMPLEMENT, SET. DRIVER, OUTPUT. NO. WAIT. ENABLE. HANDSHAKE) 

66  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

67  CALL  I . V ( I , DISABLE. HANDSHAKE, SET. THRESHOLD. SET. START. BIT , SET . NUM . BITS ) 

68  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

69  CALL  I . V ( I , SET . LOGIC . SENSE , J , J , J  ) 

70  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

71  CALL  I .C{ I , POSITIVE, NEGATIVE, TWOS, UNSIGNED) 

72  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

73  CALL  I .Cl  I .OC.TTL.RO ,R1 ) 

79  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

75  CALL  I . C ( I ,R2,R3,R9,R5) 

76  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

77  CALL  I . C ( I , R6 , R7 , R8 , R9  ) 

7  8  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

79  CALL  I . C( I .RIO , Rll ,R12 ,R13  ) 

80  IF  PCIB.ERROO  THEN  ERROR  ,  PCI  B  .  BASERR 

81  CALL  I.C(I,R19,R. O.R16.J) 

82  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

83  1-6 

69  CALL  I  V( I, SET. FREQUENCY, SET. AMPLITUDE. SET. OFFSET, SET. SYMMETRY) 

85  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

86  CALL  I ,V( I . SET. B"RST. COUNT, J.U.J) 

87  IF  PCIB.ERROO  THEN  ERROP  PCIB.  BASERR 

89  CALL  I ,C(I , SINE, SOUAPE, TRIANGLE, CONTINUOUS) 

89  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

90  CALL  I . C ( I , GATED , BURST , J , J  ) 

91  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 
<>2  1  =  7 

9  3  CALL  I . V ( I , AUTOSCALE , CALI  PRATE . SET . SENS  I TI V I TY , SET . VERT . OFFSET I 
Q9  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

95  CAI.L  I  ,V(  I  , SET. COUPLING, SET. POLARITY, SET. SWEEPSPEF.D, SET. DELAY) 

96  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

97  CALL  I . V ( I , SET. TRIG. SOURCE, SET. TRIG. SLOPE, SET. TRIG. LEVEL, SET. TRIG. MODE) 

98  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

CALL  I . V ( I , GET. SINGLE. WF . GET. TWO. WF . GET. VERT. 1 NFO , GET . TIMEPASE . INFO' 

100  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

101  CALL  I . V  < I  .GET. TRIG. INF O . CALC . WFV0LT . CALC . WFTIME , CALC . WF . STATS ) 

102  IF  PCIB.ERROO  THEN  ERROR  PCIB.PASEPP 

103  CALL  I . V ( I , CALC. R I SET I ME, CALC . FALl.TIME , CALC . PERIOD , CALC . FREQUENCY  ) 

109  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

105  CALL  I ,V( I , CALC. PLUSWIDTH, CALC. MINUSWIDTH. CALC. OVERSHOOT. CALC. PRESHOOT) 

106  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

107  CALL  I . V ( I , CALC. PK . TO. PK , SET. TIMEOUT, SCOPE. START, MEASURE. SINGLE. WF  ) 

108  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

109  CALL  I . VI  I , MEASURE . TWO . WF . J ,0 , J ) 
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110  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERP 

11  1  CALL  J  .  C  (  1  ,  R  ION  AMU  ,  RIOONANO  .  R1M1CRO  ,  RIOMICL’O  ) 

112  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

113  CALL  I . C ( I , R1 00MI CRO .R1MILL1 , R10H1LLI . R100HI LLI  ) 

114  IF  PCIB.ERROO  THEM  ERROR  PC  I B  .  BASF.RR 

115  CALL  I.C(I,R1,R10, R20NANO , R200MANO ) 

116  IF  PCIB.ERROO  THEM  ERROR  PCIB. BASERR 

11”  CALL  I . C ( I , R2MICRO , R20MICRO , R200MI CRO , R2MI LLI ) 

118  IF  PCIB.ERROO  THEM  ERROR  PCIB. BASERR 

119  CALL  I  .  C  ( I ,  R20MI  LLI  .  R200MILLI  ,  R2  ,  P.20  ) 

120  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

121  CAT  L  I . C ( I , R50NANO , R500MANG , K5MICR0 . R50MICRO ) 

122  IF  PCIP..ERROO  THEN  ERROR  PCIB.  BASERR 

123  CALL  I . C ( I , R500MI CRO , R5MILLI , R50MI LLI .R500MILLI 1 

124  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

125  CALL  I.C(I,R5,R50,CHAN.At  CHAN . B ) 

126  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

127  CALL  I . C f I , EXTERNAL, POSITIVE , NEGATIVE , AC ) 

128  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

129  CALI.  I ,C( I , DC.TRIGGEPED, AUTO. TRIG , AUTO. LEVEL ) 

130  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

131  CALL  I ,C( I .XI ,X10. STANDARD, AVERAGE) 

132  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR! 

133  1  =  8 

134  CALL  I  . V( 1 .OPEN. CHANNEL. CLOSE. CHANNEL. J . J  ) 

135  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

1 36  CALL  C.S 

137  IF  PCIB.ERROO  THEM  ERROR  PCIB.BASERP 

138  1 1  =  FCI  E? .  DI  Rf  *  " \PCI  B  .  PLD" 

139  CALL  L . P ( IJ  ) 

140  IF  PC  ID.ERRO0  THEN  ERROR  PCIB.  BASERR 
1  4  1  1  $  =  "DM1'  .01"  :I=3:J=0:K=0:L=1 

142  CALL  DEFINE ( DMM . 0 1 . I t , 1 , 1 , K , L ) 

143  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASF.RR 

144  I  $  =  "Func . Gen . 0 1 " : I~6: J  =  0:K=1 :L~1  1 

145  CALL  DEFINE ( Func . Gen . 01 , IS , I . J . K , L ) 

146  IF  PCIB.EPP.OO  THEN  ERROR  PCiB. BASERR 
14 /  I  $ -“Scope. 01" : I  -  7 : J  =  0 : K  =  2 : L= 1 

148  CALL  DEFINE!  Scope  .  01  ,  !  $  .  I  .  J  , K  .  1. ) 

149  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASF.RR 

150  1$= "Counter. 01": 1  =  1 : J  =  0 : K  =  3 : L= 1 

151  CALL  DEFINE ( Counter . 0 1 , 1 J , I , J , K , L ) 

152  IF  PCIB.ERROO  THFN  EFPOR  PCIB.  BASERR 

153  I  $  =  "Dig .  In. 01":  1  =  '*  :J  =  0:K  =  4  :  L=  1 

154  CALL  DEFINE! Dig . In . 01 . I t , I . J ,K, L ) 

155  IF  PCIB.ERROO  THEN  ERROR  P'"  J  p  .  PASEFR 
l^e  I  $  =  "Dig  .  Out  .01":  I  =  4  :  .3  -  1 :  K  =  •  .  L  =  1 

3  57  CAL?,  DEFINE!  Dig  .out .  01 . 3 1  .  I  ,  J  ,K  ,  L  ) 

158  IF  PCIB.ERROO  THEN  ERROR  PC  I  b  .  BASERR 

159  It- "Relay .Act . 0 1 " : I =8 : J = 0 : K= 5 : L  1 

160  CALL  DEFINE! Relay . Act . 01 , I  $ , I , J , V , L ) 

161  IF  PCIB.  ERR.'.  0  THEN  ERROR  PC  I  B  .  P  ASF.R.R 

162  I  $  =  " Re  1 ay . Mux . 0 1 " : I  =  2  :  J  ~  0  :  K  =  6  :  L  =  1 

163  CALL  DEFINE! Felay .Mux. 01 , 1 1 , I , J ,K ,1.  ) 

164  IF  PCIB.ERROO  THFN  ERROR  RCIB. BASERR 
8  00  1*  =  F.NV IKONS  (  "PANELS"  )  +  ”\PANELS  .  EXF." 

801  CALL  L . S ( I  * ) 

899  GOTO  2 
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900  IF  ERR = PC IB . BASERR  THEN  GOTO  903 

901  PRINT  "BASIC  error  It "  ;  ERR ;  "  occurred  in  line  "  :  EP.L 

902  STOP 

903  TMPERR=PC1B. ERR: IF  THPEilR  -  C  THEN  1MFERR=PCI F . GLBERR 

904  PRINT  "PC  Instrument  error  f";TMPERR;"  detected  at  line  ";ERL 

905  PRINT  "Error:  ”;PCIB.ERR$ 

906  IF  LEFTS ( PC IB . NAMES , 1 ) <  >CHR$ ( 32 )  THEN  PRINT  "Instrument:  "; PCI B . NAMES 

907  STOP 

908  COMMON  PCI B . DIRS , PCI B . SEG 

909  COMMON  LD . FI LE , GET . MEM , PANELS , DEF . ERR 

910  COMMON  PC IB . BASERP . PC i h . FRR  pc  £ B . ERRS , PC IB. NAMES , PC  I P . ^LPERP 

911  COMMON  READ. REGISTER .READ. SELF  TO, DEFINE, INITIALI ZE . SYSTEM , ENABLE . SYSTEM . DI SA 
BI.E. SYSTEM. I NI T I AL I ZE . POWER . ON , ME ASURE . OUTPUT , START . HALT , ENABLE . I NT . TR I GGER . DI SA 
BLE . INT . TRIGGER . ENABLE . OUTPUT , DISABLE . OUTPUT , CHECK . DONE , GET . STATUS 

912  COMMON  SET . FUNCTION . SET . RANGE , SET . MODE . WRITE . CAL , READ . CAL , STORE . CAL , DELAY , SA 
VE . SYSTEM , SET . GATETIME . SET . SAMPLES , SET . SLOPE . SET . SOURCE , ZERO . OHMS . SET . SPEED , SET . 
COMPLEMENT . SET . DRIVER , OUTPUT . NO . WA IT , ENABLE . HANDSHAKE , D I  SABLE . HANDSHAKE 

913  COMMON  SET . THRESHOLD . SET . START .BIT. SET . NUM . BITS , SET . LOGIC . SENSE . SET . FREQUENC 
Y . SET . AMPLITUDE , SET . OFFSET , SET . SYMMETRY . SET . BURST . COUNT . AUTOSCALE . CALIBRATE . SET . 
SENS ITI V I TY , SET . VERT . OFFSET , SET . COUPL I NG , SET . POLAR  I TV , SET . SWEEPSPEED 

9 1  4  COMMON  SET  .  DELAY  .  SET  .  TR  I G  .  SOURCE  .  SET  .  TRI G  .  SLOPE  ,  SET  .  TR  I G  .  LEVEL  .  SET  .  TP.  I G  .  MODE 
, GET . SINGLE . WF , GET . TWO . WF , GET . VERT . INFO, GET . TIMEBASE . INFO . GFT .TRIG. INFO , CALC . WFV 
OLT , CALC . WFTIME , CALC . WF . STATS , CALC .RISETIME , CALC . FALLTIME . CALC . PERIOD 

9 1 5  COMMON  CALC . FREQUENCY , CALC . PLUSWIDTH , CALC .MINUSWIDTH . CALC . OVERSHOOT . CALC . PRE 
SHOOT , CALC . PK . TO . PK , S  ET . TIMEOUT . SCOPE . START , MEASURE . S I NGLE . WF , MEASURE . TWO . WF . OPE 
N . CHANNEL , CLOSE . CHANNEL 

916  COMMON  FREQUENCY, AUTO. FREQ, PERIOD , AUTO . PEP , INTERVAL. RATIO, T07ALI ZE, R100MILLI 
, R1 , RIO, RIOO.RIKILO.R IOMEGA , R1 OOMEGA .CHAN.A.CHAN.B, POSITIVE .NEGATIVE, COMN , SEPARA 
TE . DC VOLTS , AC VOLTS , OHMS .R200MILLI , R2 ,  P.2  0 , R200 , R2KILO , R20KILO , R200KI LO 

9  17  C'_**MON  R2MEGA  .  R2  OMEGA  ,  AUTOM,  R2 . 5  ,R12 . 5  ,  POSITI  VE  .  NEGATI  VE  ,  TWOS  ,  UNSI GNED  ,  OC  .  TT 
L,RO,R:R2,R3.R4,R5.R6.R7,R8,R9,RjO,Rir.Ri2,R13 , R14 , R1 5 , R 1 6 , SI NE , SQUARE . TRI ANGLE 
, CONTINUOUS, GATED, BURST, RIONANO .  P.IOONANO , R1MICRO , RIOMICRO , P 100MI CRO 

918  COMMON  R1MILLI , R10MILLI , R1O0MILLI , R1 ,R10 , R20NANC* , R200NANO , R2M I CPO , R20MICRO , R 
2O0MICRO.R2MILLI , R20MILLI .P200MILLI , R2 . R20 . R50NANO , R500NANO , R5MICRO , R50MICRO. R50 
OMICRO.R5MILLI .R50MILLI .RSOOMILLI , R5 ,  P.50 . CHAN . A . CHAN . B , EXTERNAL , FOS IT  I VE 
9 1  Q  COMMON  NEGATIVE  ,  AC  ,  PC  ,  TR  I  GGEP.ED .  AUTO  .  TRI  G  .  AUTO  .  LEVEL ,  X 1  ,  X 1  b  ,  STANDARD  ,  AVERAGE 
920  COMMON  DMM  .  01  ,  Func  .  Gen  .  01  .  Scope  .  01  .  Counter  .  01  ,  Dip, .  1  n  .  0 1  ,  Di g  .  Ou t  .  0 1  .  Rel  ay .  Act 
. 0 1 , Re  1  ay . Mux . 0 1 

999  'End  PCIB  Program  Shell 

1000  REM  This  step  initialzcs  the  HP  system 
1010  CLS 

1020  OPTION  EASE  1 

10  30  DIM  P( 5 ) , FA<  50 , 5 ) . FF( 50 , 5 ) ,XPT( 50 ) , YFT( 50 ) , X ( 50 ) , Y ( 50  )  , YAW l 50 ) 

1040  REM 

1050  CALL  INITIALIZE. SYSTEM ( PGM3HEL . HFC > 

1060  PEM 

1070  PEM  SET  FUNCTI ONON  THE  'DMM  ,  'RELAY  MUX  ,  'RELAY  ACTUATOR’ 

1080  PEM 

10  90  CALL  SET . FUNCTI OIK  DMM. 01 , DCVOLTS ) 

1100  CALL  SET . RANGE ( DMM .01, AUTOM ) 

1110  CALL  DISABLE.  INT.  TP  I CC-ER  (  DMM  .  01  ) 

1120  CALL  ENABLE. OUTPUT (PELAY .MUX. 01 ) 

1130  call  ENABLE. OUTPUT! RELAY. ACT. 01 ) 

1140  .EM  *************  PROGRAM  TRAVERSE  ******************* 

1150  PEM 

1160  REM  OPEN  THE  c-M  port  AND  INITIALIZE  THE  MOTOR  SETTINGS 

1170  OPEN  "corol: 1200.  1 . r s , cs , ds , cd"  AS  II 

1190  FEM  SET  MOTOR  DI  7  VALUES 

1190  DATA  2000,2000,2!  ,2,2,2,0.000125,0.000125.0.000125  . 
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1200  READ  V1.V2.V3.R1 ,R2,R3,C1.C2,C3 

1210  REM  DEFINE  CHARACTERS  FOR  DATA  REDUCTION  ALGORITHM 
1220  RN2$  =  "RENAME  A : RAW . DAT  " 

1230  HEAD1S  =  "  #  X  Y  PI  P2  P3  P4  P5  Y  A 

W  " 

12  40  FORMAT*  =  "##  ##.##  ##.##  I##.###  ###.*##  *##.##11  ###.«##  ###.###  ### 

.  ##" 

1250  PRINT 

1260  PRINT  "*************t***********t***t******t**t**»" 

1270  PRINT  "**  USER  MUST  SELECT  'CAPS  LOCK'  FUNCTION  ♦+" 

1280  PRINT  "a******************************************" 

1290  REM  DISPLAY  MOTOR  DEFAULT  SETTINGS 

1300  PRINT  "  ******♦*♦****♦*******♦**♦*******+*+» 

1310  PRINT  "  INITIALIZED  VALUES  FOR  ALL  MOTOR  SETTINGS:” 

1320  PRINT  "  VELOCITY  =  1000  STEPS/SEC" 

1330  PRINT  "  RAMPIMOTOR  ACCELERATION)  =  2  (6000  STEPS/SEC ‘ 2 > " 

1340  PRINT  "  DEFAULT  INCREMENTAL  UNITS  ARE  INCHES" 

1350  PRINT  "  ******************************»*»»».. 

1360  PRINT 

1370  PRINT  "NOTE I !  USE  MANUAL  CONTROL  TO  INITIALIZE  PROBE  POSITION  BEFORE” 

1380  PRINT  "  SELECTING  COMPUTER  CONTROLLED  MOVEMENT. 

1390  PRINT 

1400  INPUT  "MANUAL  CONTROL  OR  COMPUTER  CONTROL  (ENTER  'MAN'  or  'CP'!":CON* 

1410  IF  CON*  =  "CP"  THEN  3490 

1420  REM  OPTION  TO  CHANGE  DEFAULT  SETTINGS  OF  VELOCITY  OR  ACCELERATION  RAMP 
1430  PRINT 
1440  PRINT 

1450  PRINT  "  DO  YOU  WANT  TO  CHANGE  THE  VELOCITY  CR  ACCELERATION  RAMP" 

1460  PRINT  "  DEFAULT  SETTINGS?  (Y  or  N)" 

1470  PRINT 

1480  PRINT  "IF  'NO',  THIS  PROGRAM  WILL  THEN  LET  YOU  DEFINE  THE" 

1490  PRINT  "DISTANCE  YOU  WANT  TO  MOVE  (IN  INCHES).  IF  'YES'," 

1500  PRINT  "YOU  CAN  CHANGE  ANY  OR  ALL  OF  THE  DEFAULT  SETTINGS  FOR  ANY  MOTOR." 
1510  PRINT 

1520  PRINT  : 

1530  PRINT 

1540  INPUT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THE  DEFAULT  SETTINGS?  (Y  or  N ) " ; D* 

1  550  IF  D*  =  "Y"  THEN  1590 
1560  IF  D*  =  "N"  THEN  2220 
1570  REM 

1580  REM  *♦**  OPERATOR  SELECTED  MOTOR  VARIABLES  ***** 

1590  PRINT 
1600  PRINT 

1610  INPUT  "WHICH  DEFAULT  VALUE?  (ENTER  ' 1 ' FOR  VELOC  OR  '2'  FOP  ACCEL  RAMP ) " ; L 
1620  ON  L  GOTO  1690, 1V30 

1630  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  DEFAULT  VELOCITY?  (Y  OR  N)" 

1640  INPUT  V* 

1650  IF  V*= "Y"  THEN  1690 

1660  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  DEFAULT  ACCELERATION  RAMP?  (Y  or  N)" 

1670  IF  R* - "Y"  THEN  1990 
1680  IF  R* -  "tl"  THEN  1450 
1690  PRINT 
1700  PRINT 

1710  INPUT  "WHICH  MOTOR  VELOCITY  DO  YOU  WISH  TO  CHANGE?  (1,2,  or  j)";J 
1720  ON  J  GOTO  1730,1830,1880 
1730  PRINT 
1740  PRINT 

1750  INPUT  "ENTER  DESIRED  VELOCITY  O^  MOTOR  #1";V1 
1760  PRINT 
1770  PRINT 
1780  TRINT 

1790  PRINT  "DO  YOU  WANT  TO  CHANGE  VELOCITY  OF  ANOTHER  MOTOR?  (Y  OR  N ) " 
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1800  INPUT  VS 

1810  IF  VS  =  " Y "  THEN  1690 

1820  IF  V$="N"  THEN  1430 

1830  PRINT 

18<i0  PRINT 

1850  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  2";V2 

1860  PRINT 

1870  GOTO  1780 

1880  PRINT 

1890  PRINT 

1900  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  »3";V3 

1910  PRINT 

1920  GOTO  1780 

1930  PRINT 

1940  PRINT 

1950  INPUT  "WHICH  MOTOR  ACCEL  RAMP  DO  YOU  WANT  TO  CHANGE?  (1,  2,  or  3)";K 
1960  ON  K  GOTO  1970,2060,2120 
1970  PRINT 
1980  PRINT 

1990  INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  *1";R1 
2000  PRINT 
2010  PRINT 

2020  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  ACCEL  RAMP  OF  ANOTHER  MOTOR?  ( Y  or  N 

2030  INPUT  RMS 

2040  IF  RMS  =  "Y"  THEN  1930 

2050  IF  RMS  =  "N"  THEN  14  50 

2060  PRINT 

2070  PRINT 

2080  INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  *2";R2 

2090  PRINT 

2100  PRINT 

2110  GOTO  2000 

2120  PRINT 

2130  PRINT 

2140  INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  «3”:R3  : 

2150  PRINT 
2160  PRINT 

2170  GOTO  2000  . 

2180  REM 

2190  REM  DEFINE  DISTANCE  TO  MOVE  MOTOR 
2200  PRINT 
2210  PRINT 
2220  PRINT 

2230  REM  INITIALIZE  MOTOR  INCREMENTS  TO  ZERO 
2240  11-0 
2250  12-0 
2260  13-0 
2270  PRINT 
2280  TRIMT  ' 

2290  .'Hr 
2300  TP 
2  315  - R i ! .  T  ' 

2320  PRINT  ' 

2330  TRINT  ' 

2340  PRINT  ’ 

2350  PRINT  ' 

2360  PRINT  ' 

2370  PRINT  ’ 

2380  PRINT 
2390  PRINT 
2400  INPUT  1 

2410  ON  L  GOTO  2420,2680,2970 


*  ♦ 
*  * 
*  * 
*  * 
*  * 
♦  * 
*  * 
*  * 


DEFINE  WHICH  MOTOR  YOU  WANT  TO  MOVE 


Note  I  I  I  A  POSITIVE  ('  +  ')  INCREMENT  TO  A  MOTOR. 
MOVES  TRAVERSER  AWAY  FROM  THAT  PARTICULAR  MOTOR 


MOTOR 

MOTOR 

MOTOR 


HI  MOVES  7 HE  PROBE  UPSTREAM  AGAINST  THE  FLOW 
12  MOVES  THE  PROBE  TOWARD  THE  ACCESS  WINDOW 
»3  MOVES  THE  PROBE  VERTICALLY  DOWNWARD 


*  t 
1 1 

*  4 

*  4 

*  » 
*  ♦ 
*  ♦ 
*  * 


*»44444*4**4*4-444444444»*******4*4********»4*4*4*444*4444**»» 


'WHICH  MOTOR  DO  YOU  WANT  iO  MOVE?  (1,2,  or  3  ) " ; L 
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2420  PRINT 
2430  PRINT 

2440  PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #1?" 

2450  PRINT  "  *********  ( ENTER  DISTANCE  IN  INCHES)  *»**»****'• 

2460  INPUT  II 
2470  PRINT 

2480  PRINT”  *********************************" 

2490  PRINT 

2500  PRINT  "SUMMARY  OF  OPERATOR  INPUTS:" 

2510  PRINT  "  MOTOR  II  VELOCITY  =  ":V1 

2520  PRINT  "  ACCELERATION  RAMP  =  ";R1 

2530  PRINT  "  INCREMENTAL  DISTANCE  =  1 1 ;" INCHES " 

2540  PRINT"  *********************************" 

2550  PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or  N)" 

2560  PRINT 

2570  PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.  ENTER  'Y'  TO  RETURN" 
2580  PRINT  "TO  VARIABLE  SELECTION  SUBROUTINE." 

2590  INPUT  VJ 

2600  IF  VJ="Y"  THEN  1430 

2610  GOSUB  3410 

2620  PRINT 

2630  PRINT  "DO  YOU  WANT  TO  MOVE  ANOTHER  MOTOR  ALSO?  (Y  or  N)?" 

2640  INPUT  CJ 

2650  IF  CJ="Y"  THEN  2220 

2660  IF  CJ="N"  THEN  3260 

2670  PRINT 

2680  PRINT 

2690  PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  12?" 

2700  PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  *********" 

2710  INPUT  12 
2720  PRINT 
2730  PRINT 

2740  REM  DISPLAY  OPERATOR  SELECTED  MOTOR  VARIABLES 
2750  PRINT"  *********************************" 

2760  PRINT 

2770  PRINT  "SUMMARY  OF  OPERATOR  INPUTS:" 

2780  PRINT  "  MOTOR  12  VELOCITY  =  ";V2 

2790  PRINT  "  ACCELERATION  RAMP  =  ";R2‘ 

2800  PRINT  "  INCREMENTAL  DISTANCE  =  " ; 1 2 ; " INCHES" 

2 B 10  PRINT"  *********************************" 

2820  PRINT 
2830  PRINT 

2840  PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or  N)" 

2850  PRINT 

2860  PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.  ENTER  ’ Y ’  TO  RETURN" 
2870  PRINT  "TO  VARIABLE  SELECTION  SUBROUTINE." 

2880  INPUT  VJ 

2890  IF  Vt= "Y"  THEN  1430 

2900  GOSUB  3410 

2910  PRINT 

2920  PRINT  "DO  YOU  WANT  TO  MOVE  ANOTHER  MOTOR  ALSO?  (Y  or  N)?" 

2 Q 3 0  INPUT  CJ 

2940  IF  CJ="Y"  THEN  2220 

2950  IF  CJ  =  "N"  THEN  3260 

2960  PRINT 

2970  PRINT 

2980  PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  13?" 

2990  PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  *********" 

3000  INPUT  13 
3010  PRINT 
3020  PRINT 

3030  REM  DISPLAY  OPERATOR  SELECTED  MOTOR  VARIABLES 


74 


3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
35  50 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 


PRINT"  ***********»************»********" 

PRINT 

PRINT  "SUMMARY  OF  OPERATOR  INPUTS:" 

PRINT  ”  MOTOR  It  3  VELOCITY  =  ";V3 

PRINT  "  ACCELERATION  RAMP  =  ";R3 

PRINT  "  INCREMENTAL  DISTANCE  =  " ; I  3 ; "INCHES" 

PRINT 

PRINT"  *********************************" 

PRINT 

PRINT 

PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  ( Y  or  N ) " 

PRINT 

PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.  ENTER  *Y'  TO  RETURN" 

PRINT  "TO  VARIABLE  SELECTION  SUBROUTINE." 

INPUT  V* 

IF  V$="Y"  THEN  1430 

GOSUB  3410 

PRINT 

PRINT  1 

INPUT  "DO  YOU  WANT  TO  INPUT  ANOTHER  MANUAL  MOTOR  MOVEMENT  (Y  or  N)";MS 
IF  M$="Y"  THEN  2210 
PRINT 

PRINT  "DO  YOU  WANT  TO  INPUT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT?" 

PRINT  "  *********  NOTE!!!  *********  " 

PRINT  "  ALL  PREVIOUS  MOTOR  INCREMENT  INPUTS  HAVE  BEEN  ZEROIZED." 

PRINT  "PROGAM  WILL  LET  YOU  CHOOSE  MANUAL  OR  CP  -  CONTROLLED  MOVEMENT.” 
PRINT  "***»■*  (IF  'NO'.  THE  PROGRAM  WILL  END).  ****♦" 

PRINT 

INPUT  "DO  YOU  WANT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT  (Y  or  N)";N$ 

IF  Ni = "Y"  THEN  3500 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
END 

REM  *******  MOTOR  MOVEMENT  SUBROUTINE  ********* 

PRINT  *1.  : PR I  NT  HI.  "E" : "Cl  =  " ; Cl ; " :  C2=  " ; C2 ; " : C3  =  " : C3 

PRINT  *1.  "1 1  =  ” : 1 1 ; ” : Vl= " ; VI ; " :R1  =  " ;R1 ; 

PRINT  #1,  "  :  1 2  =  "  ;  1 2  :  ”  :  V2  =  "  :  V2  :  "  :  R2  =  "  : P.2 
PRINT  #1.  " I  3  =  " ; 1 3 ; " : V3  =  " ; V3 ; " : R3= " : R3 ; " : 

RETURN 

REM  ******************************************** 

REM  ******************************************** 

PRINT 

REM  *******  COMFUTER  CONTROLLED  MOVEMENT  ******* 

PRINT 

PRINT  "THE  PRESSURE  DATA  WILL  BE  WRITTEN  TO  FILES  ON  DRIVE  'A' 

PRINT 

PRINT  "YOU  WILL  BE  ASKED  TO  INPUT  FILE  NAMES  FOR  THESE." 

PRINT 
INPUT 
PkII’T 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
FRINT 

REM  SET  INITIAL  MOVEMENT  DISTANCE  AND  NUMBER  OF  DATA  POINTS  TO  ZERO 
HT”  0 


THE  PROGRAM  HAS  ENDED." 


’IS  A  FORMATTED  DISK  IN  DRIVE  'A’?  PRESS  'ENTER'  TO  CONTINUE";D$ 


♦  ♦♦♦a*********************************:**'’ 

**  NOTE  II!  **" 
**  COMPUTER  CONTROLLED  MOVEMENT  **" 
**  IS  PROGRAMMED  WITH  A  **" 
**  DEFAULTED  NEGATIVE  MOTOR  INCREMENT  **" 
**  (i.e.  MOTOR  13  WILL  MOVE  UPWARD  **" 
**  BY  ENTERING  A  (+)  DISTANCE).  **" 
*********************************»******M 
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3710  WD=  0 
3720  DIST-0 
3730  XPT  =  0 

3  7  40  YPT  =  0 
3750  N  =  0 
3760  PRINT 
3770  PRINT 

3780  INPUT  "WHAT  IS  THE  DIMENSION  (  X  ,  Y  )  ( IN  INCHES)  THAT  YOU  WANT  TO  MEASURE 
;WD,HT 
3790  PRINT 

3800  INPUT  "WHAT  IS  THE  STEP  (IN  INCHES)  THAT  YOU  WANT  TO  MOVE.";DIST 

3810  YPT=INT!HT  /DIST)  +  1 

3820  XPT= I NT ( WD  /DIST)+  1 

3830  N=XPT* YPT 

3840  PR I NT 

3850  PRINT  "THERE  ARE  ";XFT:"  *  ";YPT;"  =  ";N;"  POINTS  TO  BE  MEASURED  " 

3860  PRINT 

3870  INPUT  "ARE  THE  NUMBER  OF  POINTS  IS  OK.(Y  OR  N ) " ; CS 
3880  IF  C$="N"  THEN  3780 
3890  CLS 
3900  N=XPT 

3910  IF  (N  <  1)  OR  ( N  >  99)  GOTO  3780 

3920  REM  **♦  GENERATING  STRING  STRING  SEGMENTS  FOR  DATA  FILE  NAMES 
3930  BS  =  MI D$  ( STR$ ( 1  )  ,  2):  REM  **  STRING  NUMBER  "1" 

3980  ES  =  MID$(STR$(N) ,  2):  REM  **  ENDING  STRING  NUMBER  ”N" 

3950  XS  =  "XXXXXX" 

3960  EXS  =  ".DAT" 

3970  CLS 

3980  PRINT  "DATA  FILES  WILL  BE  INCREMENTED  FROM:" 

3990  PRINT 

6000  PRINT  (X*  +  BJ  +  EXJ  )  :  "  To  <X$  +  E$  *  EXS) 

6010  PRINT 
6020  PRINT 

6030  INPUT  "ENTER  DATA  FILE  NAME  (6  CHARACTERS  MAX  --  NO  EXTENSION )": F2S 
6060  PRINT 
6050  PRINT 

6060  IF  LEN( F2S )  >  6  OR  LEN (  F2S  )  <  1  GOTO  8030 
6070  CLS 

6080  PR  I  NT  N;  "DATA  FILES  WILL  BE  GENERATED  AND  INCREMENTED1  AS  FOLLOWS:" 

6090  PRINT 
4100  PRINT 

4110  PRINT  (F2S  +  BS  +  EXS);  "  To  ";  ( F2S  *  EJ  *  EXS) 

4120  PRINT 
4130  PRINT 

4140  INPUT  "ARE  THE  NUMBER  OF  POINTS  AND  FILE  NAMES  OK . ( Y  OR  N)";  CS 

4150  IF  CS  =  "N"  GOTO  37RQ 

4160  IF  CS  =  "Y"  GOTO  4180 

4170  GOTO  4140 

4180  CLS 

4190  PRINT 

4200  PRINT 

4210  REM  SET  INITIAL  POSITION  DATA 

4220  X( 1 ) =-DIST 

4230  Y  (  1  )  =  -  DIST 

4240  FOR  IX=  2  TO  XFT+1 

4250  X ( IX ) =0 

4260  NEXT  IX 

4270  FOR  JY  =  2  TO  YPT*-1 

4280  Y ( JY  )  =0 

4290  NEXT  JY 

4300  FOR  1=1  TO  XPT 

4302  11=0 

4304  12=0 

4306  13=0 

4310  FOR  J=1  TO  YPT 

4  3  20  REM  MOTOR  CP  -  CONTROLLED  MOTOR  MOVEMENT 
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(.330  11  =  0 
4340  12-0 
4350  13=0 

4360  REM  EACH  POINT  TAKE  10  TIMES  READINGS 
4370  X( 1+1 >=X( I ) +DIST 
4380  XPT ( J ) =X ( 1+1) 

4390  Y( J+l )=Y( J ) +DIST 
4400  YPT ( J ) = Y ( J+l ) 

4405  INPUT  "  ADJUST  THE  WHEEL  TO  MAKE  THE  P2  =  P3 , INPUT’ THE  YAW  ANGLE" ; YAW< J  ) 

4408  PRINT 

4410  INPUT  "  PRESS  'ENTER'  TO  START  THE  MEASUREMENT" ; MOVE* 

4420  REM 

4430  REM  READ  FIVE  CHANNELS  AND  DISPLAY  THE  DATA 

4440  REM 

4450  ST£PPER=4 

4460  SWITCH  =  3 

4470  HOMER=8 

4480  DELAY 1  =  .1 

4490  DELAY2  =  1 

4500  REM  SET  THE  S.V  PORT  TO  14 
4510  FOR  I L= 1  TO  3 
4520  THYME  =  TIMER 

4530  CALL  OUTPUT ( RELAY . ACT . 0 1 , STEPPER ) 

4540  CHKTIME  =  TIMER 

4550  IF  CHKTIME  <  (THYME  +  DELAY 1  )  GOTO  4540 
4560  CALL  OPEN. CHANNEL ( RELAY . ACT. 01 , SWITCH ) 

4570  CLS 
4580  NEXT  IL 
4590  PRINT 

4600  PRINT  "  NOW  IS  POINT  ";J 

4610  REM  START  MEASURE  FROM  PORT  4  TO  PORT  8 
4620  FOR  JJ=1  TO  5 

4630  CALL  OUTPUT (RELAY. ACT. 01 .STEPPER) 

4640  CHKTIME  =  TIMER 

4650  IF  CHKTIME  <  (THYME  +  DELAY2 )  GOTO  4640 
4660  REM  EACH  PORT  SAMPLE  10  TIMES 
4670  FOR  11=1  TO  10 
4680  ROUT= 1 

4690  CALL  OUTPUT (RELAY. MUX. 01, ROUT) 

4700  CALL  MEASURE (DMM. 01, VOLTS)  .  , 

47 1U  PA ( II ,JJ  )=VOLTS 

4720  NEXT  II  I 

4730  CALL  OFEN . CHANNEL) RELAY . ACT . 01 . SWITCH ) 

4740  IF  JJ=  5  THEN  4760 
4750  NEXT  JJ 

4760  REM  HOME  THE  S.V  PORT  TO  *48 
4770  CALL  OUTTUT ( PELAY . ACT . 0 1 . HOMER ) 

4780  CALL  OPEN . CHANNEL ( RELAY . ACT . 0 1 , HOMER ) 

4  7  90  RF.M 

4800  REM  DISPLAY  THE  SAMPLE  DATA 
4810  REM 

4820  PRINT  HEAD1S 
4830  FOR  IS=  1  TO  10 

4840  PRINT  USING  FORMATS ; IS , XPT( J ) . YPT( J ) . PA( IS . 1 ) . PA < IS , 2 ) . PA(  IS , 3 ) , PA ( IS , 4 ) , PA 
( IS , 5  )  , YAW ( J ) 

4850  NEXT  IS 
4860  REM 

4870  REM  AVERAGE  THE  DATA 
4880  REM 

4890  FOR  JA  =  1  TO  5 

4900  TOTAL  =  J 

4910  FOR  IA  =  1  TO  10 

4920  TOTAL  =  TOTAL  ♦  PA(IA.JA) 

4930  NEXT  IA 

4940  AVERAGE  =  TOTAL  /10 
4950  F ( JA ) = AVERAGE 
4960  NEXT  JA 
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4970  PRINT 

4980  PRINT  "THE  AVERAGES  ARE:  ” 

5000  PRINT  HEAD1S 
5010  FOR  JD= 1  TO  5 
5020  PP(J,JD)=P< JD) 

5030  NEXT  JO 

3040  PRINT  USING  FORMATS ; J , XPT ( J ) . YPT ( J ) , PP ( J , 1 ) . PP ( J , 2 ) , PP ( J , 3 > . PP ( J . 4 ) , PP t J . 5 ) 

, YAW( J ) 

5045  PRINT 

5048  PRINT  USING  "THE  NULLING  ERROR  IS  +» . 1 1 1 1 " ; PP< J , 3 ) -PP( J , 2 > 

5049  PRINT 

5050  PRINT  "DO  YOU  WANT  RE-MEASURE  AGAIN  (Y  /  N)" 

5060  PRINT 

5062  PRINT  "IF  'Y*  WILL  RE-SAMPLE  AGAIN." 

5064  PRINT 

5070  INPUT  "IF  'N'  WILL  MOVE  THE  TRAVERSER  STEP  UPWARD  (WAIT  7  SEC  )";CS 
5075  PRINT 

5080  IF  C$="Y"  THEN  4405 
5082  IF  CS="N"  THEN  5090 
5084  GO  TO  5070 
5090  IF  J=YPT  THEN  5160 
5100  REM 

5110  REM  MOVE  THE  TRAVERSER  STEP  UPWARD. 

5120  REM 
5130  I  3  =  - DIST 
5  i. 40  GOSUB  3410 
5150  NEXT  J 

5160  REM***  STORE  DATA  BEFORE  NEXT  SAMPLE*** 

5170  OPEN  "A:\RAW.DAT"  FOR  OUTPUT  AS  »2 
5180  PRINT  *2  , HEAD1S 
5190  FOR  I D= 1  TO  YPT 

5200  PRINT  12  .USING  FORMATS : ID , XPT( ID ) , YPT ( ID ) , PP ( ID . 1 ) . PP (  ID , 2 ) , PP ( I D , 3 ) . PP ( I D 
, 4 ) . PP< ID. 5 ) . YAW( ID) 

5210  NEXT  ID 

5220  CLOSE  »2  .  '  ' 

5230  REM  ***  GENERATING  INCREMENTED  DATA  FILE  NAME 
5240  IF  (I  >  10)  OR  (I  =  10)  THEN  IS  =  MI DS ( STRS ( I ) ,  2) 

5250  IF  (I  <  10)  THEN  IS  =  ( MI DS ( STRS ( 0 ) .  2)  +  MIDS ( STRS ( I ) ,  2))  , 

5260  FI2S  =  (F2S  +  IS  +  EXS ) 

5270  PRINT  1 

5280  PRINT  "  WRITING  DATA  FILE  FI2S 

5290  DF2S=RN2S^FI2S 

5300  REM  **  RENAME  DATA  FILE 

5310  SHELL  DF2S 

5320  REM 

5330  REM  MOVE  THE  TRAVERSER  TO  THE  NEXT  SAMPLE  POSITION 
5340  REM 
5350  PRINT 

5360  IF  I  =XPT  THEN  5430 

5370  INPUT  "THEN  PRESS  'ENTER'  FOR  NEXT  COLUMN  SAMPLE (  90  SEC)  MOVES 

5390  12= -DIST 

5400  I  3  =  HT 

5410  GOSUB  3410 

5420  NEXT  I 

5430  CLS 

5440  PRINT  "ALL  MOVEMENTS  COMPLETE" 

5450  PRINT 
5460  PRINT 

5470  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER  MOVEMENT  (Y  OR  N)?" 

5480  PRINT 

5490  PRINT  "IF  ' Y  *  .  THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL  CONTROL  SUBROUTINE." 

5500  PRINT  "IF  'N',  THE  PROGRAM  WILL  END." 

5510  PRINT 

5520  INPUT  "ANOTHER  MOVEMENT" ;RS 
5530  IF  RS  =  "Y"  THEN  1370 
5540  IF  RS  =  "N"  THEN  3370 
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APPENDIX  B.  CALP  PROGRAM 


1  DEF  SF.G :  CLEAR  ,&HFEOO:GOTO  4  'Begin  PCIB  Program  Shell 

2  GOTO  1000  '  User  program 

3  GOTO  900  '  Error  handling 

4  I  =  SrUFEOO  1  Copyright:  Hewlett-Packard  1984,1985 

5  PCIB . DIRt -ENVIRON? ( "PCIB" ) 

6  It=PCIB.DIRt+ "\PCI BI LC . BLD" 

7  BLOAD  I  * , I 

8  CALL  I  (PCIB.  DIR*  ,  1%,  J7.)  :PCIB.SEG=I% 

9  IF  J%=0  THEN  GOTO  13 

10  PRINT  "Unable  to  load."; 

11  PRINT  "  (Error  )" 

12  END 

13  ’ 

14  DEF  SEG=PCIB.SEG:O.S=5:C.S=10: I .  V  =  1  5 

15  I .C=20:L.P=25:LD.FILE=30 

16  GET . MEM  =  35:L.S=40:PANELS=45:DEF. ERR- 50 

17  PC  I B . ERRJ  =  STRI NG$ ( 64,32)  ;  PCIB. NAMEt  =STRI NCi (16,32) 

18  CALI  DEF. ERR( PCIB. ERR , PCIB. ERR$ , PCIB. NAME* , PCIB. GLBERP.  )  :  PCI B . 8ASERR= 2 55 

19  ON  ERROR  GOTO  3 

20  J=-l 

21  IJ=PCIB.DIR*+"\PC1B.SYN" 

22  CALL  O. S( 1$  ) 

23  IF  PCIB.ERRO0  THEN  ERROR  PCIB.  BASERR 

24  1=0 

25  CALL  I . V ( I , READ. REGI STEP, READ. SELFID, DEFINE, INITIALIZE. SYSTEM) 

2  6  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

2  7  CALL  I . V ( I , ENABLE . SYSTEM , DISABLE . SYSTEM , IN  I TI ALI ZE , POWER . ON ) 

28  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 
2Q  CALI.  I . V( I .MEASURE, OUTPUT, START. HALT) 

30  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR  ! 

31  CALL  I  I .ENABLE. INT. TRIGGER, DISABLE. INT. TRIGGER. ENABLE. OUTPUT, DISABLE. OUTPU 
T) 

32  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

33  CALL  I .V( I .CHECK. DONE. GET. STATUS, SET. FUNCTION, SET. RANGE) 

34  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

35  CALL  I ,V( I .SET. MODE, WRITE. CAL , READ. CAL. STORE. CAL) 

36  IF  PCIB.ERROO  THEN  F.PROR  PCIB. BASERR 

37  CALL  I. V( I, DELAY, SAVE. SYSTEM, J.J) 

38  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

39  1  =  1 

40  CALL  I . V ( I .SET. GATETIME. SET. SAMPLES, SET. SLOPE, SET. SOURCE) 

41  IF  PCIB.ERROO  THEN  EPROR  PCIB. BASERR 

42  CALL  I . C( I .FREQUENCY, AUTO. FREQ, PERIOD, AUTO. PER) 

43  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

44  CALL  I .C( I , INTERVAL, RATIO, TOTALIZE, R100MILLI ) 

4  5  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

46  CALL  1 .C( I ,R1,R10,R100,R1KILO) 

47  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

48  CALL  I . C ( I . R10MEGA . R100MEGA , Cl  IAN . A , CHAN . B ) 

4  9  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

50  CALL  I .C< I .POSITIVE, NEGATIVE, COMN, SEPARATE) 

51  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

52  1  =  2 

53  1  =  3 


79 


54  CALL  I  .  Vt  I  ,  ZERO. OHMS, SET. SPEED.  J  ,  j  ) 

55  IF  PCIB.ERROO  THEM  ERROR  PC  I  D  .  BASERR 

56  CALL  I . C ( I , DCVOLTS , ACVOL  > . OHMS , R200MI LLI ) 

57  IF  PC I B . ERR<  >0  THEN  ERROR  PC  I B . BASERR 

58  CALL  I . C ( I , R2 , F20 . R200 , R2KI LO ) 

59  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

60  CALL  I .C( I , R20KI LO , R200K1 LO , R2MEGA , R 2 OMEGA ) 

61  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

62  CALL  I . C ( I , AUTOM , R2 . 5 , R1 2 . 5 , J  ) 

63  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

64  1  =  4 

6  5  CALL  I . V ( I , SET . COMPLEMENT , SET . DR  I VER , OUTPUT . NO . WAIT , ENABLE . HANDSHAKE ) 

66  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

67  CALL  I . V( I , DISABLE. HANDSHAKE, SET. THRESHOLD. SET. START. BIT, SET. NUM. BITS  ) 

68  IF  PCIB.ERROO  THEN  ERROR  PC  I  B  .  BASERR 

69  CALL  I . V ; I , SET . LOGIC . SENSE, J . J , J  ) 

70  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

71  CALL  I. C(I, POSITIVE. NEGATIVE, TWOS. UNSIGNED) 

72  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

73  CALL  I .C( I , OC.TTL.RO ,R1 ) 

74  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

75  CALL  I.C(I,R2,P3,R4,R5) 

76  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

77  CALL  I . C ( I , R6 . R7 . RR , R9  ) 

78  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

79  CALL  I . C( I ,P10 , Rll , R12 . R13 > 

80  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

81  CALL  I . C( I , R14 , R15 , RI6 . I  i 

82  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

83  1  =  6 

84  CALL  I . V ( I .SET. FREQUENCY, SET. AMPLITUDE, SET. OFFSET, SET. SYMMETRY) 

8  5  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

8  6  CALL  I  .Vd.SET.  BURST.  COUNT.  J.J.J) 

87  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

88  CALL  I .C< I .SINE, SQUARE, TRIANGLE, CONTINUOUS)  , 

89  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

90  CALL  I .C( I .GATED. BURST. J.J) 

91  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

92  1  =  7 

93  CALL  I  .  V  (  I  .AUT0.SCALE. CALIBRATE, SET. SENSITIvn  Y.SET.  VERT. OFFSET) 

94  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASEPR 

9  5  CALL  I . V ( I , SET. COUPLING, SET. POLARITY . SET . SWEEPSPEED . SET . DELAY  ) 

96  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

97  CALL  I . V ( I .SET. TRIG. SOURCE, SET. TRIG. SLOPE, SET. TRIG. LEVEo SET. TRIG. MODE) 

98  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

99  CALL  I  . V ( I , GET. SINGLE. WE, GET. TWO.WF, GET. VERT. INFO , GET . TIMFBASF. . INFO) 

100  IF  PCIB.ERROO  THEN  EPROR  PCIB. BASERR 

101  CALL  I . V ( I , GET. TRIG. INFO . CALC . WFVOLT , CALC . WFTIME , CALC . WF . STATS ) 

102  IE  PCIB.ERROO  THEN  ERROR  PCIB.PASEPR 

103  CALL  I  .V<  I  ,  CALC.  RISETIMF.,  CALC.  r ALLTIME,  CALC.  PERIOD,  CALC.  FREQUENCY  ) 

104  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

105  CALL  I .V( I . CALC. PLUSWIDTII. CALC. MINUSWIDTH, CALC. OVERSHOOT, CALC. PRESHOOT) 

106  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

107  CALL  I .V( I . CALC. PK. TO. PK, SET. TIMEOUT, SCOPE. START, MEASURE. S I NGLE.WF) 

108  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

109  CALL  I .V( I .MEASURE. TWO.WF, J ,J ,J ) 
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110  IF  PCIB.ERROO  THEN  FRROR  PCIB. BASERR 

111  CALL  I  .  C  (  1  ,  R  ION  A  NO . RIOONANO ,  R1M1  CP.O ,  R1  0M 1  CRO  ) 

112  IF  PC  I B . ERR  O  0  THEM  ERROR  PCIB. BASERR 

113  CALL  1 . C( I , R1 0OMI CRO , R1MILLI , R10M1 LI.I  , R100MILLI ) 

114  IF  PCIB.ERROO  THEN  ERROR  PCI B .  BASERR 

115  CALL  1 . C ( I . R1 , RIO , R20UANO , R200UANO ) 

116  IF  PCIB.ERROO  THEM  ERROR  PCIB. BASERR 

117  CALL  I . C ( I , R2MICRO . R20MI CRO ,R200MICRO,R2MILLI ) 

118  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

119  CALL  I .C( I .R20M1LLI .R200MILLI .R2.R20) 

120  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

121  CALL  I . C( I , R50NANO, R500NANO, R5MICRO, R50MI CRO ) 

122  IF  PCIB.ERROO  THEM  ERROR  PCIB. BASERR 

123  CALL  I . C ( I , R500MICRO . R5MILLI , R50MI LLI , R^OOMILLI ) 

124  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

125  CALL  1 . C ( I , R5 ,R50 , CHAM . A . CHAN . B ) 

126  IF  PCIB. ERR< 20  THEM  ERROR  PCIB. BASERR 

127  CALL  I. C( I. EXTERNAL. POSITIVE. NEGATIVE, AC) 

128  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

129  CALL  I . C ( I , DC , TRIGGERED , AUTO . TRIG , AUTO . LEVEL ) 

130  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

131  CALL  I .C< I .XI ,X10, STANDARD. AVERAGE) 

132  IF  PCIB.ERROO  THEN  ERROR  PCIB.  BASERR 

133  1-8 

134  CALI,  I  .  V  (  I  , OPEN  .  CHANNEL. ,  CLOSE  .  CHANNEL  ,  J  ,  J  ) 

135  IF  PCIB.ERROO  THEN  ERROR  PC1B.BASF.RR 

136  CALL  C.S 

137  IF  PCIB.ERROO  THEM  ERROR  PCIB.  BASERR 

138  If'PCIB.DlRJ* "\PCI B . PLD" 

139  CALL  L.P< If  ) 

140  IF  PCIB.ERROO  THEN  ERPOR  PCIB. BASERR 

141  I  $  =  " PMM . 01" : I  -  3 : J -0 ;K  =  0 : L- 1 

142  CALL  DF.F  I  HE  (  PMM  .01, If, I, J.K.L) 

143  IF  PCIB.  ERR oq  THEN  ERROR  PCIB. BASERR 

144  I  f  =  "F  uric  .  Geri .  0  1  "  :  I  ~  6  :  J  =  0  :  K  =  1 :  L  =  1 

145  CAI.L  DEFINE!  FUNC .  GEN  .01, It, I, J.K.L) 

146  IF  PCIB.ERROO  THEN  ERROR  PC  IB.  BASERR 

147  1$^=  "Scope.  01":  1*7  :J-0:K-2:L=1 

148  CALL  DEFINE ( SCOPE ,01,It,I,J,K,L) 

149  IF  FC I B  .  ERR''  >0  THEN  ERROR  PCIB.  BASERR 

150  I  $  =  "Coi’nter  .01":  T  ■  1 :  J  =  0  :  K  =  3  :  L=  1 

151  CAI.L  DEF 1  ME  (  COUNTER  .01, If, I, J.K.L) 

152  IF  PCI D . ERR<  >0  THEN  ERROR  PCIB. BASERR 

153  IJ-’T  . In. 01" : 1-4:J=0:K=4:L*1 

154  CAL^  -EFINEI DIG.  IN. 01 , I  $ . I , J ,K,L) 

155  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

156  It -"Dip .Out .  01”: 1  =  4 : J=1 :K  =  4  :L=1 

157  CALL  DFF 1NE ( DIC . OUT .01,I$,I,J.K,L) 

158  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

159  If = "Re lay  Act .01 " :I=8:J=G:K=5:L=1 

160  CALL  DEFINE! RELAY. ACT. 01. It. I . J.K.L) 

161  IF  PCIB.ERROO  THEM  ERROR  PCIB.  BASERR 

162  If =  "Pel ay .Mux . 01 " : I =2 : J  =  0 :K=6 : L  = 1 

163  CALL  DEFINEIRELAY. MUX. 01,1$. I. J.K.L) 

164  IF  PCIB.ERROO  THEN  ERROR  PCIB. BASERR 

800  1 t  =  F.NV  I  DONf  (  "PANELS”  )  +"\PAHELS  .  EXE" 

801  CALL  L. S< If  ) 

899  GOTO  2 
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900  IF  ERR=PCIB .  BRSF.RR  THEN  GOTO  903 

901  FRINT  "BASIC  error  »";ERR;"  occurred  in  line  " ; ERL 

902  STOP 

903  TMPF.RR  =  PC.B.  ERR:  IF  TMPERR  =  0  THEil  TMPERR=  PCI  B  .  GLBERR 

904  PRINT  "PC  Instrument  error  #";TMPERR;"  detected  at  line  ";ERL 

905  PRINT  "Error:  ";PCIB.ERR* 

906  IF  LEFTS ( PCI B . NAMES , 1 ) < >CUR# ( 32 )  THEN  PRINT  "Instrument:  ";PCIB. NAMES 

907  STOP 

908  COMMON  PCI B . DIRS , PC  1 B . SEG 

909  COMMON  LD . FILE , GET . MEM , PANELS , DEF . ERR 

9 1 0  COMMON  PC  I B . BASERR , PC I B . ERR , PC I B . ERRS , PC I D . NAMES , PCI B . GLBERR 

911  COMMON  READ. REGISTER, READ. SELFID, DEFINE, INITI ALI ZE . SYSTEM , ENABLE . SYSTEM . DISA 
BLE  SYSTEM, INITI ALI ZE , POWER . ON , MEASURE , OUTPUT . START , HALT , ENABLE . I NT . TRIGGER , DISA 
BLE . I NT . TR I GGER , ENABLE . OUTPUT .DISABLE. OUTPUT , CHECK . DONE , GET . STATUS 

9 1 2  COMMON  SET . FUNCTION . SET . RANGE , SET . MOPE .WRITE. CAL , READ . CAL . STORE . CAL , DELAY . SA 
VE . SYSTEM , SET . GATETIME , SET . SAMPLES , SET . SLOPE , SET . SOURCE , ZERO . OHMS . SET . SPEED , SET . 
COMPLEMENT , SET . DRIVER . OUTPUT . NO . WA IT , ENABLE . HANDSHAKE , DISABLE . HANDSHAKE 

9 1 3  COMMON  SET  . THRESHOLD ,  SET  .  START  .BIT,  SET  . NUM .BITS,  SET  .  LOGIC  .  SENSE ,  SET  .  FREQUF.NC 
Y , SET . AMPLITUDE . SET . OFFSET , SET . SYMMETRY , SET . BURST . COUNT , AUTOSCALE , CALIBRATE , SET . 
SENS ITIVITY, SET. VERT. OFFSET , SET. COUPLING, SET. POLARITY, SET. SWEEPSPEED 

914  COMMON  SET. DELAY, SET. TRIG. SOURCE. SET. TRIG. SLOPE, SET. TRIG. LEVEL. SET. TRIG. MODF. 
, GET . S I NGLE . WF , GET . TWO . WF , GET . VERT . I NFO , GET . TI MEBASE . I NFO , GET . TR I G . I NFO , CALC . WFV 
OLT, CALC. WFT1ME, CALC. WF. STATS, CALC. RISETIME, CALC. FALLTIME, CALC. PERIOD 

915  COMMON  CALC . FREQUENCY , CALC . PLUSWIPTH , CALC . MINUSWI DTH . CALC . OVERSHOOT . CALC . PRE 
SHOOT , C ALC . PK . TO . PK , SET . TIMEOUT , SCOPE . START . MEASURE . S I NGLE . WF , MEASURE . TWO . WF , OPE 
N.CHANNEr  .CLOSE. CHANNEL 

916  COMMON  FREQUENCY, AUTO. FREQ, PERIOD. AUTO. PER, INTERVAL, RATIO, TOTAL I ZE , R100MI LLI 
, R1 . RIO, R100. R1KI LO,P.  IOMEGA, R100MEGA. CHAN. A, CHAN. B, POSITIVE. NEGATIVE, COMN.SEPARA 
TE , DC VOLTS , AC VOLTS , OHMS , R200MI LLI . R2 , R20 . R200 . R2K I LO , R20KI LO . R200KILO 

9 1 7  COMMON  R 2MEGA , R20MEGA , AUTOM , R2 . 5 , P 1 2 . 5 , POS I TI VE , NEGAT I VE , TWOS . UNS I GNED , OC , TT 
L,RO,Rl,R2,R3,R4,R5,R6,R7,R8,R9,R10.R13 , R1 2 . R 1 3 , R 1 4 , R1 5 , R1 6 , SINE , SQUARE . TRI ANGLE 
, CONTINUOUS , GATED , BURST , RIOMANO, R ICON A NO . R1MICRO , PIOMICRO , RIOOMICRO 

918  COMMON  R1MILLI .R10MILLI ,R100M1LLI , P 1 , RIO , R20NANO , R200NANQ ,  P.2MICRO , R20MI CRO , R 
200MICRO.R2MILLI .R20MILLI , R200MILLI ,  P.2 , R20 . R50NANO, RSOONANO , R5MI CFO , P 5 PM I CRO , R50 
OM I  CP.O  ,  R5MI  LLI  .R50MILLI  .R500MILLI  ,  R5  ,  P.50  ,  CHAN  .  A  ,  CIIAN  .  B  ,  EXTERNAL  .  POSIT  I  VE 

9 1 9  COMMON  NEGATIVE . AC , DC , TR I GGERED , AUTO .TRIG, AUTO . LEVEL , X 1 , X 1 0 , STANDARD , AVERAGE 

920  COMMON  DMM .01. FUNC . GEN . 0 1 , SCOPE .01. COUNTER . 01, DIG. tn. 01, DIG. OUT .01. RELAY . ACT 
. 01 .RELAY  MUX .01 

999  ’End  PCIB  Program  Shell 

1000  REM 

1010  REM  This  step  initial zes  the  HP  system 
1020  CLS 

1030  OPTION  BASE  1 

1040  DIM  P( 10 ) , PA( 50 , 6 ) , PP( 50, 6 ) , XPT( 40 ) , CAL( 40 ) 

1050  CALL  INITIALIZE. SYSTEM) PGMSHEL . HPC ) 

3060  REM 

3070  REM  All  PC  devices  now  have  an  initial  state 
1080  REM  Set  function  on  the  DMM  and  Relay  MUX 
1090  REM 

1100  CALL  SET. FUNCTION! DMM. 01 , DC VOLTS) 

1110  CALL  SET. RANGE (DMM. 01, AUTOM) 

1120  CALL  DISABLE. INT . TRIGGER ( DMM . 01 ) 

1130  CALL  ENABLE. OUTPUT (RELAY. MUX. 01 ) 

1140  FORMAT*  =  "M  M.MM  *  It  .  »  »  «  I  M.MM  M.MM  M  .  I  Ml  (III. MM" 

1200  FOR  1  =  1  TO  10 
1210  CAL ( I  )  =0 . 0 
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1220  NEXT  I 
1510  REM 

1520  REM  READ  THE  VOLTAGE  OF  48TH  CHANNEL  AND  DISPLAY  THE  DATA 
1530  REM 

15<.0  PRINT  "  CHOOSE  6  POINTS" 

1550  PRINT 

1550  PRINT  "THE  CALIBRATION  WILL  BE  STORES  IN  'CAL. DAT'" 

1560  REM 

1570  REM  Begin  sampling  loop 
1580  REM 

1600  FOR  J  =  1  TO  1 
1610  PRINT 

1630  FOR  JJ=1  TO  6 

1631  INPUT  "INPUT  THE  CALIBRATION  PRESSURE" ; CAL ( JJ ) 

1632  INPUT  "PRESS  'ENTER'  TO  START  MEASUREMENT” ; MOVES 
1640  FOR  11=1  TO  10 

1650  ROUT= 1 

1660  CALL  OUTPUT (RELAY. MUX. 01, ROUT) 

1670  CALL  MEASURE ( DMM . 01, VOLTS) 

1680  I’A (II.  JJ  )  VOLTS 

1690  NEXT  II 

1700  IF  J  J  =  6  THEN  1740 

1730  NEXT  JJ 

171,0  REM 

1750  REM  DISPLAY  THE  SAMPLE  DATA 
1760  REM 

1780  FOR  IS=  1  TO  10 

1790  PRINT  USING  FORMATS : IS , PA( 1 S , 1 ) , PA( IS , 2 ) , PA( IS , 3 ) . PA ( IS , 4 ) . PA ( IS , 5 ) , PA ( IS , 6 
) 

1800  NEXT  IS 
1810  REM 

1820  REM  AVERAGE  THE  DATA 

1830  REM  1 

1840  FOR  JA  =  1  TO  6 

1860  TOTAL  =  0  < 

1860  FOR  IA  =  1  TO  10 

1870  TOTAL  =  TOTAL  +  PA(IA.JA) 

18  80  NEXT  I A 

1890  AVERAGE  =  TOTAL  /10 
1900  F ( JA ) = AVERAGE 
1920  NEXT  JA 
1930  PRINT 

1940  PRINT  "THE  AVERAGE  ARE:  " 

2000  FOR  JP=1  TO  6 
2010  PP( J , JD)=P< JD) 

2020  NEXT  JD 

2055  PRINT  USING  FORMATS ; J , PP ( J , 1 ) , FP ( J , 2 ) , FP ( J , 3  )  , PP ( J , 4 )  , PP ( J , 5 )  , DP ( J , 6  ) 

2070  PRINT 

2080  INPUT  "DO  YOU  WANT  P.E  MEASURE  AGAIN  ?  (Y  /  N ) "  ;  CS 
2090  IF  CS  =  " Y ”  THEN  1580 

2101  REM***  STORE  DATA  BEFORE  NEXT  SAMPLE*** 

2102  OPEN  "A:\CAL.DAT"  FOR  OUTPUT  AS  12 

2106  FOR  1D= 1  TO  6 

2107  PRINT  *2, USING  FORMATS ; I D . PP ( J , I D ) , CAL ( 1 D ) 

2108  NEXT  ID 

2109  CLOSE  #2 
2210  NEXT  J 
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APPENDIX  C.  CONVERT  PROGRAM 


ft**************************************************************** 

*  THIS  PROGRAM  CONVERTS  THE  VOLTAGE  OF  TRANSDUCER  INTO  PHYSICAL  * 

*  PRESSURE, VELOCITY, YAW  ANGLE  AND  PITCH  ANGLE.  THOSE  DATA  ARE  * 

*  USED  FOR  PLOT  PROGRAM  LATER.  * 

MUtMItttttttmittltttMOtlMtltttlttMMtllHMMMItnXttt 

CHARACTER* 12  FNAME 
CHARACTER* 12  NAME 
CHARACTER* 2  A(S0> 

CHARACTER*80  ST 
REAL  K, INTR 

INTEGER  COLS , RWS , DTPTS 


DATA  A/’Ol ’ . 

'02  '  . 

'03  ' 

,  '04' 

.  '05'  ,  ’ 

1  06'  , 

’  07  ' 

,  '08'  , 

'09' 

* 

'  10'  , 

'll’, 

*12' 

.  ’  13' 

.  '14'  .  ' 

1  16’  . 

’  16’ 

,  '  17  ’  , 

'  18  ' 

* 

'  19  '  . 

'  20'  , 

'21' 

.  ’  22  ’ 

.  ’  23  '  . 

'  24  '  , 

'  25  ' 

,  '  26  '  , 

'  2  7  ' 

* 

'  28  '  . 

'  29'  , 

'  30* 

.  '  31 ' 

,  '32'  , 

'  33  ’  . 

'  34  ' 

.  '  35’  . 

'  36  ' 

* 

'  37  '  , 

•38'. 

'  39  ' 

,  '  40' 

,  '41'  , 

’  42  '  , 

'43' 

,  '44  ’  , 

’  45  ' 

* 

'46'  , 

’  47  '  . 

'  48  ’ 

,  ’  49' 

,  ’  50’  / 

WRITE 

(*, ’ ( A\ ) ’ ) 

’  *  OF  COLS 

(AWAY  FROM  MSL )  =  ’ 

READ 

(  ♦  ,  ’  (  15  )  ’  ) 

COLS 

WRITE 

(*,  ’ (A\)’  ) 

’  I  OF  DATA 

PTS  IN  A  COL  (UP/DOWN)  = 

READ 

<  »  ,  '  <  15  > '  ) 

RWS 

WRITE 

(  ♦  ,  ’  (A\ ) '  ) 

’  DATA  FILE 

NAME?  (IER001A2XX.DAT) 

READ 

(  *  ,  '  (A12  )•’  ) 

NAME 

WRITE 

(*  ,  '  (A\) '  ) 

'  PI  (F4.2) 

_  » 

READ 

(  * ,  ’  (F4 . 2 ) ' 

)  PI 

WRITE 

(*,  '  (A\) '  ) 

’  PF  (F4.2) 

=  ' 

READ 

(  * , ' ( F4 . 2 ) ' 

)  PF 

WRITE 

(*. ' <A\) ' ) 

’  TI  (F3.1) 

_  ( 

READ 

(  *,  ' (F3. 1 ) ’ 

)  TI 

WRITE 

(*,  ' ( A\ ) '  ) 

’  TF  (F3.1) 

_  |  ' 

READ 

(  *  ,  '  (  F3. 1 ) ' 

!  TF 

WRITE 

(*,  ' <A\) ’  ) 

II 

v£> 

u, 

* 

f 

READ 

(  * ,  '  ( F6. 4 ) ' 

)  K 

WRITE 

(* ,  '  <A\)’  ) 

’  SLOPE  FOR 

DELTAP  (F9.6)  =  ’ 

READ 

(♦,  '  (F9.6) ' 

)  SLOrE 

WRITE 

( * , ' <A\) ’ ) 

’  INTERCEPT 

FOR  DELTAP  (F9.6)  =  ' 

READ 

(  * ,  ' (F9. 6)  ' 

)  INTR 

WRITE 

( * ,  ’  ( A\  )  '  ) 

’  OM1  FACTOR 

!  (F4.2)  =  • 

READ 

(  * ,  ' ( F4 . 2  )  ' 

)  QM1FAC 

WRITE 

(  *  ,  '  <A\  i '  ) 

’  X  OFFSET  = 

1 

READ 

(»  ,  * (F5. 2 ) ' 

)  XOFF 

WRITE  ( * ,  '  ( A\  )  1  )  '  Y  OF  FSET  =  1 

READ  (♦, ’ (FS.2) ’ )  YOFF 

*  CONVERT  THE  PRESSURE  UNIT  FROM  inlle;  TO  psf 
PATM=(PI+PF>*35. 3631 
R= 17 1 6 . 5 
E-0.0123 

T= (TI +  TF 1/2 . *460 
R0=PATM/(R*T) 

DTPTS  =  RWS*  COLS 
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*  OPEN  A  NEW  FILE  TO  STORE  THE  REDUCED  PftTft 

OPEN ( 2 , FI  LE  = 'RESULT. DAT' , STATUS = ' OLD ' ) 

WRITE (2,222)  DTPTS 
222  FORMAT  (15) 

*  OPEN  ft  SEQUENTIAL  OF  DATA  FILE 

DO  20  1=1, COLS 
NAME ( 7:8) =A (I) 

FNAME=NAME 

OPEN (  1 .  FI  I,E  =  FNAME  ) 

READ( 1 , 100 , END=20 )ST 
100  FORMATt  fl65 ) 

15  READ ( 1 , 1000, END=30 )NO,X,Y, VI , V2 , V3 ,V4 , V5 .BETA 

1000  FORMAT (12, F7 . 2 , F6 . 2 , 5F9 . 3 , F8 . 2  ) 

*  CONVERT  THE  VOLTAGE  TO  PRESSURE  IN  LBF/FT**2 

P 1  =  DELTA P  (  VI  ,  SLOPE  ,  INTR  )  +  2 . 04  7  5  +  P ATM 
P2  =  DELTAP  (  V2  ,  SLOPE  ,  1 NTR  )  *  2 . 04  7  5  +  PATM 
P3=DELTAP(V3, SLOPE, INTR ) *  2 . 0  4  7 5+FATM 
P4 = DELTAP (  V4 ,  SLOPE  ,  INTR  )  *  2  .  047  5  +  PATM 
P5= DELTAP ( V5 , SLOPE, INTR >* 2 . 047 5+ PATM 

*  CALCULATE  THE  PITCH  ANGLE  IN  DEGREES 

F  =  I  P4-P5  )/(  P1-P2  ) 

ALPHA=FPITCH ( P ) 

*  CALCULATE  THE  VELOCITY  IN  FT/SEC 

YSLOF =FYSLOP ( ALPHA ) 

VELM-SORT( ( 2* YSLOP* ( F1-P2 ) )/(RO+K) ) 

VEL.=VELM*  (  l  +  E) 

*  CALCULATE  THE  LOCAL  DYNAMIC  PRESSURE 

0M1=0M1FAC*2 . 0475/K 
QM  =  RO+  VEL*  *  2/2  . 

Q 1  =  QM 1  * ( 1  +  2  *E  ) 

Q  =  QM  + ( 1+2*E> 

*  CALCULATE  THE  YAW  ANGLE  IN  DEGREES 

YAW-FYAW; BETA + 5 . 0  ) 

*  CALCULATE  THE  VELOCITY  COMPONENTS 

BETAR=YAW*. 017453 

ALPHAR*  (ALPJiA-17.  9  4  2  )  *  .  017453  1 

VELY=VEL*SIN( ALPHAR ) 

VELX=VEL+COS  (  ALPHAR  )  *SIN  (  BF.TAR  ) 

*  CALCULATE  THE  TOTAL  PRESSURE  IN  LBF/IN**2 

FTC- FFT ( ALPHA  ) 

PT1 =  F1 -Q* PTC 
PT=PT1/144 . 

CPT= ( TT1 - PATM-Q1 )/Ql 

*  CALCULATE  THE  STATIC  PRESSURE  IN  LBF/IN**2 

PS  1  =  PT1 -O 
PS  =  PS 1/14  4 

cfs=(psi-patm>/qi 

*  WRITE  ZERO  VALUES  IF  VELOCITY  IS  TOO  HIGH  (BAD) 

I F ( VEL . GT . 200.0)  THEN 
2  =  0.0 

WRITE  ( 2 , 2000 ) -X  +  XOFF , Y+YOFF , Z , Z ,Z,Z,Z,Z,Z,Z,Z 
ELSE 

WRITE ( 2 , 2000 ) -X+XOFF , Y+YOFF , VEL , VELX , VELY , YAW , 
C  ALPHA-17. 942, FT.CPT, PS, CPS 
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2000 


FORMAT ( 11F10 . 3 ) 

END  IF 
GO  TO  15 
30  CLOSE ( 1 ) 

20  CONTINUE 

CLOSE ( 2 ) 

STOP 

END 

************************************************************* 

*  THIS  FUNCTION  CONVERTS  THE  VOLTAGE  TO  PHYSICAL  PRESSURE 

FUNCTION  DELTAP (X, SLOPE, INTR! 

REAL  INTR 

P£LTAP=X  *SI,OPE+ 1 NTR 
END 

***** *  ************************** 

*  THIS  FUNCTION  CALCULATES  THE  PITCH  ANGLE 

FUNCTION  FPITCH(X) 

FPITCH=  3 . 759+53. 7568 +X-1 . 3C85*X* *2-1. 6583*X» *3 
*  -0. 8061*X**<*  +  16. S115*X**5 

END 

************************************************************ 

*  THIS  FUNCTION  CALCULATES  THE  VELOCITY  PRESSURE  COEFFICENT 

FUNCTION  FYSLOP ( X ) 

IF(/..LT. -iOJTHEN 

FYSLOP =0 .981-0. 0102+X-3.000E-4+X **2-2. 500E-6+X**3 
ELSE  I F ( (X.GE. -10) .AND. (X.LE. 10) (THEN 
FYSLOP  =  0 .98-0.006*X+2.000E-4*X**2 
ELSE 

FYSLOP  =  0. 9801-0. 0035‘X-l.  l(i3E-4*X* *2  +  5 . 833E-6*X*» 3 
END  IF 
END 

************************************************************* 

*  THIS  FUNCTION  CALCULATES  THE  YAW  ANGLE 

FUNCTION  FYAW(X)  ' 

IF( (X.GE.O ) .AND. (X.LE. 180) )  THEN 
FYAW= -  X 
ELSE 

FYAW=  360  -  X 
END  IF 
END 

************************************************************* 

*  THIS  FUNCTION  CALCULATES  THE  TOTAL  PRESSURE  COEFFICIENT 

FUNCTION  FPT(X) 

IFIX.LE.-30)  THEN 
FPT= -0 . 0 1 


ELSE  IF( (X.GT. -30) .AND. 
FPT=  0 .02  +  1.00E-3  +  X 

(X.LT.-20) 

)  THEN 

ELSE  IF( (X.GE. -20) .AND. 

FPT -0 

ELSE 

FPT =0. 03-1. 00E-3*X 

(X.LE. 30) ) 

THEN 

END  IF 
END 
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APPENDIX  D.  RESULT  00. DAT 


8 4 

1 

.000 

.000 

122.209 

-7.000 

18.214 

15.002 

.017 

1 

(..875 

2 

-.060 

.000 

.500 

122.155 

-7.000 

18.388 

15.002 

.015 

1 

4.875 

3 

-  .061 
.000 

1.000 

122.267 

-7.000 

18.254 

15.002 

.017 

1 

4.875 

4 

-  .062 
.000 

1.500 

122 . 199 

-7.000 

18.297 

15.002 

.016 

1 

4.875 

5 

-  .061 
.000 

2.000 

122.310 

-7.000 

18.433 

15.002 

.018 

1 

4.875 

6 

-  .061 
.000 

2.500 

122.155 

-7.000 

18.388 

15.002 

.016 

1 

4.875 

7 

-.061 

.000 

3.000 

122 . 293 

-7.000 

13.302 

15.002 

.018 

1 

4.875 

8 

-  .  061 
.000 

3  .  500 

122 .359 

-7.000 

18.287 

15.002 

.018 

1 

4.875 

9 

-  .062 
.000 

4 . 000 

122.287 

-7.000 

18 . 357 

15.002 

.017 

1 

4.875 

10 

-  .062 
.000 

4 . 500 

125 . 471 

-5.000 

17.651 

15.002 

.017 

1 

4.868 

11 

-  .  118 
.  000 

5.000 

125.733 

-5.000 

17.773 

15.002 

018 

1 

4.868 

12 

-.122 

.000 

5.500 

125.881 

-5.000 

17.614 

15.002 

.  019 

1 

4.867 

1 

-  .  124 
.500 

.000 

125.918 

-5.000 

17.7B4 

15.002 

.017 

1 

4.867 

2 

-.127 
.  500 

.  500 

125.733 

-5.000 

17.773 

15.002 

.018 

1 

4.868 

3 

-  .122 
.500 

1.000 

126.131 

-5.000 

17.814 

15.002 

.015 

1 

4.866 

4 

-.132 

.500 

1 . 500 

125.773 

-5.000 

17.714 

15.002 

.016 

1 

4.867 

5 

-.125 
.  500 

2.000 

126.015 

-5.000 

17.764 

15.002 

.015 

1 

4.867 

6 

-  .131 
.  500 

2.500 

125 .812 

-5.000 

17.655 

15.002 

.013 

1 

4.867 

7 

-.128 

.500 

3.000 

125.624 

-5.000 

17  669 

15.001 

.012 

1 

4.867 

8 

-  .  126 
.500 

3 . 500 

125.731 

-5.000 

17.570 

15.002 

.015 

1 

4.867 

9 

-  .  126 
.500 

4 . 000 

125.777 

-5.000 

17 . 688 

15.002 

.014 

1 

4.867 

10 

-  .  127 
.500 

4 . 500 

125 . 944 

-5.000 

1  a.  1 

15.002 

.013 

1 

4.867 

11 

-.131 

.500 

5.000 

126.524 

-5.000 

17  . 

15.002 

.015 

1 

4.865 

12 

-.140 
.  500 

5.500 

125.918 

-5.000 

17.784 

15.002 

.015 

1 

4.867 

-  .  129 
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1 

4.865 

1.000 
-  .  142 

.000 

126.500 

2 

4 . 866 

1.000 
-  .  136 

.  500 

126 . 304 

3 

4.866 

1.000 

-.134 

1.000 

126 . 189 

4 

4.866 

1.000 
-  .  134 

1.500 

125.976 

5 

4.867 

1.000 
-  .  126 

2.000 

125 . 656 

6 

4.067 

1.000 

-.127 

2.500 

125.659 

7 

4.867 

1.000 

-.126 

3.000 

125 . 571 

8 

4.867 

1.000 

-.131 

3.500 

125.821 

9 

4. 867 

1.000 

-.130 

4.000 

125.881 

10 

4.867 

1.000 
-  .  129 

4.500 

125.617 

11 

4.867 

1 . 000 
-  .  125 

5.000 

125.492 

12 

4.866 

1.000 
-  .  137 

5 . 500 

126.317 

1 

4.866 

1.500 
-  .  138 

.  000 

126.253 

2 

4 . 867 

1.500 

-.131 

.500 

125 . 896 

3 

4.866 

1.500 

-.133 

1.000 

125 .926 

4 

4.866 

1.500 
-  .  132 

1.500 

125 .860 

5 

4.867 

1.500 
-  .  126 

2.000 

125.405 

6 

4.866 

1.500 
-  .  135 

2.500 

12S.941 

7 

4.866 

1.500 
-  .  139 

3.000 

126.016 

8 

4.866 

1.500 
-  .  136 

3 . 500 

126.063 

9 

4. 866 

1.500 
-  .  135 

4.000 

125 .824 

10 

4.866 

1.500 
-  .  133 

4 . 500 

125 . 881 

11 

4.866 

1.500 

-.134 

5.000 

125 .853 

12 

4 . 866 

1.500 

-.133 

5.500 

125.784 

1 

4.866 

2.000 
-  .  133 

.000 

125.709 

2 

4.867 

2.000 
-  .  125 

.500 

125.242 

3 

4.866 

2.000 

-.138 

1 . 000 

125  .  934 

4 

4 . 866 

2.000 

-.137 

1.500 

125  .  925 

5 

4 . 866 

2.000 

-.136 

2.000 

125 . 833 

6 

4 .866 

2.000 

-.136 

2 . 500 

125.844 

-5.000 

17 . 408 

15.001 

.012 

1 

-5.000 

17.474 

15.002 

.015 

1 

-5.000 

17 . 625 

15.002 

.  015 

1 

-5.000 

17.822 

15.001 

.  Oil 

1 

-5.000 

17.662 

15.001 

.012 

1 

-5.000 

17.636 

15.001 

.012 

1 

-5.000 

17 . 604 

15.001 

.011 

1 

-5.000 

17 . 602 

15.001 

.011 

1 

-5.000 

17.614 

15.001 

.013 

1 

-5.000 

17 .722 

15.001 

.009 

1 

-5.000 

17.723 

15.001 

.010 

1 

-5.000 

17 . 396 

15.002 

.014 

1 

-5.000 

17 . 611 

15.001 

.011 

1 

-5 . 000 

17.509 

15.001 

.012 

1 

-5 . 000 

17.732 

15.001 

.011 

1 

-5.000 

17.543 

15.001 

.011 

1 

-5.000 

17.895 

15.001 

.008 

1 

-5.000 

17.627 

15.001 

» 

.009 

1 

-5.000 

17.536 

15.001 

.007 

1 

-5.000 

17.424 

15.001 

.010 

1 

-5.000 

17 . 576 

15.001 

.007 

1 

-5.000 

17 .614 

15.001 

.010 

1 

-5.000 

17 . 594 

15 . 001 

.008 

1 

-5.000 

17 . 635 

15.001 

.008 

1 

-5.000 

17.727 

15.001 

.007 

1 

-5.000 

17 .725 

15.001 

.007 

1 

-5 . 000 

17.679 

15.001 

.006 

1 

-5.000 

17.529 

15.001 

.006 

1 

-5.000 

17 . 726 

15.001 

.006 

1 

-5.000 

17.876 

,15.001 

.006 

1 
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7 

4.866 

2.000 

-.135 

3.000 

125.890 

-5.000 

17 .765 

15.001 

.008 

1 

8 

4 . 865 

2.000 
-  .  141 

3  .  500 

126 . 253 

-5.000 

17.611 

15.001 

.009 

1 

9 

4.866 

2.000 
-  .  132 

4.000 

125.720 

-5.000 

17 . 649 

15.001 

.  008 

1 

10 

4.867 

2.000 
-  .  127 

4.500 

125.289 

-5.000 

17. 613 

15.001 

.006 

1 

11 

4.866 

2.000 

-.136 

5.000 

125.982 

-5.000 

17.770 

15.001 

.009 

1 

12 

4.867 

2.000 

-.130 

5.500 

125.549 

-5.000 

17.762 

15.001 

.007 

1 

1 

4.866 

2.500 

-.139 

.000 

125.971 

-5.000 

17 . 848 

15.001 

.005 

1 

2 

4.866 

2.500 

-.137 

.  500 

125.873 

-5.000 

17 . 667 

15.001 

.006 

1 

3 

4.866 

'2.500 
-  .  136 

1.000 

125.760 

-5.000 

17 . 589 

15.001 

.005 

1 

4 

4.866 

2.500 

-.137 

1.500 

125.886 

-5.000 

17.791 

15.001 

.006 

1 

5 

4.867 

2.500 
-  .  127 

2.000 

125 . 235 

-4.000 

17.778 

15.000 

.005 

1 

6 

4.866 

2.500 
-  .  135 

2.500 

125 . 696 

-4.000 

17.603 

15.000 

.005 

1 

7 

4.867 

2.500 
-  .  124 

3.000 

125.117 

-4.000 

17.726 

15.000 

.005 

1 

8 

4.867 

2.500 
-  .  130 

3.500 

125.394 

-4.000 

17.539 

15.000 

.004 

1 

9 

4.867 

2.500 
-  .  124 

4 . 000 

125.169 

-4.000 

17.793 

15.001 

.006 

1 

10 

4.866 

2.500 

-.135 

4.500 

125.871 

-4.000 

17.464 

15.001 

.007 

1 

11 

4.868 

2.500 

-.119 

5.000 

124.912 

-4.000 

17.848 

15.001 

.006 

1 

12 

4.868 

2.500 
-  .  118 

5 . 500 

124.873 

-4.000 

17.676 

15.001 

.007 

1 

1 

4.867 

3.000 
-  .  125 

.000 

125.254 

-4.000 

17.646 

15.001 

.006 

1 

2 

4.867 

3.000 
-  .  125 

.500 

125.013 

-4.000 

17.800 

15.000 

.002 

1 

3 

4.867 

3.000 
-  .  129 

1.000 

125.231 

-4.000 

17 . 804 

15.000 

.002 

1 

4 

4.866 

3.000 
-  .  137 

1.500 

125.775 

-4.000 

17.484 

15.000 

.  004 

1 

5 

4.867 

3.000 
-  .  124 

2.000 

125.213 

-4.000 

17.706 

15.001 

.007 

1 

6 

4.867 

3.000 
-.  128 

2.500 

125.353 

-4.000 

17 . 599 

15.001 

.005 

1 

7 

4.867 

3.000 
-  .  130 

3.000 

125.550 

-4.000 

17 . 532 

15.001 

.007 

1 

8 

4.868 

3 . 000 
-  .  121 

3  .  500 

124.923 

-4.000 

17.768 

15.000 

.005 

1 

9 

4.866 

3.000 

-.135 

4 . 000 

125.558 

-4 . 000 

17 . 479 

15.000 

.002 

1 

10 

4.867 

3.000 

-.123 

4.500 

125.077 

-4.000 

17.787 

15.001 

.006 

1 

11 

4.868 

3.000 
-  .  119 

5.000 

124.878 

-4.000 

17.649 

15.001 

.006 

1 

12 

4.867 

3.000 
-  .  130 

5.500 

125.433 

-4.000 

17 .480 

15.001 

.005 

1 
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APPENDIX  E.  RESULT  OADAT 


299 


1 

(*.516 

.000 
-  .  604 

•  COO 

150.248 

-16.000 

15.148 

14.700 

-  .041 

1 

2 

4  .  502 

.000 
-  .122 

.  250 

155 .911 

-15.000 

14.475 

14 . 701 

-.039 

1 

3 

it  .  491 

.000 
-  .817 

.500 

160.315 

-15.000 

13 . 866 

14.701 

-  .037 

1 

4 

4. 480 

.  000 
-.911 

.750 

164 . 509 

-14 . 000 

13.393 

14 . 701 

-  .037 

1 

5 

4 . 465 

.  000 
-1.033 

1 . 000 

169.716 

-14.000 

12 . 697 

14.701 

-  .039 

1 

6 

4 . 448 

.  000 
-1. 177 

1 . 250 

175 . 393 

-13 . 000 

10.631 

14 . 700 

-  .047 

1 

7 

4 . 420 

.000 

-1.414 

1 .  500 

180 . 980 

-8 . 000 

2.980 

14.688 

-.147 

1 

8 

4 . 357 

.  000 
-1.955 

1.750 

185.168 

9 . 000 

-5.242 

14 . 637 

-  .  581 

1 

9 

4.293 

.000 

-2.493 

2 . 000 

179.544 

39 . 000 

-23.961 

14 . 556 

-  1 . 262 

1 

10 

4.318 

.000 

-2.281 

2.250 

181 . 551 

73.000 

-13 .886 

14.587 

-  .  999 

1 

11 

4.349 

.  000 
-2.016 

2  .  500 

199 .011 

84 .000 

-1 .061 

14.673 

-  .  274 

1 

12 

4.347 

.000 

-2.035 

2.750 

192. 122 

81 . 000 

-  .  585 

14 . 649 

-  .  480 

1 

13 

4.307 

.  000 
-2.379 

3 . 000 

190 . 920 

61.000 

-17 .689 

14 . 604 

-  .  856 

1 

14 

4.250 

.000 

-2.857 

3 . 250 

186  .  191 

24.000 

-3.723 

14 .533 

-  1 .456 

1 

15 

4.325 

.000 

-2.227 

3 . 500 

190 . 574 

4 . 000 

3  .  291 

14.621 

-  .  713 

i 

1 

16 

4  .  388 

.000 

-1.687 

3.750 

180 .810 

-8.000 

4.141 

1  4  .,679 

-  .  219 

1 

17 

4.433 

.000 
-1. 310 

4.000 

1 78.573 

-11.000 

3.863 

14.693 

-  .  102 

1 

18 

4.457 

.  000 
-1.101 

4.250 

171.008 

-13.000 

4 . 608 

14 . 696 

-  .077 

1 

19 

4 . 480 

.000 
-  .  °11 

4 . 500 

162.822 

-13.000 

6.181 

14 . 696 

-  .075 

1 

20 

4.499 

.  000 
-  .752 

4.750 

155.653 

-13.000 

8 .050 

14.697 

-  .074 

1 

21 

4.510 

.000 
-  .  656 

5.000 

151.143 

-12.000 

9.417 

14 . 696 

-  .  075 

1 

22 

4.513 

.000 
-  .  631 

5.250 

150.027 

-12.000 

10.473 

14 . 697 

-  .073 

1 

23 

4.518 

■  OOt) 

-  .  585 

5 . 500 

147.773 

-  12 . 000 

11.461 

14 . 697 

-  .073 

1 

1 

4 .515 

.250 
-  .  609 

.000 

147 .058 

-12.000 

16.164 

14 . 692 

-  .  112 

1 

2 

4.506 

.  250 
-  .688 

.  250 

152 . 564 

-16.000 

15.207 

14 . 696 

-  .077 

1 

3 

4.494 

.  250 
-  .  789 

.  500 

157 . 323 

-16.000 

14 . 221 

14 . 696 

-  .076 

1 

4 

4.481 

.  250 
-.901 

.750 

159.896 

-17.000 

14.424 

14.690 

-  .  131 

1 

5 

4.463 

.250 

-1.056 

1.000 

169.024 

-16.000 

12 . 905 

.  14 . 696 

-  .078 

1 
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6 

4.442 

7 

4  .  407 

8 

4 . 359 
9 

4 . 366 
10 
4.325 
11 
4 . 362 
12 
4.361 

13 
4.281 

14 
4 . 259 

15 
4.318 

16 
4 .385 

17 

4 . 436 

18 
4.464 

19 
4.488 

20 

4 . 505 
21 

4 . 512 
22 
4 . 514 
23 
4.516 
1 

4.512 

2 

4.506 

3 

4 . 497 

4 

4.483 

5 

4.466 

6 

4.445 

7 

4. 405 

8 

4 . 368 
9 

4 . 328 
10 
4.340 
11 
4.383 
12 
4.387 


.250 

1.250 

176.146 

-  16 . 000 

11 . 158 

14 . 695 

-  .  087 

-1.235 
.  250 

1 . 500 

183110 

-  10.000 

5 . 264 

14.681 

-  .  206 

-1.527 

.250 

1.750 

184.142 

2 . 000 

1.110 

14.636 

-  .  589 

-1. 936 
.250 

2.000 

164 . 120 

25.000 

12 . 257 

14 . 586 

-  1 . 009 

-1.674 

.250 

2.250 

187.888 

66.000 

-2.918 

14.613 

-  .779 

-2.223 
.  250 

2 . 500 

189.832 

84 . 000 

1 . 471 

14 . 656 

-  .  418 

-1.913 
.  250 

2.750 

191 .492 

81 . 000 

-3.333 

1 4 . 661 

-  .  376 

-1.915 

.250 

3 . 000 

192 . 988 

63 . 000 

-10.674 

14 . 586 

-1.014 

-2.592 
.  250 

3.250 

193 .453 

22 . 000 

-23.030 

14.565  . 

-1.189 

-2.780 

.250 

3 . 500 

203 . 195 

1.000 

-5.046 

14.655 

-  .  626 

-2 . 285 
.250 

3.750 

192.065 

-9.000 

.  336 

14.686 

-  .  159 

-1.713 

.250 

4.000 

178 .709 

-11. 000 

3 . 538 

14.896 

-  .075 

-1.286 
.  250 

4.250 

169. 131 

-  11 . 000 

4 . 692 

14.698 

-  .064 

-1 . 045 
.  250 

4 . 500 

159.931 

-12.000 

6.521 

14 . 697 

-  .071 

-  .842 
.250 

4.750 

153.467 

-  12 . 000 

7.862 

14 . 698 

-  .065 

-  .  696 
.  250 

5 . 000 

151-014 

-12.000 

9.189 

14 . 698 

-  .063 

-  .  642 
.250 

5 . 250 

149.874 

-12 . 000 

10.150 

14 . 698 

-  .063 

i 

-  .619 
.  250 

5.500 

149. 154 

-12 . 000 

11.166 

14 . 698 

-  .061 

-.602 
.  500 

.  000 

150.077 

-  14 . 000 

16.013 

14 . 696 

-  .070 

-  .  637 
.  500 

.250 

152.885 

-13 . 000 

15.585 

14 . 697 

-  .067 

-  .  685 
.  500 

.  500 

156 . 663 

- 14 . 000 

15.429 

14.698 

-  .065 

-  .  764 
.  500 

.750 

161.989 

-14.000 

14 . 704 

14 . 698 

-  .063 

-  .880 
.  500 

1.000 

168 .212 

-15.000 

14 .154 

14.697 

-.069 

-1.028 
.  500 

1.250 
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-  .  101 

1 

2 
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2 . 500 
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1 . 000 
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3 . 000 

123 . 766 
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-  .  280 

17 

2. 500 
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3.000 
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14 . 694 

-  .  090 

-2.000 

11.045 

14.694 

-  .097 

-2.000 

11 . 369 

14.694 

-  .  097 

-4.000 
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14 . 694 
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14 . 694 

-  .  067 

1 

23 

4.614 

.  000 
.024 

5 . 500 

99 . 244 

-12.000 
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3.250 

84 . 599 

15 

4.633 

2.250 

.236 

3  .  500 

86.013 

16 

4 .630 

2.250 
.  210 

3.750 

•  87 . 663 

17 

4 . 629 

2.250 
.  191 

4 . 000 

88.991 

18 

4.627 

2.250 
.  169 

4 . 250 

90. 110 

19 

4.625 

2.250 

.145 

4  .  500 

91 . 566 

20 

4 . 624 

2  .  250 
.  141 

4.750 

91.610 

21 

4 . 625 

2.250 

.144 

5.000 

91.555 

22 

4.623 

2.250 
.  130 

5 .250 

92 . 289 

23 

4.624 

2 .250 
.  134 

5 . 500 

91 . 893 

1 

4.625 

2 . 500 
.152 

.  000 

91 . 560 

2 

4.625 

2  .  500 
.  148 

.  250 

92 . 009 

3 

4.625 

2 . 500 
.  150 

.  500 

91.723 

4 

4.626 

2 . 500 
.159 

.750 

91.319 

5 

4 .627 

2.500 

.176 

1 . 000 

90.127 

6 

4.628 

2.500 

.186 

1.250 

89.567 

7 

4 . 630 

2.500 
.  203 

1  .  500 

88.310 

8 

4.632 

2 . 500 
.  231 

1 .750 

86 . 670 

-7.000  19.373  14.  695  -.059  1 

-6.000  20.099  14.695  -.059  1 

-5.000  20.561  14.694  -.066  1 

-5.000  21.612  14.694  -.064  1 

-4.000  22.387  14.694  -.065  1 

-2.000  23.316  14.694  -.066  1 

1.000  23.870  14.694  -.066  1 

3.000  24.368  14.694  -.068  1 

6.000  24.278  14.693  .  -.073  1 

9.000  22.324  14.693  -.077  1 

10.000  19.162  14.694  -.061  1 

9.000  16.490  14.693  -.079  1 

7.000  13.975  14.693  -.081  1 

4.000  11.915  14.693  -.083  1 

2.000  10.789  14.693  -.083  1 

.000  10.384  14.693  -.080  1 

-2.000  10.481  14.693  -.083  1 

i 

-4.000  10.769  14.693  -.083  1 

-4.000  11.171  1 4 . 6‘9  2  -.086  1 

-5.000  11.587  14.692  -.084  1 

-6.000  12.029  14.692  -.085  1 

-6.000  12.500  14.692  -.088  1 

-7.000  19.061  14.693  -.077  1 

-6.000  19.338  14.694  -.072  1 

-6.000  19.779  14.693  -.075  1 

-5.000  20.322  14.693  -.073  1 

-4.000  21.200  14.693  -.076  1 

-3.000  21.728  14.693  -.076  1 

-2.000  22.420  14.693  -.079  1 

.000  22.892  14.693  -.077  1 
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9 

<*.633 

2 . 500 
.  240 

2.000 

85.969 

2.000 

22.707 

14 . 693 

-  .  080 

j 

1 

10 

4.635 

2 . 500 
.  262 

2.250 

84.224 

4 .000 

22.530 

14 . 692 

-  .  085 

1 

11 

4.635 

2.500 
.  264 

2 . 500 

84.007 

5.000 

21.010 

14 . 692 

-  .086 

1  ! 

12 

4.636 

2 . 500 
.  271 

2.750 

84.061 

6.000 

19.271 

14 . 693 

-  .079 

i 

1  ! 

13 

4.634 

2.500 
.  251 

3.000 

84.715 

5.000 

16.808 

14 . 692 

-  .088 

i 

i 

14 

4.634 

2.500 
.  246 

3.250 

84.589 

4 . 000 

15.321 

14.692 

-  .095 

i 

15 

4.631 

2 . 500 
.219 

3 . 500 

86.336 

2.000 

13.493 

14 .692 

-.094 

i 

16 

4 . 630 

2.500 
.  205 

3.750 

87.139 

1.000 

12.460 

14.691 

-  .096 

i 

17 

4.628 

2 . 500 
.182 

4 .000 

88.613 

.000 

12.040 

14.692 

-  .095 

1 

18 

4 . 626 

2.500 
.  164 

4.250 

89.726 

-2.000 

11 . 616 

14 . 692 

-  .094 

i 

19 

4 . 626 

2 . 500 
.  160 

4 . 500 

90.245 

-3 . 000 

11.725 

14.692 

-  .090 

i 

20 

4.625 

2 . 500 
.  145 

4.750 

90.971 

-4 . 000 

12.015 

14 . 692 

-  .  093 

i 

21 

4.624 

2  .  500 
.135 

5 . 000 

91.511 

-5.000 

12.219 

14 . 692 

-  .094 

i 

22 

4 . 624 

2 . 500 
.  137 

5 . 250 

91 . 362 

-6 . 000 

12.652 

14 . 692 

-.095 

1 

i 

23 

4 . 624 

2.500 
.  136 

5 . 500 

.  91.328 

-6 . 000 

12.983 

14.691 

-  .096 

i 

1 

4.625 

2.750 

.155 

.  000 

91.392 

-  7 . 000 

19.116 

14.693 

-  .  077 

i 

2 

4.626 

2.750 
.  158 

.  250 

91 . 334 

-6.000 

19.617 

14 . 693 

-  .075 
i 

i 

3 

4.626 

2.750 

.163 

.  500 

91.078 

-6.000 

19.948 

14 . 69,3 

-  .073 

i 

4 

4.628 

2.750 
.  185 

.  750 

89.948 

-4 . 000 

20.765 

14.694 

-  .070 

i 

5 

4  .  628 

2.750 
.  181 

1 . 000 

90.279 

-  3 . 000 

20.865 

14 . 694 

-  .069 

i 

6 

4 . 629 

2  .  750 
.  194 

1.250 

89.358 

-2.000 

21 . 300 

14.694 

-  .071 

i 

7 

4.630 

2  750 
.  209 

1  .  500 

88.115 

-2 . 000 

21.736 

14.693 

-  .  077 

i 

8 

4.631 

2.750 
.  219 

1.750 

87  .  428 

-  1 . 000 

21  .  959 

14.693 

-  .077 

i 

9 

4.633 

2.750 

.239 

2 . 000 

86.165 

1 . 000 

21.878 

14.693 

-  .  078 

i 

10 

4 . 634 

2.750 
.  248 

2.250 

85.369 

4 . 000 

21  .  404 

14.693 

-  .081 

i 

1  1 
4.634 

2.750 
.  253 

2.500 

85.128 

5 . 000 

20.256 

14 . 693 

-  .  080 

i 

12 

4 . 634 

2.750 

.255 

2.750 

85.040 

5 . 000 

18.772 

14 . 693 

-  .079 

i 

13 

4.633 

2.750 

.242 

3 . 000 

85.167 

3 . 000 

17.104 

14 . 692 

-  .090 

i 

14 

4.633 

2.750 

.236 

3.250 

85.256 

2 . 000 

15 . 700 

14 . 692 

-  .095 

i 

15 

4.631 

L 

2.750 

.216 

3  .  500 

86.455 

2.000 

108 

14.368 

14.691 

-  .096 

i 

1  t 

4.630 

2  .  /  50 
.  208 

3.7  50 

U  7  .  2  4  0 

.  uuu 

1 J . 3  OU 

1  4  .  b  9  2 

-  .  U 1 

1 

17 

<4.628 

2.760 
.  180 

4 . 000 

88 .755 

-1.000 

12.647 

14.692 

-  .095 

1 

18 

4.627 

2.750 

.175 

4.250 

88  .  990 

-2.000 

12  573 

14.691 

-  .  096 

1 

19 

4.626 

2.750 
.  161 

4 . 500 

89.990 

-3 . 000 

12.328 

14 . 692 

-  .  094 

1 

20 

4 .625 

2.750 

.152 

4.750 

90.463 

-4 . 000 

12.409 

14 . 592 

-  .094 

1 

21 

4.625 

2.750 
.  147 

5.000 

90. 749 

-  5.000 

12.439 

14.692 

-  .  095 

1 

22 

4 . 624 

2  .  750 
.136 

5 .250 

91 . 233 

-5.000 

12.568 

14.691 

-.098 

1 

23 

4.624 

2.  750 
.  142 

5 . 500 

90.735 

-5.000 

13.073 

14.691 

-  .  101 

1 

1 

4 . 626 

3.000 
.  159 

.000 

91 . 045 

-6 . 000 

19.015 

14.693 

-.078 

1 

2 

4.627 

3 . 000 
.169 

.250 

90.844 

-5 . 000 

19.513 

14.694 

-  .071 

1 

3 

4.627 

3.000 

.177 

.  500 

90.071 

-5 . 000 

20.110 

14.693 

-  .077 

1 

4 

4.629 

3 . 000 
.191 

.750 

89.285 

-4 . 000 

20.535 

14.693 

-  .075 

1 

5 

4.629 

3.000 
.  194 

1 . 000 

89.378 

-3 .000 

20.613 

14 . 694 

-  .071 

1 

6 

4.630 

3.000 
.  201 

1.250 

88.641 

-2.000 

20.975 

14 .693 

-  -  .075 

1 

7 

4.631 

3.000 

.215 

1  .  500 

87 . 780 

-  1 . 000 

21.311 

14 . 693 

-.076 

1 

8 

4.632 

3 . 000 
.  227 

1.750 

86.667 

.  000 

21.519 

14.693 

-  .081 

1 

9 

4.633 

3 . 000 
.  238 

2 . 000 

86. 102 

1 . 000 

21.138 

14. 693 

-  .  079 

1 

10 

4.634 

3 . 000 
.  250 

2.250 

85.416 

3 . 000 

20.609 

14 . 693 

-  .07.9 

1 

1  1 

4.634 

3.000 
.  249 

2  .  500 

85.428 

2.000 

19 . 804 

14 . 693 

< 

-  .079 

1 

12 

4.634 

3.000 

.250 

2.750 

85.133 

3 . 000 

18.556 

14 . 693 

-  .  083 

1 

1  3 
4.633 

3 . 000 
.242 

3.000 

85.059 

2.000 

17.413 

14.692 

-  .0  32 

1 

14 

4.632 

3 . 000 
.  229 

3.250 

85.748 

1 .000 

16.023 

14 . 692 

-  .  094 

1 

15 

4.631 

3 . 000 
.221 

3  .  500 

85 . 929 

1 . 000 

15.21-9 

14.691 

-  .  099 

1 

16 

4.630 

3 . 000 
.  205 

3  .  750 

87 . 250 

.  000 

14 . 095 

14 . 69? 

-  .094 

1 

17 

4.629 

3 . 000 
.  193 

4 . 000 

88.061 

.  000 

13 . 537 

14 . 692 

-  .  093 

1 

18 

4.628 

3 . 000 
.183 

4.250 

88.435 

-2.000 

13.227 

14.69) 

-  .  097 

1 

19 

4.627 

3.000 

.173 

4  ,  500 

89.064 

-2 . 000 

13.298 

14.691 

-  .097 

1 

20 

4.627 

3 . 000 
.171 

4.750 

88.967 

-3 . 000 

13.290 

14.691 

-  .  101 

1 

21 

4.626 

3.000 
.  164 

5.000 

89.440 

-3.000 

13.249 

14.691 

-  .099 

1 

22 

4.626 

3.000 
.  157 

5 .250 

69.756 

-4 . 000 

13.409 

14.691 

-  .  101 

1 

23 

4.626 

3 . 000 
.157 

5 . 500 

89.789 

-4.000 

13.653 

14.691 

-  .  101 

1 
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1 

4.585 

.  000 
-  .  225 

.  000 

135 . 701 

■12. 000 

14.346 

14.737 

.  063 

1 

2 

4.581 

.000 
-  .  262 

.  250 

140 .986 

-13.000 

12.734 

14 . 745 

.  128 

1 

3 

4.576 

.000 
-  .  301 

.  500 

143.144 

-14.000 

12.003 

14.745 

.  132 

1 

4 

(t.576 

.000 
-  .  304 

.750 

143.325 

-15.000 

10.993 

14.745 

.132 

5 

(t  .  570 

.  000 
-  .  351 

1 . 000 

145.356 

17.000 

9 . 953 

14 .745 

.  127 

1 

6 

4.554 

.  000 
-  .  488 

1  .  250 

143.712 

-  15 . 000 

4.871 

14.725 

-  .  043 

1 

7 

4.529 

.000 
-  .  699 

1 . 500 

146 .850 

4 . 000 

-  11 . 861 

14.715 

-  .  128 

1 

8 

4.501 

.  000 
-  .  942 

1.750 

1 5 J . 629 

15. 000 

-  1 9.567 

14.695 

-  .  291 

1 

9 

4.476 

.  000 
-1.148 

2 . 000 

151 . 278 

44 . 000 

-29.743 

14.665 

-  .  547 

1 

10 

4.477 

.  000 
-1.145 

2 . 250 

162.101 

70 . 000 

-12.751 

14.694 

-  .  307 

1 

1  1 

4  .  481 

.  000 
-1.112 

2  .  500 

167 .065 

77 . 000 

2 . 944 

14.713 

-  .  160 

1 

1  2 
4.445 

.  000 
-1.418 

2.750 

175.726 

73 . 000 

9.831 

14.699 

-  .  258 

1 

1  3 
4.411 

.  000 
-1.694 

3 . 000 

172.112 

55.000 

16 . 570 

14 . 656 

-  .  627 

1 

1  4 

4.423 

.000 

-1.597 

3.250 

168 . 349 

26.000 

21.499 

14.657 

-  .  615 

1 

15 

4.468 

.  000 
-1.720 

3  .  500 

168 . 201 

8 . 000 

21.068 

14.701 

-  .  241 

1 

16 

4.491 

.  000 
-1.026 

3.750 

149 . 545 

-6.000 

8 .895 

14.675 

-  .  462 

1 

1  7 

4.575 

.  000 
-  .  736 

4.000 

157 .828 

-  17 . 000 

8.079 

14.731 

.007 

1 

18 

4.645 

.  000 
-  .  ^ 

4  .  250 

154.562 

- 16.0C0 

12 . 284 

14 . 742 

.  107 

1 

1  9 
4.558 

.  000 
-.452 

4  .  500 

149. 529 

-14.000 

12.413 

14.743 

.112 

1 

20 

4.561 

.  000 
-  .  427 

4.750 

148 . 561 

-  14 . 000 

12 . 446 

14.744 

.  116 

1 

21 

4.567 

.  000 

-  .  384 

5 . 000 

146 . 364 

-12.000 

12.679 

14.743 

.  114 

1 

22 

4.572 

.  000 
-  .  336 

5 . 250 

143.843 

-  11 . 000 

12.963 

14.743 

.113 

1 

23 

4.574 

.  000 
-  .  321 

t  ton 

143.081 

-11 . 000 

13 . 261 

14.743 

•  111 

1 

1 

4  .  585 

.250 
-  .  228 

.000 

132.867 

-12.000 

15.605 

14.731 

.007 

1 

2 

4.576 

.  250 
-  .  305 

.  250 

142.775 

-13.000 

13.341 

14.744 

.  121 

1 

3 

4.570 

.  250 
-  .  355 

.  500 

145.758 

-15.000 

12.454 

14.745 

.131 

1 

4 

4.562 

.  250 
-  .  424 

.  750 

148 . 977 

-16.000 

11 . 239 

14.745 

.  128 

1 

5 

4.557 

.  250 
-  .  462 

1 . 000 

150.799 

-  18 . 000 

9.884 

14.745 

.128 

1 
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6 

4.537 

.2  50 
-  .  633 

1.250 

154.569 

-16. 000 

6.363 

14.734 

.038 

1 

7 

4.502 

.250 
-  .  930 

1  .  500 

162.436 

-  2 . 000 

-5.137 

14.7  20 

-  .  084 

1 

6 

4 . 456 

.250 

-1.321 

1.750 

162 . 148 

12.000 

-16. 843 

14.673 

-  .  483 

1 

9 

4 .404 

.  250 
-1.760 

2 . 000 

162 . 590 

42 . 000 

-27 . 952 

14 . 622 

-.911 

1 

10 

4.443 

.  250 
-1 .435 

2 . 250 

164.872 

74.000 

-13.870 

14 . 667 

-  .  533 

1 

1) 

4.458 

.260 

-1.306 

2 . 500 

173.461 

02 . 000 

.  397 

14 . 706 

-  .  201 

1 

12 

4.432 

.  2  50 
-1.521 

2.750 

175 .977 

76 . 000 

2 . 980 

14 . 688 

-  .  355 

1 

13 

4 . 409 

.  250 
-1.722 

3.000 

149 . 982 

60 . 000 

2.154 

14 . 594 

-1  .  149 

1 

34 

4.377 

.  250 
-1.989 

3 . 250 

168 .375 

19.000 

7.057 

14.611  ■ 

-1.011 

1 

15 

4.427 

.  250 
-1.565 

3.500 

175.525 

3.000 

10.704 

14.601 

-  .410 

1 

16 

4.477 

.  250 
-1.143 

3.750 

167 . 877 

-10.000 

6 . 064 

14 . 709 

-  .  172 

1 

17 

4.523 

.  250 
-  .751 

4 . 000 

159.261 

-17 . 000 

9.054 

14.732 

.023 

1 

18 

4.546 

.  250 
-  .  556 

4 .250 

154 . 163 

-  16 . 000 

11 . 100 

14.742 

.  106 

1 

19 

4.563 

.250 
-  .  428 

4 . 500 

148 . 675 

-15 . 000 

11.766 

14.744 

.118 

1 

20 

4 . 568 

.250 
-  .  373 

4.750 

346.031 

-14 . 000 

12 . 087 

14 . 744 

.118 

1 

21 

4.568 

.250 
-  .  368 

5 . 000 

145.914 

-  12 . 000 

12.199 

14.744 

.121 

1 

22 

4.573 

.  250 
-  .  330 

5.250 

143.882 

-  11 . 000 

12.613 

14 . 744 

.118 

1 

23 

4.576 

.250 
-  .  302 

5  .  500 

142.481 

-10.000 

12.946 

14.744 

.  118 

1 

1 

4.582 

.  500 
-  .  252 

.000 

139.198 

-12.000 

14 .892 

14  .  7%2 

.  103 

1 

2 

4.577 

.  500 
-  .  297 

.250 

142 . 098 

-13.000 

13.863 

14.743 

.115 

1 

3 

4  .  569 

.  500 
-  .  361 

.  500 

145.454 

-14.000 

13 . 243 

14.744 

.119 

1 

4 

4  .  554 

.  500 
-  .  491 

.750 

151.754 

-16.000 

12.04.3 

14.744 

.  120 

1 

5 

4.542 

.  500 
-  .  595 

1.000 

156. 454 

-IB . 000 

10.963 

14 . 744 

.117 

1 

6 

4  .  524 

.  500 
-  .  740 

1.250 

160.952 

-  18.000 

7.847 

14 . 73B 

.072 

1 

7 

4.484 

.  500 
-1.085 

1 . 500 

168 . 533 

-9.000 

1  .  301 

14.718 

-  .098 

1 

8 

4.447 

.  500 
-1.395 

1.750 

163.098 

5.000 

-3.047 

14 . 667 

-  .534 

1 

9 

4  .  396 

.  500 
-1.826 

2.000 

143.480 

35. 000 

-26.721 

14 . 566 

-1.387 

1 

10 

4.426 

.  500 
-1.575 

2 .250 

171.707 

72.000 

-4 . 330 

14.669 

-  .513 

3 

1  1 

4.447 

.  500 
-1.397 

2.500 

182 . 794 

80.000 

3.619 

14.723 

-  .060 

1 

17 

4.429 

.  500 
-1.551 

2.750 

183.419 

78 . 000 

-2.722 

14 . 706 

-.198 

1 

Ill 


20 

.750 

4.575 

-  .  312 

21 

.750 

4.575 

-  .  313 

22 

.750 

4.576 

-.306 

23 

.  750 

4.578 

-  .  286 

1 

1.000 

4 . 580 

-  .  274 

2 

1.000 

4.578 

-  .  280 

3 

1.000 

4 . 568 

-  .  3  /3 

4 

1 . 000 

4 . 557 

-  .465 

5 

1 . 000 

4 . 542 

-  .  590 

6 

1.000 

4.522 

-  .759 

7 

1 . 000 

4.494 

-1.000 

8 

1 .000 

4.473 

-1.175 

9 

1 . 000 

4.467 
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2 . 500 
.023 

2 . 000 

124 . 946 

3.000 

26.810 

14.743 

.115 

1 

10 

<1.621 

2 . 500 
.077 

2.250 

121.848 

0 .000 

26.05/ 

14.743 

.115 

1 

1  1 
<*.627 

2.500 
.  125 

2.500 

118. 980 

12.000 

23.961 

14.743 

.  115 

1 

12 

<*  .  628 

2.500 

.138 

2.750 

118 . 270 

13.000 

20.871 

14 . 744 

.116 

1 

13 

<*.625 

2.500 
.  1 14 

3.000 

119.701 

12.000 

16.567 

14 . 744 

.116 

1 

1<* 

<* .  621 

2.500 
.  079 

3.250 

121 . 593 

10.000 

12.813 

14.743 

.113 

1 

15 
<* .  615 

2.500 
.  027 

3 . 500 

124 .659 

5 . 000 

10.150 

14.743 

.  114 

1 

16 

<*.611 

2 . 500 
-  .007 

3.760 

126 . 673 

2 . 000 

8 . 682 

14.743 

.  115 

1 

17 

<* .  60<* 

2.500 
-  .067 

4 . 000 

130.052 

-1 . 000 

8 . 662 

14.744 

.  116 

1 

10 

<*  .  600 

2.500 
-  .  104 

4 . 250 

132.067 

-3.000 

9.495 

14.743 

.116 

1 

19 
<* .  696 

2.500 
-  .  136 

4.500 

133.770 

-4.000 

10.173 

14 . 743 

.116 

1 

20 
<*  .  S9<* 

2 . 500 
-  .  155 

4.750 

134.656 

-5 . 000 

10. 938 

14.743 

.113 

1 

21 

4 . 593 

2 . 500 
-  .  162 

5.000 

134.945 

-6 . 000 

11 . 654 

14.743 

.  112 

1 

22 

4.591 

2.500 
-  .  176 

5 . 250 

135.749 

-7 . 000 

12.251 

14.743 

.113 

1 

23 

4 . 591 

2 . 500 
-  .  179 

5 . 500 

135 . 960 

-7 . 000 

12  .  703 

14.743 

.114 

1 

1 

4 . 594 

2.750 
-  .  151 

.  000 

134 . 168 

-8 . 000 

18 . 668 

14.743 

.  108 

1 

2 

4.592 

2.750 
-  .  169 

.  250 

135.286 

-8.000 

18.770 

14 . 743 

■  111 

i 

1 

3 

4. 596 

2.750 
-  .  132 

.  500 

133. 362 

-7 . 000 

19.800 

14. 743 

.112 

1 

4 

4 . 596 

2.7  50 
-  .  131 

.750 

133.330 

-7.000 

20.206 

14 . 7*43 

.112 

1 

5 

4 . 599 

2.750 

-.111 

1.000 

132 . 317 

-6.000 

21.125 

14 . 743 

.113 

1 

6 

4.606 

2 . 750 
-  .048 

1 . 250 

128.822 

-4.000 

22.820 

14.743 

.113 

1 

7 

4 . 609 

2 . 750 
-  .026 

1 . 500 

127.548 

-2.000 

23 . 608 

14.743 

.112 

1 

8 

4.611 

2.750 
-  .008 

1.750 

126.496 

.000 

24 . 278 

14.743 

•  111 

1 

9 

4.614 

2.750 

.022 

2.000 

124.958 

3 . 000 

24 . 495 

14 . 743 

.114 

1 

10 

4.621 

2.750 
.  074 

2 . 250 

121.745 

7 . 000 

24 . 064 

14.743 

•  111 

1 

11 

4.624 

2.750 
.  103 

2 . 500 

120.189 

10.000 

22.501 

14 .743 

.114 

1 

12 

4.625 

2.750 

.110 

2.750 

119.673 

10 . 000 

20.362 

14 . 743 

.112 

1 

13 

4.624 

2.750 
.  103 

3.000 

120.237 

9.000 

17.457 

14.743 

.114 

1 

14 

4.619 

2.750 

.062 

3.250 

122.532 

7 . 000 

14 . 598 

14.743 

.112 

1 

15 

4.614 

2 .750 
.022 

3 . 500 

124.683 

4 . 000 

12 . 629 

14.743 

* 

.  109 

1 

118 


1  6 

4 . 609 

2.7  50 
-  .025 

3.750 

127 . 432 

1 . 000 

13.250 

14.743 

■  111 

1 

17 

4 . 605 

2 .750 
-  .059 

4 . 000 

129.310 

-  1.000 

10.685 

14.743 

•  111 

1 

18 

4.602 

2.750 
-  .082 

4.250 

130 . 544 

-2 . 000 

10 . 684 

14.743 

.  109 

1 

19 

4 . 599 

2.750 
-  .  109 

4 . 500 

132.113 

-3 . 000 

11.232 

14.743 

.  112 

1 

20 

4 . 594 

2.750 
-  .  149 

4.750 

134 . 199 

-4.000 

11.598 

14.743 

•  111 

1 

21 

4.594 

2.750 
-  .  150 

5 . 000 

134.250 

-5.000 

12.084 

14 . 743 

•  111 

1 

22 

4 . 593 

2 .750 
-  .162 

5 .250 

134.819 

-6.000 

12.688 

14.743 

.  109 

1 

23 

4.592 

2.750 
-  .  171 

5 . 500 

135.451 

-6.000 

12.928 

14.743 

.  113 

1 

1 

4 . 595 

3.000 

-.140 

.000 

133.856 

-6.000 

18 . 782 

14.743 

.114 

1 

2 

4 . 599 

3 . 000 
-  .  107 

.250 

131 .943 

-5 . 000 

19.493 

14 .743 

.111 

1 

3 

4.598 

3.000 
-  .  114 

.  500 

132.324 

-5 . 000 

19.619 

14.743 

■  111 

1 

4 

4.599 

3 . 000 
-  .  112 

.  750 

132.137 

-6.000 

20.044 

14.743 

.  109 

1 

5 

4 . 603 

3 . 000 
-  .075 

1.000 

130. 158 

-4.000 

21.078 

14 . 74  3 

.110 

1 

6 

4 . 605 

3.000 
-  .  055 

1 . 250 

128 .937 

-3.000 

21  .836 

14 . 74? 

.  107 

1 

7 

4.611 

3 . 000 
-  .  010 

1 . 500 

126. 693 

-1 .000 

22.696 

14.743 

.113 

1 

8 

4.613 

3 . 000 
.007 

1.750 

125.627 

1 . 000 

23.022 

14.743 

.  Ill 

1 

9 

4.617 

3.000 
.  046 

2.000 

123 . 270 

3.000 

23.437 

14.743 

.  108 

i 

1 

10 

4 .620 

3.000 

.069 

2 .250 

121.979 

5.000 

22 . 902 

14.743 

.  109 

1 

11 

4 . 623 

3.000 
.  0q5 

2 . 500 

120.487 

7.000 

21.763 

14.743 

1 

•  111 

1 

12 

4.623 

3.000 

.091 

2.750 

120. 773 

8.000 

19.838 

14.743 

.111 

1 

13 

4 . 620 

3.000 

.067 

3.000 

122 . 215 

6.000 

17.752 

14.743 

.112 

1 

14 

4.615 

3.000 

.030 

3.250 

124 . 307 

4 . 000 

15.308 

14.743 

.113 

1 

15 

4.613 

3.000 

.014 

3  .  500 

125.427 

2.000 

13.911 

14 .743 

.114 

1 

16 

4.609 

3 . 000 
-  .020 

3.750 

127.281 

1.000 

12.724 

14.743 

.113 

1 

17 

4 . 606 

3 . 000 
-  .040 

4 . 000 

328.661 

.  000 

32 . 058 

14.743 

.  109 

1 

18 

4.602 

3  .  UOU 
-  .082 

4.250 

130.582 

-  2 . 000 

11 . 948 

14.743 

•  111 

1 

19 

4 . 598 

3 . 000 
-  .  116 

4  .  500 

132.377 

-  3 . 000 

12.096 

14.743 

.110 

1 

20 

4 . 598 

3 . 000 
-  .  121 

4.750 

132.726 

-4 . 000 

12.532 

14 . 743 

•  111 

1 

21 

4 . 595 

3 . 000 
-  .  139 

5 . 000 

133.751 

-5.000 

12.591 

14.743 

.112 

1 

22 

4 . 595 

3 . 000 
-.140 

5 . 250 

133.585 

-5.000 

13.246 

14 . 743 

.  108 

1 

23 

4 . 595 

3 . 000 
-  .  145 

5  .  500 

133 . 965 

-6.000 

13.427 

.14.743 

.  110 

1 

119 
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1 

4.614 

.000 
.  303 

.000 

96.737 

-11.000 

13.372 

14 . 690 

.  172 

1 

2 

4  .  612 

.000 
.  279 

.250 

97.806 

-12.000 

12.936 

14.690 

.  168 

1 

3 

4.610 

.  000 
.  259 

.500 

99.054 

-13.000 

12.116 

14 . 690 

.  170 

1 

4 

4 . 609 

.000 
.  243 

.  750 

99.880 

-14.000 

11.094 

14 . 690 

.  169 

1 

5 

4 . 607 

.000 
.  225 

1.000 

100 .846 

-16.000 

10.170 

14.690 

.  169 

1 

6 

4.599 

.000 
.  133 

1.250 

100.238 

-15.000 

8 . 460 

14.681 

.065 

1 

7 

4.573 

.  000 
-  .  160 

1.500 

98.974 

-3 . 000 

-4 . 648 

14.653 

-  .  250 

1 

8 

4.553 

.000 
-  .  386 

1.750 

98. 912 

15.000 

-13.092 

14.633 

-.478 

1 

9 

4.528 

.  000 
-  .  673 

2.000 

94 . 907 

37 . 000 

-22.132 

14 . 602 

-  .  837 

1 

10 

4 . 538 

.000 
-  .  565 

2  .  250 

99.870 

65 . 000 

-13.838 

14.619 

-  .  639 

1 

11 

4 . 550 

.  000 
-  .  426 

2 . 500 

104 . 963 

73 . 000 

2  .  199 

14 . 640 

-  .403 

1 

12 

4.539 

.  000 
-  .  547 

2.750 

108 .612 

71.000 

10. 560 

14. 636 

' - . 452 

1 

13 

4.513 

.000 
-  .  840 

3 . 000 

103 . 958 

59.000 

15.672 

14 . 602 

-  .837 

1 

14 

4.  509 

.000 
-  .  888 

3 . 250 

102 .772 

30.000 

19.279 

14 . 595 

-  .  907 

1 

15 

4.538 

.  000 
-  .  561 

3 . 500 

103 .466 

11.000 

21.003 

14.625 

-  .  567 

i 

1 

16 

4.551 

.000 
-  .  417 

3.750 

101.685 

-4 .000 

12.623 

14  .,635 

-  .457 

1 

17 

4.582 

.000 
-  .064 

4.000 

107 . 602 

-15.000 

10.328 

14.676 

.011 

1 

18 

4.593 

.000 
.  070 

4 . 250 

104 . 468 

-16.000 

13.372 

14.683 

.083 

1 

19 

4 . 600 

.000 
.  144 

4 . 500 

103.733 

-14.000 

12.835 

14.688 

.143 

1 

20 

4 . 603 

.000 
.  180 

4 .750 

101 . 574 

-13.000 

12.751 

14.687 

.138 

1 

21 

4 . 606 

.  000 
.210 

5 . 000 

99.819 

-12.000 

12.832 

14.687 

.135 

1 

22 

4 . 607 

.000 
.  226 

5.250 

98 .659 

-11.000 

12.899 

14.687 

.130 

1 

23 

4 . 609 

.000 
.  252 

5 . 500 

96 . 949 

-10.000 

13.331 

14.686 

.  125 

1 

1 

4.612 

.  250 
.  284 

.000 

96 . 502 

-  11 . 000 

14.210 

14.688 

.  148 

1 

2 

4.610 

.  250 
.  261 

.250 

98 . 184 

-  13 . 000 

13 . 427 

14 . 689 

.157 

1 

3 

4 . 608 

.250 
.  240 

.  500 

99 . 601 

-14.000 

32.825 

14 . 690 

.  161 

1 

4 

4 . 605 

.  250 
.  201 

.  750 

102.050 

-16.000 

11.671 

14 . 690 

.  168 

1 

5 

4.601 

.  250 
.  160 

1.000 

103.829 

-17.000 

10.458 

*  14 . 689 

.161 

1 

120 


6 

4.593 

.  250 
.063 

1 . 250 

104 . 395 

-  17 . 000 

8.573 

14 . 682 

.074 

1 

7 

4 . 566 

.250 
-  .  242 

1 . 500 

103 .435 

-  7 . 000 

-  .  823 

14.653 

-  .  249 

1 

B 

4.537 

.250 
-  .  574 

1.750 

103.086 

-  1 . 000 

-7.788 

14.624 

-  .  507 

1 

9 

4 . 506 

.  250 
-  .  923 

2.000 

94.559 

39 . 000 

-22 . 304 

14.579 

-  1 . 093 

1 

10 

4.518 

.250 
-  .786 

2.250 

100.821 

60 . 000 

-13 . 193 

14 . 601 

-  .  842 

1 

11 

4.541 

.  250 
-.524 

2.500 

110.243 

77.000 

.  516 

14.640 

-.396 

1 

12 

4.529 

.250 
-  .  662 

2.750 

113.049 

74 . 000 

5 . 440 

14 . 633 

-  .  475 

1 

13 

4.490 

.  250 
-1.102 

3 . 000 

101 . 824 

61 . 000 

3.994 

14.575 

-1.139 

1 

14 

4.474 

.  250 
-1.285 

3 . 250 

101.485 

26.000 

4.513 

14.558 

-1.329 

1 

15 

4.513 

.250 
-  .843 

3.500 

106.355 

6.000 

11.286 

14.606 

-  .  792 

1 

16 

4.544 

.  250 
-  .  490 

3.750 

107 .919 

-9.000 

8  .  127 

14.639 

-  .408 

1 

17 

4 . 580 

.250 
-  .  086 

4.000 

106 . 988 

-17 .000 

10.581 

14.673 

-  .023 

1 

18 

4 . 596 

.250 
.  101 

4.250 

104 . 264 

-18.000 

12 . 185 

14 . 685 

„  .110 

1 

19 

4 . 602 

.  250 
.  170 

4 . 500 

102 . 121 

-16.000 

12.377 

14 . 687 

.138 

1 

20 

4 . 604 

.  250 
.  188 

4 . 750 

101.333 

-14.000 

12.111 

14 . 6  8  B 

.  141 

1 

21 

4 . 606 

.  250 
.  214 

5.000 

99.732 

-  12 . 000 

12.258 

14 . 687 

.137 

1 

22 

4 . 609 

.  250 
.242 

5.250 

97 . 907 

-11.000 

12.518 

14.607 

.132 

1 

23 

4.611 

.250 

.266 

5  .  500 

96 . 575 

-10.000 

12.897 

14.687 

.132 

1 

1 

4 . 606 

.  500 
.211 

.000 

91.565 

-11.000 

16.192 

14 . $74 

-  .010 

1 

2 

4 . 609 

.  500 
.  244 

.250 

97 . 493 

-12.000 

14.454 

14.686 

.126 

1 

3 

4 . 606 

.  500 
.  211 

.  500 

100. 384 

-14.000 

13  '  25 

14.688 

.  147 

1 

4 

4 . 602 

.  500 
.  163 

.750 

103.031 

-16.000 

12.529 

14.688 

.  148 

1 

5 

4  .  598 

.  500 
.118 

1 .000 

105.289 

-18.000 

11.701 

14 . 688 

.  147 

1 

6 

4 . 588 

.  500 
.010 

1  .  250 

107 . 485 

-18.000 

10.302 

14.683 

.083 

1 

7 

4 . 564 

.  500 
-  .  266 

1  .  500 

108.919 

-10.000 

4.025 

14.661 

-  .  164 

1 

8 

4.531 

.  500 
-  .  645 

1 .750 

99.871 

5.000 

-1 . 323 

14.612 

-  .  719 

1 

9 

4.488 

.  500 
-1.131 

2 . 000 

92.828 

33.000 

-15.353 

14 . 558 

-1.331 

1 

10 

4.514 

.  500 
-  .  835 

2.250 

105.376 

77 . 000 

-8.410 

14 . 605 

-  .804 

1 

11 

4.544 

.  500 
-.493 

2  .  500 

114.442 

78 . 000 

1.004 

14.651 

-  .277 

1 

12 

4.533 

.  500 
-  .  619 

2.750 

115.768 

75 . 000 

.447 

14.642 

-  .  375 

1 

121 


13 

.  500 

3 . 000 

106.614 

61 . 000 

-  8 . 569 

14.579 

-1.088 

1 

4  .  <.87 

-1.144 

14 

.  500 

3.250 

105.033 

23 . 000 

-14.855 

14 . 559 

-1.317 

1 

4.469 

-1.341 

15 

.  500 

3  .  500 

110.667 

1 . 000 

-1.026 

14 . 609 

-  .  753 

1 

4 . 509 

-  .  890 

16 

.  500 

3.750 

111.552 

-12.000 

3 . 680 

14. 649 

-  .  302 

1 

4.547 

-.458 

17 

.  500 

4.000 

107 . 804 

-18.000 

8.075 

14.672 

-.038 

1 

4.577 

-  .  117 

18 

.  500 

4 .250 

105 . 384 

-18.000 

10.495 

14.685 

.110 

1 

4 . 594 

.079 

19 

.  500 

4 . 500 

100 . 400 

-16. 000 

11.451 

14.687 

.  1  30 

1 

4 . 604 

.  194 

20 

.  500 

4.750 

99 . 509 

-14 . 000 

11.722 

14.687 

.127 

1 

4 . 606 

.  207 

21 

.  500 

5 . 000 

98 . 404 

-  12 . 000 

11.801 

14.686  . 

.124 

1 

4 . 607 

.  225 

22 

.  500 

5 . 250 

97.393 

-  11 . 000 

12.154 

14.686 

.121 

1 

4  .  608 

.  241 

23 

.  500 

5 . 500 

96 . 005 

-  10 . 000 

12.782 

14.686 

.  118 

1 

4.610 

.  262 

1 

.  750 

.  000 

97 . 198 

-12.000 

15.176 

14.688 

.145 

1 

4.611 

.  268 

*■* 

.  750 

.  250 

98.757 

-13 . 000 

14.856 

14 . 689 

.151 

1 

4 . 609 

.  245 

3 

.750 

.  500 

100.395 

-  14 . 000 

14.467 

14.689 

.  154 

1 

4.607 

.219 

4 

.  750 

.750 

•103.847 

-  15 . 000 

13. 774 

14.689 

.156 

1 

4.601 

.  1  S  5 

i  5 

.  750 

1.000 

106 . 623 

-16.000 

13.358 

14 . 689 

.151 

1 

4 . 596 

.095 

6 

.  750 

1  .  250 

109 . 165 

-19.000 

12.820 

14.683 

.082 

1 

4 . 585 

-  .  024 

1 

7 

.750 

1  .  500 

109 . 984 

-12.000 

9 . 363 

14 . 663 

-  .  142 

1 

4  .  S  6  4 

-  .  265 

8 

.750 

1.750 

106 .534 

.000 

7 . 824 

14i 630 

-  .  512 

1 

4.537 

-  .  566 

9 

.750 

2.000 

97.337 

26.000 

1 .073 

14.581 

-1.075 

1 

4.503 

-  .  954 

10 

.  750 

2.250 

103 . 220 

61.000 

.  406 

14 . 603 

-  .817 

1 

4.516 

-  .  806 

1  1 

.750 

2  .  500 

117.028 

74 . 000 

3.171 

14.662 

-.148 

1 

4  .  550 

-  .  420 

12 

.  750 

2.750 

117. 409 

74 . 000 

-  .  160 

14.658 

-  .  198 

1 

4.545 

-  .  478 

13 

.  750 

3 . 000 

112.016 

56 . 000 

-15.666 

14.605 

-  .  799 

1 

4.502 

-  .  964 

14 

.750 

3.250 

111  .  196 

23 . 000 

-22 . 328 

14 . 597 

-  .  884 

1 

4*496 

-1.032 

15 

.750 

3 . 500 

115.114 

-  1  .000 

-10.267 

14.626 

-  .  564 

1 

4.517 

-  .  794 

16 

.  750 

3.750 

112 . 534 

-12.000 
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